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b(x)  Dirac  delta  function 

A  dilatation 

V2  Laplacian 

scalar  potential  function 

<f>  vector  of  the  nodal  values  of  the  electric  potential 
W  vector  potential  function 

0  angular  coordinate  in  cylindrical  coordinate  system 

wavenumber  multiplier,  Lame’s  constant 
X  wavelength,  Lame’s  constant 

p  density 

p;  density  of  fluid  in  region  I 

ps  density  of  fluid  in  region  III 

pf  density  of  infinite-length  elastic  tube  in  region  II 

a  coordinates  of  surface  point 

t  stress 

to  angular  frequency  ( 2nf) 

|  integration  variable 
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INVESTIGATION  OF  A  CYLINDRICAL 
NONACOUSTIC- WAVENUMBER 
CALIBRATION  ARRAY 


INTRODUCTION 

This  thesis  describes  an  investigation  of  a  prototype  Cylindrical  Wavenumber  Calibrator  (CWC). 
The  CWC  is  a  cylindrical,  tubular  device  designed  to  generate  nonacoustic-wavenumber  pressure  fields 
in  water  within  the  CWC.  It  can  therefore  be  used  to  simulate  the  effects  of  flow  noise  on  a  device 
placed  in  the  water.  The  research  consists  of  both  theoretical  and  experimental  investigations  of  the 
CWC. 

Wavenumber  Calibration 

There  has  been  considerable  interest  in  recent  years  in  the  acoustic  effects  of  the  turbulent 
boundary  layer  (TBL)  produced  by  the  motion  of  a  vehicle  in  water.  There  is  concern  about 
interference  by  the  pressure  field  in  the  TBL  with  the  proper  operation  of  transducers  used  in  sonar 
arrays.  Sonar  arrays  are  systems  of  transducers,  often  distributed  over  areas  of  a  vessel  or  in  lines 
towed  behind  a  vessel,  designed  to  detect  acoustic-pressure  fluctuations  in  the  water.  There  is  also 
concern  about  the  excitation  of  structures  of  the  vehicle,  which  results  in  undesirable  acoustic  fields 
both  internally  and  externally. 

The  "flow  noise"  produced  by  the  TBL  is  generally  characterized  by  pressure  fluctuations  that  are 
transported  in  the  flow  direction  at  a  speed  somewhat  less  than  the  free-stream  velocity.  Since  these 
fluctuations  are  transported  at  speeds  less  than  the  sound  speed  in  water,  they  have  nonacoustic 
wavenumbers  (i.e.,  wavenumbers  higher  than  acoustic  wavenumbers  at  the  same  frequency  in  water) 
and  are,  therefore,  evanescent. 

Several  methods  have  been  employed  to  reduce  the  undesirable  effects  of  the  flow  noise.  One 
method  used  is  to  slow  the  motion  of  the  vehicle  through  the  water,  thereby  reducing  the  amplitude  of 
the  flow  noise.  This  method  is  effective  but  can  cause  serious  tactical  problems. 

Another  technique  for  reducing  the  undesirable  effects  of  the  flow  noise  is  to  spatially  separate 
the  acoustic  sensors  from  the  TBL.  This  is  sometimes  accomplished  by  using  a  sonar  dome  or  a  layer 
of  acoustically  transparent  material.  Since  the  convective  field  is  evanescent,  a  relatively  small 
separation  distance  can  have  a  significant  effect.  This  technique  has  limitations  at  low  frequencies, 
however,  where  the  beneficial  effect  is  reduced. 

A  third  technique  for  reducing  the  flow  noise  problem  is  spatial  averaging.  This  is  accomplished 
by  using  sonar  sensors  that  are  longer  than  the  flow-noise  wavelength. 

A  major  problem  encountered  in  studying  the  effects  of  flow  noise  on  acoustic  arrays  and 
underwater  structures  is  the  fact  that  the  structure  of  the  TBL  tends  to  be  poorly  defined  and  somewhat 
random  in  nature.  Because  of  this,  it  is  difficult  to  reproduce  the  same  flow  noise  from  one 
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measurement  to  the  next.  To  gain  better  understanding  of  the  effects  of  flow  noise,  it  would  be  useful 
to  be  able  to  plot  or  calibrate  the  response  of  an  array  or  structure  as  a  function  of  the  wavenumber  of 
the  excitation  for  frequencies  of  interest  [1],  This  could  be  achieved  if  a  device  existed  that  could 
generate  a  pressure  field  consisting  of  essentially  a  single  nonacoustic-wavenumber  component.  The 
structure  being  analyzed  could  then  be  placed  in  the  pressure  field  and  its  behavior  monitored  for 
different  wavenumbers. 

This  thesis  concerns  an  investigation  of  a  cylindrical  fluid-filled  array  designed  to  generate  the 
desired  nonacoustic-wavenumber  pressure  fields  within  the  array.  The  work  consisted  of  both 
theoretical  and  experimental  investigations  of  the  array. 

Literature  Review 

Effects  of  Flow  Noise 

Over  the  years  numerous  attempts  have  been  made  to  better  understand  flow  noise  and  its  effects 
on  acoustic  measurements  in  general  [2,3].  There  has  been  specific  interest  in  the  effects  of  flow 
noise  on  sonar  sensors  and  underwater  structures. 

A  characteristic  wavevector-frequency  spectrum  of  TBL  wall  pressure  [4]  is  shown  in  Fig.  1.  In 
this  figure,  kx  is  the  wavenumber  in  the  flow  direction  with  kx  and  k%  lying  in  the  plane  of  the  wall. 

In  general,  a  qualitative  description  of  the  spectral  density  of  the  flow  noise  versus  the  streamwise 
wavenumber  at  a  single  frequency  shows  a  convective  domain  or  "convective  ridge"  located  at  the 
hydrodynamic  wavenumber  kh  =  co/w,  where  a)  is  the  radian  frequency,  and  u  is  the  convection 
velocity,  which  is  a  large  fraction  of  the  free-stream  velocity.  The  spectral  density  decreases  with 
decreasing  wavenumber  below  the  convective  ridge.  The  region  below  the  convective  ridge  and  above 
the  acoustic  wavenumber  ka  =  co/c,  where  c  is  the  sound  speed  in  the  fluid,  is  sometimes  called  the 
low-wavenumber  region  of  the  flow  noise, 

Corcos  [5]  discussed  the  ability  of  a  transducer  to  resolve  the  spatial  details  of  a  turbulent 
pressure  field  with  two  opposite  applications  in  mind.  One  involved  the  use  of  a  transducer  with  a 
sensing  area  small  enough  to  give  a  good  account  of  the  fine  spatial  structure  of  the  turbulent  field  for 
studying  the  field.  The  other  involved  the  use  of  transducers  with  sensing  areas  large  enough  to  be 
insensitive  to  the  flow  noise  within  a  given  frequency  band  (i.e.,  larger  than  the  flow-noise 
wavelength)  and  yet  remain  sensitive  to  incoming  sound  waves  at  the  same  frequency.  The  latter  case 
is  the  one  generally  desired  for  sonar  sensors.  Corcos  proposed  a  theoretical  model  that  could  be  used 
to  determine  the  response  of  a  transducer  to  a  turbulent  pressure  field  using  a  size-correction  factor 
and  gave  asymptotic  formulas  for  large  transducers  that  yield  estimates  of  the  degree  to  which  a  flush- 
mounted  sonar  receiver  is  able  to  discriminate  against  flow  noise.  Gilchrist  and  Strawderman  [6] 
presented  experimental  results  that  showed  good  agreement  with  Corcos’  predictions  for  the  size 
correction  value. 

Heaps  [7]  extended  Corcos’  analysis  to  treat  the  case  where  the  transducer  is  separated  from  the 
turbulent  layer  by  a  medium  free  of  turbulent  sources  and  through  which  pressure  fluctuations  are 
transmitted  acoustically.  A  theory  was  developed  to  compute  the  attenuation  of  spectral  density  as  a 
function  of  transducer  size  and  distance  from  the  turbulent  flow.  The  theory  predicts  that  (for  small 
Mach  numbers)  the  attenuation  is  effected  more  by  a  change  in  the  separation  distance  than  by  a 
similar  change  in  the  radius  of  a  circular  transducer. 

There  has  also  been  considerable  interest  in  recent  years  in  the  effects  of  flow  noise  on 
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structures.  If  structures  on  a  vehicle  moving  through  the  water  are  excited  by  the  flow  noise,  they  can 
then  radiate  into  the  vehicle  as  well  as  into  the  surrounding  medium.  The  evanescent  energy  in  the 
flow  noise  can,  therefore,  be  converted  into  acoustic  energy  that  can  propagate  over  large  distances. 
Chandiramani  [8]  presented  a  method  for  estimating  the  vibrational  response  of  fluid-loaded  structures 
when  excited  by  low-speed  TBL  flow.  An  approximate  method  was  also  suggested  for  estimating 
vibrational  characteristics  of  spatially  nonhomogeneous  structures  excited  by  surface  interactions. 

Chandiramani  [9]  later  modified  an  exact  theory  [10]  of  transmission  of  sound  across  an  isotropic 
flat  plate,  involving  longitudinal  and  shear  waves,  to  be  applicable  to  convective  wavenumbers  and 
presented  closed-form  results. 
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Wavenumber  Filters 

One  problem  encountered  in  measuring  the  wavenumber  content  of  the  low-wavenumber  region 
of  the  flow  noise  is  the  fact  that  the  typically  low  levels  of  the  wavenumber-frequency  spectral  density 
are  easily  contaminated  by  the  background  acoustic  noise.  Maidanik  and  Jorgensen  [11]  proposed  a 
solution  to  this  problem  in  the  form  of  a  wave-vector  filter.  The  filter  would  consist  of  an  array  of 
equally-spaced  transducers  aligned  parallel  to  the  flow  direction.  It  was  shown  that  with  proper 
choices  of  sensor  size,  spacing,  and  polarity,  a  filter  could  be  constructed  to  pass  a  single  narrow  band 
of  wavenumbers.  A  particular  example  given  shows  14  dB  of  attenuation  outside  the  passband. 

Blake  and  Chase  [12,13]  constructed  a  four-element  microphone  array  and  tested  it  as  a  wave- 
vector  filter  in  a  wind  tunnel.  Unfortunately,  they  were  unable  to  conclusively  establish  whether  the 
dominant  array  noise  in  the  passband  was  spurious  acoustic  noise  or  flow  noise.  Therefore,  the 
measurements  placed  only  an  upper  limit  on  the  flow-noise  contribution. 

It  should  be  noted  that,  although  reference  is  made  in  this  thesis  to  both  wavenumber  filters 
(designed  primarily  as  receiving  devices)  and  wavenumber  calibrators  (designed  primarily  as 
transmitting  devices),  if  a  device  is  reciprocal,  it  can  be  used  for  either  purpose. 

Previous  Wavenumber-Calibration  Devices 

Several  techniques  have  been  considered  previously  for  wavenumber  calibration  with  planar 
geometries.  Powers  and  Sherman  [14]  examined  three  techniques  for  hydrophones;  beam  patterns 
taken  in  air  (rather  than  water),  the  use  of  an  evanescent-wave  generator,  and  measurement  of  the 
spatial  impulse  response.  Devices  were  constructed  and  measurements  made  to  investigate  the  first 
two  techniques. 

Since  the  sound  speed  is  approximately  4.4  times  slower  in  air  than  in  water,  measuring  the  beam 
pattern  of  a  hydrophone  in  air  instead  of  in  water  simulates  wavenumbers  that  are  approximately  4.4 
times  higher  than  in  water.  Powers  and  Sherman  constructed  a  12.7x12. 7-cm  hydrophone  consisting 
of  15  piezoelectric  ceramic  cylinders  potted  in  urethane.  Beam  patterns  taken  at  20  kHz  in  water  and 
5  kHz  in  air  were  compared.  The  beam  patterns  were  found  to  be  significantly  different  even  though 
these  two  cases  had  similar  wavenumbers.  A  probable  reason  given  for  the  differences  is  the  different 
scattering  from  the  hydrophone  in  the  two  cases.  The  acoustic  impedance  of  the  urethane  potting 
material  is  a  good  match  to  water  but  a  poor  match  to  air.  This  would  cause  much  more  scattering 
from  the  potting  in  the  air  case  than  in  the  water  case.  This  difference  in  scattering  is  a  significant 
problem  when  trying  to  simulate  water  measurements  in  air. 

Powers  and  Sherman  also  constructed  and  tested  an  evanescent-wave  generator  in  air.  The  planar 
array  was  10-cm  long  and  7-cm  wide  and  contained  22  strip  elements.  Each  element  consisted  of  an 
air  backed  hemi-cylindrical  membrane  of  piezoelectric  polymer.  The  elements  were  wired  so  that  the 
drive  at  each  element  could  be  inverted  (i.e.,  shifted  by  180°).  The  array  was  driven  in  air  at  3  and  6 
kHz.  By  changing  the  frequency  and  phasing  of  the  elements,  wavenumbers  that  were  7.9,  8.8,  and 
16.4  times  the  acoustic  wavenumber  in  air  were  generated.  Plots  of  the  pressure  in  front  of  the  array 
were  obtained  using  a  small  probe  microphone.  A  spatial  Fourier  transform  was  then  performed  on 
the  pressure  data  to  determine  the  wavenumber  content.  The  results  showed  that  high  wavenumbers 
could  be  generated  independent  of  the  drive  frequency  using  a  phased  transducer  array.  It  was  also 
demonstrated  that  the  amplitude  of  these  waves  decreased  exponentially  with  increasing  normal 
distance  from  the  array  as  would  be  expected  for  evanescent  waves.  It  was  further  shown  that 
undesired  acoustic  wavenumbers  were  also  generated  that  contaminated  the  measured  field  to  a  level 
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that  would  make  the  array  unusable  for  wavenumber  calibration.  It  is  concluded  that  a  usable 
calibration  array  must  be  large,  include  many  elements,  and  must  work  underwater  for  the  calibration 
to  be  meaningful  for  hydrophones. 

Schau  et  al.  [15,16]  constructed  and  tested  a  planar  array  using  a  piezoelectric  composite  material 
(Piezo  Rubber)  [17]  for  the  active  material.  The  array  measured  approximately  12.7x12.7  cm  and  was 
mounted  on  a  0.635-cm  thick  LEXAN  backing  plate  to  keep  it  planar.  The  array  was  electroded  to 
provide  20  independently  driven  active  strips.  This  array  was  operated  in  a  water  tank  at  1  kHz  and 
the  pressure  in  front  of  the  array  was  measured  with  a  small  (0.7  cm)  hydrophone.  Evanescent  fields 
were  generated  in  front  of  the  array  but  these  fields  were  heavily  contaminated  with  the  acoustic 
wavenumber.  The  pressure  field  was  also  apparently  contaminated  by  the  excitation  of  an 
antisymmetric  Lamb  wave  [18]  in  the  array  and  backing  plate  structure  with  a  phase  velocity  close  to 
that  of  the  desired  evanescent  wave. 


A  planar  array  was  constructed  and  tested  by  Trivett  et  al.  [19]  using  piezoelectric  polyvinylidene 
fluoride  (PVDF)  [20,21]  for  the  active  material.  The  PVDF  was  bonded  to  a  2.54-cm-thick  LEXAN 
backing  plate.  This  thick  LEXAN  backing  plate  was  used  to  shift  the  phase  speed  of  the 
antisymmetric  Lamb  wave  above  the  region  of  the  evanescent  waves  while  keeping  the  array 
essentially  acoustically  transparent  at  the  frequencies  of  intended  operation.  Utilizing  the  acoustic 
transparency  of  the  array,  drives  were  calculated  using  a  model  of  the  array  consisting  of  a  simple 
point-source  integrated  over  a  driven  band.  It  was  also  noted  that  the  window  width  (the  length  in  the 
simulated  flow  direction  of  the  region  in  which  the  evanescent  waves  are  to  be  created)  must  satisfy 
certain  criteria  or  the  algorithm  used  to  calculate  the  drives  would  produce  drives  that  had  a  significant 
dipole  outside  of  the  window.  This  dipole  drive  would  produce  acoustic  radiation  that  would 
contaminate  the  desired  evanescent  field.  Since  the  wavenumber  spectrum  of  the  pressure  field  being 
generated  is  a  convolution  of  the  evanescent  wavenumber  ke  with  the  wavenumber  spectrum  of  the 
window  (i.e.,  a  rectangular  aperture),  the  wavenumber  spectrum  being  generated  is  proportional  to 


sin[(*-*e)>v] 


(1) 


where  w  is  half  of  the  width  of  the  window.  If  the  evanescent  wavenumber  and  the  window  width  are 
not  chosen  such  that  Eq.  (1)  equals  zero  for  the  acoustic  wavenumber  ka  then  drives  are  generated  that 
produce  a  dipole  outside  the  window  region. 


This  planar  array  was  operated  at  frequencies  between  500  and  2000  Hz  and  evanescent  phase 
speeds  between  30  and  150  m/s.  The  results  showed  good  agreement  with  predicted  pressure  fields  as 
long  as  the  criteria  that  there  be  at  least  ten  driven  stripes  per  evanescent  wavelength  was  satisfied. 


Outline  of  the  Thesis 


This  thesis  reports  on  an  investigation  of  a  cylindrical  wavenumber  calibration  array  (CWC) 
designed  to  produce  evanescent  pressure  fields  within  the  array.  The  research  consists  of  both 
theoretical  and  experimental  investigations  of  the  array. 

The  theoretical  investigation  of  various  aspects  of  the  CWC  is  described.  First  the  pressure  field 
in  an  infinite-length  region  surrounded  by  a  cylindrical  surface  with  a  specified  axially-symmetric 
normal  velocity  is  investigated.  The  predicted  radial  dependence  of  the  pressure  is  plotted  for  several 
evanescent  wavenumbers.  The  results  indicate  the  need  for  a  CWC  with  a  low  phase  speed  within  the 
calibrator.  The  investigation  of  infinite-length  arrays  continues  with  the  calculation  of  the  phase 
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speeds  for  propagation  inside  an  idealized  array  consisting  of  an  interior  fluid  region  surrounded  by  a 
cylindrical  elastic  tube  and  an  outer  fluid  region.  A  PVDF  CWC  is  investigated  and  found  to  have 
predicted  phase  speeds  less  than  one-fifth  the  speed  of  sound  of  water.  A  finite-length  CWC  with 
discrete  band  drives  is  then  investigated.  A  method  is  described  for  using  the  electroacoustic  transfer 
matrix  of  the  system  to  compute  electrical  drives  that  will  produce  the  desired  nonacoustic- 
wavenumber  pressure  fields  within  the  CWC.  The  transfer  matrix  at  the  desired  frequency  can  be 
determined  by  measurement  or  by  using  a  computer  model  of  the  CWC.  Three  computer  models  of 
the  CWC  are  presented  that  can  be  used  to  determine  the  transfer  matrix.  The  acoustic-radiation 
program  Combined  Helmholtz  Integral  Equation  Formulation  (CHIEF)  [22]  and  the  finite-element 
program  ATI  LA  [23]  are  used  for  this  investigation.  These  models  are  combined  in  a  computer 
program  that  determines  the  CWC’s  electrical  drives  and  can  also  be  used  to  map  the  predicted 
pressure  fields  within  the  array.  These  predictions  show  a  good  evanescent  pressure  field  throughout 
most  of  the  interior  of  the  CWC. 

I  will  also  explain  the  experimental  setup  used  for  the  investigation  of  the  CWC.  First  a 
prototype  array  consisting  of  a  PVDF  tube  with  40  independently  driven  bands  is  described.  Then  an 
electronic  drive  system  is  described  that  allows  the  40  bands  to  be  driven  with  independently 
controlled  amplitudes  and  phases.  A  measurement  setup  is  then  presented  that  allows  precise 
positioning  of  a  hydrophone  within  the  water-filled  CWC  with  either  air  or  castor  oil  on  the  outside  of 
the  array.  An  electronic  measurement  system  is  then  described  that  automatically  maps  the  amplitude 
and  phase  of  the  acoustic  pressure  within  the  CWC  at  specified  locations. 

The  results  of  the  experimental  investigation  are  reported  and  measured  results  are  compared  with 
the  theoretical  predictions.  Conclusions  and  recommendations  for  future  work  are  then  presented. 

THEORETICAL  INVESTIGATION 

The  theoretical  investigation  of  the  CWC  includes  analytical  studies  of  infinite-length  cylindrical 
geometries  and  numerical  analysis  of  a  finite  length  CWC.  The  analysis  begins  by  studying  a 
nonacoustic  or  evanescent  wave  inside  an  infinite-length  cylindrical  region  bounded  by  a  surface  with 
a  specified  normal  velocity  having  only  a  z  dependence.  The  analysis  continues  with  the  study  of  the 
propagation  speed  in  the  fluid  within  an  infinite-length  elastic  tube  surrounded  by  another  fluid.  A 
finite-length  CWC  is  then  investigated  using  the  CHIEF  acoustic-radiation  program  and  the  ATILA 
finite  element  program.  These  models  of  a  finite-length  CWC  are  integrated  into  a  computer  program 
that  determines  the  electrical  band  drives  for  the  CWC  and  predicted  pressures  within  the  CWC. 

Analysis  of  Infinite-Length  Geometries 

Cylindrical  Region  with  a  Specified  Normal  Wall  Velocity 

One  area  of  concern  regarding  a  CWC  is  the  radial  uniformity  of  the  evanescent  pressure  field 
within  the  calibrator.  To  estimate  the  field  uniformity  within  an  actual  calibrator  a  mathematical 
analysis  was  performed  on  an  idealized  device  consisting  of  an  infinitely  long,  fluid-filled,  cylindrical 
region  bounded  by  a  surface  at  r  =  a.  Considering  the  harmonic,  steady-state  situation  [i.e.,  exp(-icot) 
time  dependence],  the  pressure  p  is  a  solution  to  the  Helmholtz  equation  [24] 

V2p  +  k2p  =  0  , 

where  k  is  the  magnitude  of  the  wavevector  k  and  is  equal  to  co/c,  oj  is  the  angular  frequency  equal  to 
2jt/,/is  the  frequency,  and  c  is  the  sound  speed  in  the  fluid.  In  cylindrical  coordinates  (r,9,z),  the 
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Laplacian  V2  is  given  by  [25] 


V2  =  i± 
r  dr 


d  \  Id2  d2 

r  _ 

dr. 


r2  d02  dz: 


(3) 


For  axially  symmetric  fields,  the  cylindrical  Helmholtz  equation  then  becomes 

1  d  (r  dp\  +  +  k2  Q  ' 


r  dr  I  dr 


(4) 


dz‘ 


Imposing  the  condition  that  the  pressure  be  finite  at  r  =  0,  the  solution  to  Eq.  (4)  can  be  written  as 

p(r,z)  =  AJ0(krr )  exp (ikzz)  ,  (5) 

where  A  is  a  constant  to  be  determined  by  the  boundary  condition  at  r  =  a,  J0  is  a  Bessel  function,  and 
kr  and  kz  are  the  wavevector  components  in  the  r  and  z  directions  with 


,2  j  2  , 2 

kr  +  kz  =  k 


(6) 


In  order  to  conveniently  incorporate  a  specified-velocity  boundary  condition  at  r  =  a,  the  pressure  can 
be  written  in  the  form  of  a  Fourier  integral  as  [26] 


p(r,z)  =  J ‘A(kz)  J(j \Jk  k]  rj  exp (ikgz)  dkz  . 


(7) 


Defining  the  radial  velocity  as  vr(r,z)  and  utilizing  the  linear  approximation  of  Euler’s  equation  [27] 

pin  -  -* ,  (8) 

31  dr 

where  p  is  the  density  of  the  fluid,  the  boundary  condition  at  r  -  a  becomes 


unpv(a,z)  =  p(r,z) 
dr 


(9) 


Letting  |  replace  the  integration  variable  k#  one  has 


-l  p(r,z)  =  |A(|)  {iF-  l2  J^jk2-  %2  r)  exp(i|z)  d% 


(10) 


or 
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—  p(r,z)  =  f  -A{g)  {F  f  jjjF  |2  r)  exp(i|z)  dg  . 
dr  J 


(11) 


The  boundary  condition  then  becomes 


iwpvr(a,z)  =  J-A(§)  Jk2-  g2  Jt[jF  g2  a)  exp(t^z)  dg  . 

-OO 

Multiplying  both  sides  by  exp(-ir]z)  and  integrating  over  z  from  -oo  to  oo,  gives 

00 

rtop  jvr(a,z)  exp(-/'r]z)  dz 

-00 

00  00  _  /  _  V 

=  f  f  ~A(g)  {F-  t  j^k2-  g2  a)  exp[ig-ri)z]  dgdz 


(12) 


(13) 


Utilizing  the  identity  [28] 


Jexp(taz)  dz  -  2jiS(a) 


where  6(a)  is  the  Dirac  delta  function  [29],  and  the  fact  that 

00 

J7(S)  8(1  -ri)^-/Ol) 


(14) 


(15) 


Eq.  (13)  becomes 


icop  J vr(a,z)  exp(-rr]z)  dz  =  -2it  A(r])  \[k 2-  r|2  jJjF-  r\2  a 


(16) 


Solving  forA(ri)  gives 


A(ri)  = 


-icop  jvr(a,z)  exp(-ir)z)  dz 
-00 

2k  )Jk2-  rj2  Jl[\[F-  r\2  a 


The  normal  velocity  on  the  surface  is  specified  as 


(17) 
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vr(a,z)  =  v0  exp (ikez)  ,  (18) 

where  ke  is  the  desired  evanescent  wavenumber.  The  evanescent  wavenumber  is  related  to  the  acoustic 
wavenumber  k  by 

ke  =  fik  (a  >  1)  ,  (19) 

where  n  is  the  wavenumber  multiplier.  With  this  velocity,  one  has 


00  00 

jvr(a,z)  exp {-ikzz)  dz  =  v0  Jexp  [i(ke-k^z]  dz  . 
-00  -00 


(20) 


Again  utilizing  the  identity  given  in  Eq.  (14),  this  becomes 


00 


Jvr(a,z)  exp( 


-ikzz)  dz  =  2jtv0  b(ke-kz)  . 


(21) 


Inserting  this  result  into  Eq.  (17)  yields 


A(kz)  = 


-i(opvc  6  (ke-kz) 


The  pressure  is  then,  from  Eqs.  (7)  and  (22), 


(22) 


P(r,z)  =  f 


-i(ppvo  jjjk2-  kz  r 

kz  Ji  {ft*-  kz  fl) 


exp (ikzz)  6 (kg-kz)  dkz  . 


(23) 


Once  again  utilizing  the  characteristic  of  the  Dirac  delta  function  given  in  Eq.  (15),  and  Eq.  (19)  one 
obtains 


P(r,z)  = 


~ipcvo  Jlkh~M2  r) 

h~P2  J^h-P2  a 


exp  (ik  z)  . 


(24) 


For  evanescent  waves  n  is  greater  than  one  making  the  arguments  of  the  Bessel  functions 
imaginary.  This  necessitates  the  use  of  the  modified  Bessel  functions  I0  and  7j  [30].  The  pressure 
inside  the  cylindrical  region  is  then 
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P(r,z)  = 


~Pcvo  r, 


exP(^eZ)  • 


(25) 


Fig.  2  —  Radial  pressure  dependence  for  various  wavenumber  multipliers. 


The  pressure  shows  the  desired  z  dependence  and  a  radial  dependence  that  has  a  maximum  at  r  - 
a  decaying  to  a  minimum  at  r  =  0.  The  results  of  this  analysis  are  shown  in  Fig.  2.  The  plot  shows 
the  radial  pressure  dependence  for  various  wavenumber  multipliers  for  a  water-filled  region  at  a 
frequency  of  1  kHz  with  a  wall  radius  a  of  1.5  cm.  The  results  show  that  the  pressure  drops  only 
slightly  for  small  wavenumber  multipliers  but  decays  rapidly  for  large  multipliers. 
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For  a  calibration  device,  it  is  desirable  for  the  pressure  field  to  be  as  uniform  as  possible.  One 
reason  for  this  is  the  desire  to  have  the  highest  pressure  possible  in  the  central  portion  of  the 
calibrator,  where  the  device  being  calibrated  is  placed,  relative  to  the  inner  surface  of  the  calibrator.  A 
rapidly  decaying  field  also  restricts  the  size  of  the  device  being  calibrated  and  makes  the  positioning  of 
the  device  very  critical. 

One  way  to  improve  the  uniformity  of  the  pressure  field  is  to  reduce  the  natural  propagation 
speed  within  the  calibrator.  This  allows  the  creation  of  pressure  fields  with  specified  evanescent 
wavenumbers  using  smaller  wavenumber  multipliers.  The  smaller  multipliers  then  result  in  more 
uniform  pressure  fields. 

At  a  specified  frequency  the  propagation  speed  within  the  calibrator  depends  on  two  things:  (1) 
the  fluid  within  the  calibrator  and  (2)  the  compliance  of  the  calibrator  wall.  Because  the  devices  that 
will  be  calibrated  are  normally  operated  in  water,  the  calibrator  will  be  filled  with  water  when  in  use. 
Even  if  another  fluid  was  used,  most  convenient  fluids  have  sound  speeds  that  do  not  vary  greatly 
from  that  of  water.  The  only  other  way  to  significantly  lower  the  propagation  speed  is  to  design  a 
calibrator  with  a  fairly  flexible  wall. 

Propagation  Speed 

To  determine  the  predicted  propagation  speed  within  the  PVDF  CWC  prototype,  the  propagation 
speed  in  the  fluid  within  an  infinite-length  elastic  tube  surrounded  by  another  fluid  was  studied.  The 
details  of  this  analysis  are  given  in  Appendix  A.  To  verify  the  predicted  propagation  speeds 
experimentally,  an  effort  was  undertaken  to  measure  the  propagation  speeds  within  the  CWC.  The 
relatively  short  length  of  the  CWC  and  low  frequencies  involved  meant  that  it  was  not  feasible  to 
determine  the  propagation  speeds  with  time-of-flight  measurements.  The  first  four  resonant 
frequencies  for  the  water  column  within  the  CWC  were  measured  with  the  CWC  air  backed  and 
castor-oil  backed.  Using  these  frequencies  and  the  effective  length  [31]  of  the  CWC  the  propagation 
speed  within  the  CWC  could  then  be  calculated  for  each  case. 

The  results  of  this  analysis  for  the  PVDF  CWC  prototype  with  either  air  or  castor  oil  for  the 
outer  fluid  are  shown  in  Fig.  3.  The  figure  shows  the  propagation  speed  in  the  water  within  the  CWC 
as  a  function  of  frequency  throughout  the  operational  frequency  range  of  the  drive  electronics.  This 
frequency  range  is  below  the  cutoff  frequency  for  all  modes  in  the  water  except  the  zeroeth  mode 
The  cutoff  frequency  for  the  next  mode  was  calculated  [32]  to  be  approximately  35  kHz.  The 
measured  propagation  speeds  at  the  first  four  resonances  show  good  agreement  with  the  theoretical 
speeds.  The  reduction  in  the  propagation  speed  with  the  oil  on  the  outside  of  the  CWC  is  caused  by 
the  mass  loading  of  the  oil.  At  very  low  frequencies  this  effect  is  of  little  importance  because  inertial 
forces  are  small  compared  to  elastic  forces  [33]. 

It  can  be  seen  that  the  propagation  speed  within  the  PVDF  CWC  prototype  is  about  one-fifth  the 
sound  speed  in  open  water  (1485  m/s).  This  would  indicate  that  pressure  fields  for  a  given  evanescent 
wavenumber  could  be  created  within  the  PVDF  CWC  with  much  less  radial  pressure  dependance  than 
within  a  more  rigid  calibrator  where  the  propagation  speed  would  be  closer  to  the  speed  of  sound  in 
water. 

Determination  of  Band  Drives 

Operation  of  the  CWC  requires  the  determination  of  appropriate  sets  of  band  drives  that  will 
produce  acceptable  high-wavenumber  fields  within  the  CWC.  The  technique  that  was  used  was  to  first 
determine  the  electroacoustic  transfer  matrix  for  the  system  and  then  use  the  inverse  transfer  matrix  to 
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Fig.  3  —  Theoretical  and  measured  propagation  speeds  in  the  water  column  within  the 
PVDF  CWC  with  air  or  castor  oil  outside  the  CWC. 


calculate  a  set  of  drives  that  will  produce  a  desired  field  [34]. 

Assume  that  there  are  n  driven  bands  each  being  driven  by  a  sinusoidal  signal  diy  i  =  1,  2,  .  .  n. 
Also  assume  that  the  pressure  can  be  determined  at  m  field  points  within  the  CWC  where  m  £  n.  The 
pressure  at  the  field  point  may  be  determined  by  actual  measurement  or  it  may  be  determined  using  a 
computer  model  of  the  device.  The  pressure  at  the ;th  field  point  (where;  =  1,  2,  .  .  m)  has  a 
complex  value  pj.  If  the  system  is  linear  and  superposition  holds,  one  can  define  the  electroacoustic 
transfer  matrix  by 
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p  =  Md  ,  (26) 

where  the  transfer  matrix  M  is  an  mxn  complex  matrix.  The  first  task  is  to  determine  the  elements  of 
M.  This  can  be  accomplished  by  driving  each  band  alone  and  measuring  the  resultant  pressures  at 
each  field  point  First  one  applies  drive  to  the  first  band,  while  setting  all  the  other  band  drives  to 
zero,  and  determine  the  resulting  pressure  at  each  of  the  field  points.  One  can  then  obtain  the  first 
column  of  M  from 


Mn  =  pj/rfj 
^21  =  Pl^l 


.  dj  -  0  for  j*!. 


Mml  =  Pm!d  1 


(27) 


This  process  is  repeated  by  driving  the  remaining  bands,  one  at  a  time,  to  fill  each  column  of  M. 

One  must  then  determine  the  set  of  drive  levels  that  best  approximates  the  desired  sound  field  at 
the  field  points.  This  can  be  accomplished  by  noting  that 

P  ■  MD  ,  (28) 


where  P  is  a  vector  of  dimension  m  consisting  of  the  desired  pressures  at  the  field  points,  and  D  is  a 
vector  of  dimension  n  consisting  of  drives  that  will  produce  the  desired  pressures  P.  If  m  =  n  the 
drive  levels  D •  can  be  computed  directly  from 


D  =  M~l  P  . 


(29) 


If  m  >  n  one  can  determine  the  drives  that  produce  a  pressure  field  that  best  fits  the  desired  pressure 
field  P  in  a  least-squares  sense  [35]  (i.e.,  where  the  squares  of  the  errors  are  minimized)  from  [36] 

D  =  (MT*M)'1MT*P  ,  (30) 


where  the  superscript  T  indicates  transposition  and  the  asterisk  indicates  complex  conjugation.  This 
can  be  expressed  in  the  simpler  form 


D  =  (Af)"V  ,  (31) 

where  Mf  =  AfT*M  is  the  reduced  transfer  matrix  and  P’  =  A/T*P  is  the  reduced  vector  of  pressures.  It 
should  be  noted  that  if  m  >  n,  M  is  not  square  but  M'  is  always  square  and  hence  invertible 

Although  the  drives  computed  by  this  method  will  mathematically  produce  the  pressure  field 
closest  to  the  desired  field  (in  a  least-squares  sense),  they  may  not  be  the  most  suitable  set  of  drives  in 
practice  [37].  The  computed  drive  levels  can  vary  dramatically  from  one  band  to  the  next.  In  this 
situation  small  changes  in  the  drive  levels  can  cause  large  changes  in  the  uniformity  of  the  resultant 
pressure  field.  This  places  severe  demands  on  the  dynamic  range  and  precision  of  the  electronic 
system  used  to  produce  the  drives.  An  example  of  this  condition  is  presented  later. 
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One  solution  to  this  problem  was  first  developed  for  antenna  pattern  synthesis  [38].  Equation 
(31)  can  be  expressed  as 


n 


E 

;  =  i 


(M* 


=  P. 


i  =  1,2 . 


(32) 


If  a  small  positive  real  constant  a  is  added  to  the  diagonal  elements  of  AT  one  obtains  the  approximate 
equation 


E(a8 $  +  M  .JlDj  ~  p,  ,  i  -  1,2,...  ,n  , 


(33) 


where  6-  is  the  Kronecker  delta.  Although  the  drive  levels  tend  to  decrease  as  a  increases,  the  larger 
drive  levels  decrease  more  rapidly  than  the  smaller  drive  levels,  thereby  creating  more  uniform  drives. 
Even  though  these  drive  levels  are  no  longer  the  ideal  least-squares  solution,  small  values  of  a  often 
produce  drives  that  are  more  suitable  in  practice  for  creating  the  desired  pressure  field.  A  wide  range 
of  values  of  a  will  produce  a  significant  improvement  in  the  uniformity  of  the  drive  levels  with  only 
slight  degradation  of  the  uniformity  of  the  pressure  field.  In  practice  an  a  that  is  one-thousandth  the 
average  amplitude  of  the  M*  elements  usually  gives  good  results. 


Utilizing  this  technique,  the  solution  for  D  then  becomes 

D  -  (a/  +  M')~lP'  , 


(34) 


where  I  is  a  unit  matrix.  It  should  be  noted  that  once  Af  is  determined  for  a  particular  frequency,  a 
set  of  drives  for  any  desired  field  can  be  acquired  by  resolving  Eq.  (34)  with  the  appropriate  values  for 
P  and  therefore  P\ 


Computer  Program  to  Model  the  CWC  and  Determine  Band  Drives 

As  noted  above,  the  drives  for  the  bands  can  be  calculated  from  a  transfer  matrix  for  the  system 
that  is  determined  by  either  measuring  the  pressures  within  the  CWC  as  each  band  is  driven  or  by 
using  a  computer  model  of  the  CWC  and  calculating  the  predicted  pressures  using  the  model.  An 
overall  view  of  the  prototype  PVDF  CWC  is  shown  in  Fig.  4.  The  piezoelectric  element  of  the  CWC 
consisted  of  a  PVDF  tube  with  an  inner  radius  of  1.55  cm,  a  wall  thickness  of  0.75  mm  and  an  overall 
length  of  25.2  cm.  The  PVDF  was  electroded  on  the  inside  and  outside  of  the  tube  in  such  a  way  that 
it  could  be  electrically  driven  as  40  bands,  each  having  a  width  of  6  mm,  with  an  additional  grounded 
band  at  each  end  of  the  tube.  A  FORTRAN  program  called  CWCMDL  was  written  to  model  the 
CWC  and  calculate  the  drives  using  three  different  methods  of  determining  the  transfer  matrix.  The 
first  method  utilizes  a  CHIEF  model  of  the  CWC  specifying  an  interior- wall  impedance.  The  second 
method  models  the  bands  of  the  CWC  as  thin-walled  piezoelectric  rings  and  includes  the  effects  of 
fluid  loading  by  using  a  radiation-impedance  matrix  computed  with  CHIEF.  The  third  method  utilizes 
the  ATILA  finite-element  code  coupled  with  CHIEF  to  include  the  effects  of  the  fluid. 

CWC  Models  Using  Only  CHIEF 

CHIEF  [39,40]  is  a  computer  program  written  in  FORTRAN  that  computes  the  acoustic  radiation 
from  arbitrarily-shaped  radiating  bodies.  The  program  consists  of  a  collection  of  subroutines  that  are 
called  by  a  user-written  control  program.  CHIEF  utilizes  the  fact  that  the  pressure  and  velocity  can  be 
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Fig.  4  —  Prototype  PVDF  CWC. 
related  through  the  classical  Helmholtz  integral  equation  [41,42] 


exp[-r&r(jt,cr)] 

r(x,o) 


r(x,o) 


dS(a) 


0  ,  x  interior  to  S 
2ji/?(jc),  x  on  S 
4np(x),  x  exterior  to  S 


(35) 


where  S  is  the  radiating  surface,  d/d na  is  the  outward-normal  derivative,  r(x,a)  is  the  distance  between 
the  surface  point  a  and  the  field  point  x,  w  is  the  radial  frequency,  and  p  is  the  fluid  density. 


The  operation  of  CHIEF  can  be  divided  into  two  basic  parts: 

•  With  the  normal  surface  velocity  specified,  calculate  the  surface  pressure  using  the 
Helmholtz  integral  equation  at  the  surface. 

•  Calculate  the  pressure  at  any  point  in  the  surrounding  medium  using  the  Helmholtz 
integral  for  the  surrounding  medium. 


To  solve  the  Helmholtz  integral  equation,  the  surface  is  divided  into  N  subdivisions  over  which 
the  pressure  and  normal  velocity  are  considered  to  be  constant.  The  Helmholtz  integral  equation  at  the 
surface  is  then  reduced  to  a  system  of  N  complex  linear  equations  in  N  unknowns.  CHIEF  can  handle 
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arbitrarily  shaped  bodies,  but  a  subroutine  is  provided  which  simplifies  the  process  for  many  common 
geometries  [43].  Symmetry  can  be  utilized  to  reduce  the  number  of  required  subdivisions.  CHIEF 
requires  the  numerical  solution  of  2N 2  surface  integrals.  The  system  is  then  solved  for  the  surface 
pressures  at  each  of  the  N  subregions. 

There  are  certain  frequencies  called  Characteristic  frequencies”  for  which  the  solution  of  the  Eq. 
(35)  is  not  unique  [44].  These  frequencies  are  those  for  which  there  exists  a  nontrivial  function 
satisfying  the  Helmholtz  equation  throughout  the  volume  enclosed  by  the  radiating  surface  and 
vanishing  at  the  radiating  surface. 

Schenck  proposed  a  solution  to  this  problem  which  is  utilized  by  CHIEF  [39].  Since  the  solution 
to  the  integral  equation  is  not  unique  at  certain  frequencies,  unique  solutions  can  be  obtained  by 
placing  additional  constraints  on  the  pressure.  CHIEF  requires  that  Eq.  (35)  be  satisfied  for  points 
interior  to  S  that  are  specified  by  the  user.  CHIEF  then  uses  a  Householder-reduction  least-squares 
technique  [45]  to  compute  the  surface  pressures. 

There  remains  the  problem  of  determining  how  many  interior  points  are  required  and  where  to 
position  them.  CHIEF  includes  a  point-source  check  [46]  to  help  determine  the  validity  of  its  results. 
This  check  allows  the  user  to  define  one  or  more  point  source  locations.  Normal  surface  velocities  for 
each  of  the  N  surface  subdivisions  are  assigned  that  coincide  with  the  normal  velocity  that  would  be 
generated  at  the  location  of  the  subdivisions  by  the  point  sources  if  the  vibrating  body  were  absent. 

The  surface  pressures  are  also  computed  by  CHIEF  from  Eq.  (35).  These  surface  pressures  are  then 
compared  with  the  pressures  that  would  be  generated  at  the  location  of  the  surface  subdivisions  by  the 
point  sources  in  the  absence  of  the  vibrating  body.  The  point-source  check  can  help  determine 
whether  sufficient  interior  points  are  specified  and  if  the  surface  subdivision  scheme  is  adequate  and 
correct.  However,  the  point-source  check  can  lead  to  erroneous  conclusions  if  the  velocity  distribution 
due  to  the  point  sources  varies  too  much  from  the  velocity  distribution  that  normally  exists  on  the 
surface.  These  potential  problems  indicate  that  the  results  from  CHIEF  require  careful  analysis. 

The  prototype  CWC’s  largest  internal  dimension  (i.e.,  its  length)  is  25.2  cm.  The  lowest 
characteristic  frequency  for  the  CWC’s  internal  modes  is  approximately  3  kHz,  which  is  above  its 
operating  frequency.  This  eliminates  the  problem  of  frequencies  associated  with  internal  modes  for  the 
radiating  body  and  no  interior  points  are  needed. 

The  task  of  writing  a  control  program  that  calls  the  CHIEF  subroutines  has  been  greatly 
simplified  with  the  CHIEF  Preprocessor  called  CID  (CHIEF  Interactive  Driver)  [47].  This  program 
allows  the  user  to  interactively  answer  questions  about  the  desired  CHIEF  analysis  and  then  generates 
a  control  program  for  CHIEF. 

A  recent  enhancement  to  CHIEF  is  the  capability  of  calculating  the  radiation  from  vibrating 
bodies  for  which  part  of  the  surface  has  an  arbitrary  acoustic  impedance  rather  than  a  specified  normal 
velocity  [48,49],  This  feature  is  used  for  one  model  of  the  CWC.  All  the  CHIEF  models  of  the  CWC 
use  18-fold  rotational  symmetry  about  the  axis  of  the  CWC. 

CWC  Model  with  All  Surface  Velocities  Specified 

One  CHIEF  model  of  the  CWC  that  CWCMDL  utilizes  involves  the  use  of  CHIEF  in  the  traditional 
manner  (i.e.,  specifying  the  normal  velocity  of  every  surface  element  of  the  radiating  body).  For  this 
model  the  CWC  is  divided  into  86  surfaces  as  shown  in  Table  1.  This  table  shows  a  cross  section  of 
the  wall  of  the  CWC  (not  to  scale).  The  inside  of  the  CWC  is  on  the  left  and  the  outside  is  on  the 
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Table  1  —  Surface  numbering  and  axial  positions  of  the  borders  of  the  CHIEF 
surfaces  for  the  model  with  all  the  normal  velocities  specified. 


POSITION 

CHIEF 

44 

CHIEF 

BAND 

i 

0.126 

43 

86 

GND 

0.120 

42 

85 

40 

0.114 

-0.096 

6 

49 

4 

0.102 

5 

48 

3 

0.108 

4 

47 

2 

0.114 

3 

46 

1 

0.120 

2 

45 

GND 

0.126 

1 

1 

right  The  inner  radius  is  1.55  cm  and  the  outer  radius  is  1.625  cm.  Positions  shown  are  the  locations 
of  the  divisions  between  surface  elements  in  meters  along  the  axis  of  the  CWC  (the  z  axis)  with  the 
origin  at  the  center  of  the  CWC.  Surface  element  widths  of  6  mm  were  chosen  to  coincide  with  the  6- 
mm  electrode-band  spacing.  The  electroding  is  approximated  by  double  lines.  The  band  numbering 
and  CHIEF  surface-area  numbering  are  also  shown.  Since  all  the  normal  surface  velocities  are 
specified  with  this  model,  the  description  of  the  surfaces  remains  the  same  no  matter  what  surface 
velocities  are  specified.  This  means  that  the  surfaces  and  the  matrices  for  CHIEF  only  have  to  be  set 
up  and  calculated  once  for  each  program  run.  The  field  pressures  for  different  sets  of  velocities  can 
then  be  calculated  without  rerunning  all  of  CHIEF.  This  model  is  used  to  determine  the  field 
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pressures  whenever  all  the  normal  surface  velocities  are  known  (i.e.,  after  the  normal  surface  velocities 
have  been  determined  using  another  model  of  the  CWC). 

CWC  Model  with  Interior-Wall  Impedance  Specified  (CHIEF-I) 

Another  CHIEF  model  is  used  to  calculate  the  field  pressures  when  a  single  band  is  driven  where 
a  surface  impedance  is  specified  for  the  remainder  of  the  interior  wall.  This  model  seemed  plausible 
because,  while  making  the  pressure  measurements  within  the  CWC  to  determine  the  resonant 
frequencies  and  therefore  the  propagation  speeds  with  air  or  water  outside  the  CWC,  it  was  noted  that 
there  seemed  to  be  very  little  coupling  between  the  pressures  fields  inside  and  outside  of  the  CWC.  It 
was  felt  that  a  simple  model  that  matched  the  conditions  within  the  CWC  reasonably  well  would  give 
useful  results  even  if  the  model  did  not  match  the  conditions  well  on  the  outside  of  the  CWC.  Since 
this  model  ignores  the  effect  of  acoustic  loading  on  a  driven  band  it  is  only  used  to  model  the  CWC 
when  a  single  band  is  driven  and  therefore  cannot  be  used  to  model  the  predicted  performance  of  the 
CWC  when  all  the  bands  are  driven  with  a  set  of  computed  drives. 

centerline  of  CWC 


v=v. 


z=z, 


z=z. 


v  =  0  I  I  I  I  Kfl  I  I  I  I . .  . . Ill  I  I  I  iTTTTTI  v  =  0 


v=0 

Fig.  5  —  Diagram  of  the  CHIEF  model  of  the  CWC  with  a  wall  impedance  specified 
on  the  interior  of  the  CWC  except  for  the  interior  of  the  driven  band. 


A  diagram  of  this  model  is  shown  in  Fig.  5  for  band  number  4  being  driven.  (Note:  The  first 
band  is  a  ground  band,  therefore,  driven  band  4  is  the  fifth  band  in  from  the  left.)  In  this  model: 

1.  The  normal  velocities  on  the  outside  and  the  end  surfaces  of  the  CWC  are  set  to  zero. 

2.  The  normal  velocity  of  the  inner  surface  of  the  driven  band  v0  is  specified. 

3.  The  surface  impedance  of  the  remainder  of  the  interior  wall  of  the  CWC  zi  is  specified. 
The  interior-wall  impedance  z(  is  calculated  from  [50] 


Yt.. 


zi  = 


ma 


(36) 


where  Y  is  the  Young’s  modulus,  tw  is  the  wall  thickness,  and  a  is  the  inner  radius  of  the  PVDF  tube. 

CHIEF  requires  that,  if  arbitrary  impedances  are  specified  instead  of  normal  surface  velocities  for 
any  of  the  surfaces,  they  must  be  the  last  surface  regions  described.  Since  impedances  are  specified 
for  the  interior  wall  regions  below  and  above  the  driven  band,  they  must  be  the  last  surfaces  loaded. 
The  locations  of  these  surfaces  and  the  inner  surface  of  the  driven  band  change  whenever  a  different 
band  is  driven,  so  these  surfaces  are  reloaded  and  the  matrices  recalculated  as  each  different  band  is 
driven.  This  means  that  CHIEF  has  to  perform  more  computations  to  determine  the  transfer  matrix  by 
driving  each  band  with  this  model  than  in  the  model  where  all  the  surface  velocities  are  specified. 
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The  CHIEF  model  of  the  CWC  with  specified  interior-wall  impedances  will  be  referred  to  as  the 
CHIEF-I  model. 

RING -CHIEF  Model 

One  limitation  of  the  interior-wall-impedance  model  of  the  CWC  is  the  fact  that  the  effects  of  the 
acoustic  loading  on  the  velocity  of  the  driven  band  are  ignored.  Another  model  of  the  CWC  called 
RING-CHIEF  was  created  that  models  the  CWC  bands  as  thin-walled  piezoelectric  rings  and  includes 
the  effects  of  the  acoustic  loading. 


Fig.  6  —  Thin-walled  piezoelectric  ring  used 
for  the  RING-CHIEF  model. 


The  radiation  impedance  of  a  fluid-loaded  vibrating  surface  is  the  ratio  of  the  force  exerted  by 
the  surface  on  the  fluid  to  the  normal  velocity  of  the  surface  [51].  The  total  force  on  the  vibrating 
surface  also  includes  the  effects  of  the  internal  mechanical  impedance  and  the  mutual  radiation 
impedances  from  other  surfaces.  One  can  model  the  bands  of  the  CWC  as  thin-walled  piezoelectric 
rings  with  each  ring  having  an  inner  radius  a,  a  wall  thickness  t ^  and  a  height  h  (see  Fig.  6).  A 
coordinate  system  is  used  in  which  1  is  in  the  circumferential  direction  and  3  is  in  the  radial  direction. 
Each  ring  is  electroded  on  the  inner  and  outer  surfaces  and  poled  radially.  If  one  assumes  that  there 
are  N  bands,  and  therefore  2N  total  inner  and  outer  surfaces,  the  radial  force  exerted  on  band  i  is  [52] 


+  . 


inner  ^  inner  Pouter  A  OL 


outer  ’ 


(37) 


where 

zL  = 

mechanical  impedance  of  the  ith  band 

vi  = 

radial  velocity  of  the  ith  band 

i 

Pinner  “ 

pressure  on  the  inner  surface  of  the  ith  band 

i 

P  outer  “ 

pressure  on  the  outer  surface  of  the  ith  band 

A* 

inner 

area  of  the  inner  surface  of  the  ith  band 

A  1 

outer 

area  of  the  outer  surface  of  the  ith  band 
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The  pressures  on  the  inner  and  outer  surfaces  of  the  band  can  be  expressed  as 


and 


27V 

p.‘  =  V  G'  -'v 

r  inner  /  j  j 

;  =  i 


(38) 


2  TV 

PLter  =  E  G,*N'j 

7  =  1 


(39) 


where  G  *•  ’  is  a  Green’s  function  [53]  and  v;  is  the  normal  velocity  of  surface  j.  If  the  inner  and  outer 
surface  areas  are  written  as 


and 


A! 


-  A  . 
inner  i 


A  outer  ~  Ai*N  > 


(40) 


(41) 


Eq.  (37)  can  be  expressed  as 

2N 

Fi  =  Z>,  ♦  E  (G‘JAi  ~  G‘*NJAi+N)vj  •  (42) 

7  =  1 

The  summation  can  be  divided  into  separate  parts  for  the  inner  and  outer  normal  surface  velocities 
giving 


or 


N 


Fi = zv  +  y  (g‘ja,  -  Gi*NjAi+N)vj 

j= i 


2  TV 


*  E  -  c"a,J'W'7 


(43) 


^  +  E  -  Gi*N’iAi+N)vj 

7=1 


TV 


♦  Y  (PiJ*NAi  -  ■ 

7=1 


(44) 


Since  the  bands  are  assumed  to  be  thin-walled  rings,  the  outer  wall  normal  surface  velocities  are  equal 
to  the  negative  of  the  inner  wall  normal  surface  velocities 


Vj*N  -  ~Vj 


(45) 


giving 
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Fi  =  Z>,  +  E  (G‘’JAi  ~  G‘*N'iAi*N  ~  gWA;  *  Gi+N^NA^N)Vj  (46) 

;'  =  i 

or 

N 

Fj  =  Z^v;  +  Y,  (Z‘J  ~  z‘*N’j  ~  Z‘J*N  +  Z‘*N’**N)Vj  ,  (47) 

>■» 

where  the  radiation  impedance  Z*;  =  Gl}Ar 

The  radiation  impedance  Zt}  is  defined  as  the  normal  force  exerted  on  surface  i  when  surface  j 
moves  with  unity  normal  velocity  and  all  other  surfaces  have  normal  velocities  of  zero.  This  coincides 
with  the  radiation-impedance  matrix  [Z]  that  can  be  computed  with  CHIEF  [54].  Equation  (47)  can  be 
written  as 


where 


N 


Fi  -  EZ<>, 
>  =  1 


(48) 


Z 


*; 

to* 


zi,;  _  _  Z'j+Ar  +  z‘*N'i+N  ,  j  *  i 

Z*’*  -  Z,+Ar’'  -  Z'-'^  +  Z‘*N'GN  +  Z '  ,  i  =  i 

tn  '  J 


Equation  (48)  can  then  be  expressed  in  matrix  notation  as 


(49) 


*  =  IZto/3  T  , 


(50) 


where 

F  =  vector  of  the  normal  surface  forces 

\Zjot\  =  total  impedance  matrix 
V  =  vector  of  the  normal  surface  velocities 

This  system  can  be  solved  for  the  resultant  velocities  for  a  set  of  specified  forces. 

The  strain  in  a  thin-walled  ring  can  be  expressed  as  [55] 

S1  =  sflTl  +  >  (51) 


where 

*f  = 

elastic  compliance  coefficient  with  constant  electric  field 

*1  = 

circumferential  stress 

^31  = 

piezoelectric  strain  coefficient 

e3  = 

electric  field 

Considering  the  blocked  case  (i.e.,  Sj  =  0)  and  noting  that  sfi  is  the  reciprocal  of  Young’s  modulus 
[56]  Y,  Eq.  (51)  becomes 
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7’i  =  -rf31e3r.  (52) 

This  stress  exerts  a  radial  force  on  the  ring  of  [55] 

F-t  =  - 2nhtwTi  .  (53) 

Substituting  Eq.  (52)  into  this  gives 

Fi  =  2xhtwdn  e3y  .  (54) 


If  a  voltage  Ei  is  applied  across  the  electrodes,  the  fact  that  e3  =  EJtw  gives  a  force  of 

Ft  =  2jt hd^YE;  . 


(55) 


The  internal  mechanical  impedance  of  the  thin  ring  can  be  expressed  as  [57] 


Z 


m 


t 

i  coM  - 

\ 


(56) 


where  M  =  2jtp ahtw  is  the  mass  of  the  ring  and  Cm  -  a/(2iihtwY)  is  the  mechanical  compliance  of  the 
ring. 


The  normal  surface  velocities  of  all  the  surfaces  when  a  single  band  is  driven  can  be  computed 
from  Eq.  (50)  by  setting  Fi  =  0  for  all  bands  except  the  driven  band.  The  force  on  the  driven  band  is 
given  by  Eq.  (55).  The  system  in  Eq.  (50)  is  then  solved  for  the  velocities.  A  FORTRAN  program 
was  written  to  compute  the  necessary  quantities  and  solve  the  resultant  system  for  the  surface 
velocities.  The  model  is  implemented  by  driving  the  first  band  alone  with  =  1  and  solving  for  the 
resultant  normal  surface  velocities,  which  are  then  stored  in  a  file.  This  process  is  repeated  for  each  of 
the  driven  bands  through  band  20  with  each  set  of  resultant  surface  velocities  appended  to  the  velocity 
file.  Since  the  CWC  is  symmetric,  the  resultant  surface  velocities  for  drives  to  bands  21  through  40 
can  be  obtained  from  the  results  for  bands  1  through  20.  The  file  containing  the  resultant  surface 
velocities  for  unity  drive  to  each  of  the  first  twenty  bands  is  then  used  to  calculate  all  the  surface 
velocities  for  any  set  of  band  drives.  This  is  accomplished  by  multiplying  the  complex  contribution  to 
each  surface  velocity  for  unity  drive  to  each  band  by  the  complex  drive  to  that  band.  The  complex 
velocities  are  then  summed  for  all  driven  bands  to  give  the  resultant  velocities. 

ATILA-CHIEF  Model 

In  order  to  obtain  a  more  complete  model  of  the  PVDF  CWC,  the  ATILA  finite-element  code 
was  used.  This  model  includes  structural  interaction  within  the  PVDF  tube  and  the  effects  of  the  fluid. 
ATILA  was  specifically  developed  to  aid  in  the  design  of  sonar  transducers,  but  can  be  used  for  other 
types  of  transducers  or  for  passive  structures.  If  one  first  considers  the  harmonic  analysis  of  the  CWC 
in  vacuum  (i.e.,  without  any  effects  of  fluid  loading),  the  matrix  equation  describing  the  problem  can 
be  written  as  [58] 
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[Kuu}-i»2[M]  [K^] 


IV 


u 

o 


F 

-q 


where 


\KUU] 

- 

stiffness  matrix 

[*«*] 

= 

piezoelectric  matrix 

[Kh] 

= 

dielectric  matrix 

[M] 

= 

mass  matrix 

U 

= 

vector  of  the  nodal  values  of  the  displacement  field 

- 

vector  of  the  nodal  values  of  the  electrical  potential 

F 

= 

vector  of  the  nodal  values  of  the  applied  forces 

q 

= 

vector  of  the  nodal  values  of  the  electrical  charges 

(57) 


and  superscript  T  indicates  matrix  transposition.  For  a  structure  that  is  in  contact  with  a  fluid  medium 
the  force  vector  F  can  be  written  as  [59] 


F  =  Fr  +  Fd  ,  (58) 

where  FR  is  associated  with  the  acoustic  loading  of  the  fluid  medium  and  FD  is  associated  with  the 
other  known  forces  that  drive  the  structure. 

The  effects  of  the  fluid  loading  can  then  be  included  by  using  the  radiation-impedance  matrix 
that  relates  forces  and  velocities  on  the  wetted  surfaces  by  [60] 

[Z]V  =  F  ,  (59) 


where 


V 

F 


radiation-impedance  matrix 

vector  of  the  normal  surface  velocities 

vector  of  the  normal  surface  forces 


n 

The  nodal  forces  F  can  be  expressed  in  terms  of  the  normal  surface  force  for  surface  S-  by  [61] 


F  =  - 


5| 


(60) 


or 


where 


(61) 
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*'  jWJT-n,dS,, 

'  Si 

area  of  surface  5, 
interpolation  function 
normal  unit  vector  for  surface  Si 


(62) 


The  can  be  reordered  in  a  larger  vector  Jt(  containing  all  the  degrees  of  freedom.  A  coupling 
matrix  [AT]  can  be  defined  as 


[X]  m  [jivn2,...,x  J  , 


(63) 


where  n  is  the  total  number  of  wetted  surfaces.  The  nodal  forces  are  related  to  the  nodal 
displacements  by 

Fr  =  -  i  u)  [X]  [Z]  [X]T  V  .  (64) 

The  finite-element  model  then  becomes  [62] 

[Kuu)-u2\M\  +  m[X][Z\[X?  [K^) 

After  solving  this  set  of  equations,  the  normal  velocities  of  the  wetted  surfaces  are  computed  from 

Y  =  -iw[X\TU  .  (66) 

CHIEF  can  then  be  used  to  compute  field  pressures  using  the  normal  surface  velocities. 

ATILA  incorporates  an  option  that  will  read  a  radiation-impedance  matrix  file  created  by  CHIEF 
and  compute  the  resultant  normal  velocities  of  the  wetted  surfaces.  Care  must  be  taken  to  describe  the 
radiating  body  in  CHIEF  and  ATILA  in  such  a  way  that  the  surfaces  in  CHIEF  correspond  to  the 
appropriate  elements  in  ATILA. 

The  procedure  for  calculating  field  pressures  using  the  ATILA-CHIEF  model  is  therefore: 

1.  Run  CHIEF  to  calculate  and  store  the  appropriate  radiation-impedance  matrix. 

2.  Run  ATILA  with  appropriate  boundary  conditions.  ATILA  uses  the  radiation- 
impedance  matrix  and  calculates  the  resultant  normal  surface  velocities. 

3.  Run  CHIEF  (specifying  the  normal  surface  velocities  determined  with  ATILA)  to 
calculate  the  resultant  field  pressures. 

An  ATILA-CHIEF  model  of  the  CWC  was  created  utilizing  the  axial  symmetry.  The  flow 
diagram  for  the  model  is  shown  in  Fig.  7.  The  region  inside  the  dashed  lines  is  labeled  GEOMETRY 
DEPENDENT  to  signify  that  it  only  depends  on  the  physical  parameters  of  the  system  (e.g., 
dimensions,  material  properties,  frequency,  etc.).  Therefore,  the  ATILA-CHIEF  computation  only  has 
to  be  performed  once  for  a  particular  set  of  physical  parameters.  The  results  can  then  be  used  to 
determine  all  the  surface  velocities  for  any  set  of  band  drive  voltages.  Once  the  surface  velocities  are 


U 

<6 


D 


F 

~Q 


(65) 
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GEOMETRY  DEPENDENT 


Fig.  7  —  Row  diagram  of  the  ATILA-CHIEF  CWC  model. 

known,  CHIEF  can  then  be  used  to  determine  the  field  pressures.  This  method  is  analogous  to  the 
method  used  with  the  RING-CHIEF  model.  The  model  was  implemented  by  setting  up  a  CHIEF 
model  of  the  CWC  that  matches  the  ATILA  model.  This  CHIEF  model  creates  a  radiation-impedance 
file.  ATILA  is  then  run  (accessing  the  radiation-impedance  file)  with  the  first  driven  band  having 
unity  drive  voltage  applied  to  its  outer  electrode  and  zero  volts  applied  to  the  inner  electrode  and  to  all 
the  other  outer  electrodes.  The  resultant  surface  velocities  output  by  ATILA  are  then  stored  in  a  file. 
This  process  is  repeated  (as  with  the  RING-CHIEF  model)  for  each  of  the  driven  bands  through  band 
20  with  each  set  of  resultant  surface  velocities  appended  to  the  velocity  file. 

In  order  to  give  the  finite-element  model  a  finer  grid  and  give  a  more  conservative  aspect  ratio 
for  the  fairly  thin-walled  tube  (i.e.,  to  reduce  the  ratio  of  the  longest  dimension  of  an  element  to  its 
shortest  dimension),  compatible  CHIEF  and  ATILA  models  of  the  CWC  were  created  where  each  6- 
mm-wide  band  is  divided  into  two  3-mm-wide  regions.  This  CHIEF  model  has  170  surface  areas 
instead  of  the  86  for  the  6-mm  model.  Although  this  significantly  increases  the  computation  time,  this 
CHIEF  run  only  has  to  be  performed  once  for  each  drive  frequency  with  a  specified  set  of  physical 
parameters  to  compute  the  radiation-impedance  matrix.  The  positions  of  the  ATILA  nodes  and  CHIEF 
surfaces  for  these  models  are  given  in  Table  2.  The  ATILA  model  utilizes  axial  symmetry  and  was 
generated  with  the  aid  of  the  MOSAIQUE  Automatic  Mesh  Generator  that  is  part  of  the  ATILA 
package.  The  model  consists  of  84  eight-node,  isoparametric,  quadrilateral  (QUAD08P)  elements  [63]. 
These  elements  are  designed  to  model  any  piezoelectric  material.  The  QUAD08P  elements  run  from 
the  bottom  to  the  CWC  to  the  top  with  the  First  element  having  nodes  1,  2,  3,  171,  257,  256,  255,  and 
170.  The  CHIEF  model  is  coupled  to  the  ATILA  model  by  the  use  of  170  three-node  isoparametric 
(LINE03Z)  elements  [64].  These  elements  are  used  to  couple  the  ATILA  model  with  the  radiation- 
impedance  matrix  output  by  CHIEF  and  are  superimposed  on  the  radiating  surfaces  of  the  mesh  of  the 
solid  structure  (one  for  each  CHIEF  surface  area).  The  LINE03Z  elements  run  counterclockwise 
around  the  surface  with  the  first  element  having  nodes  255,  256,  and  257.  The  ATILA  model  has  a 
total  of  423  nodes,  as  shown  in  Table  2. 

The  output  of  the  ATILA-CHIEF  model  is  average  normal  surface  velocities  for  the  two  end 
surfaces  and  the  168  surfaces  of  the  3-mm  wide  regions.  A  FORTRAN  program  was  written  to 
convert  these  into  velocities  for  the  86  surface  areas  as  used  in  the  RING-CHIEF  model  where  each 
band  surface  velocity  is  the  average  of  the  two  3-mm  regions  spanning  the  band. 
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Table  2  —  Numbering  and  axial  positions  of  ATILA  nodes  and  CHIEF  surface 
borders  for  the  ATILA-CHIEF  model. 


POSITION 

CHIEF 

NODES 

CHIEF 

BAND 

85 

“1 

0.2520 

169 

254 

423 

02505 

86 

168 

422 

84 

02490 

167 

253 

421 

GND 

02475 

87 

166 

420 

83 

02460 

165 

252 

419 

02445 

88 

164 

418 

82 

02430 

163 

251 

417 

40 

02415 

89 

162 

416 

81 

02400 

161 

250 

415 

0.0300 

21 

180 

275 

0.0285 

160 

20 

274 

10 

0.0270 

19 

179 

273 

4 

0.0255 

161 

18 

272 

9 

0.0240 

17 

178 

271 

0.0225 

162 

16 

270 

8 

0.0210 

15 

177 

269 

3 

0.0195 
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CWCMDL  Program 

The  FORTRAN  computer  program  CWCMDL  is  used  to  calculate  sets  of  band  drives  that  can  be 
used  to  generate  evanescent  pressure  fields  within  the  CWC.  The  program  uses  the  method  described 
earlier  to  determine  the  complex  band  drives  from  transfer  matrices  computed  with  either  CHIEF-I, 
RING-CHIEF,  or  ATILA-CHIEF.  The  program  can  also  compute  the  predicted  pressure  field  within 
the  CWC  from  a  specified  set  of  band  drives  using  either  RING-CHIEF  or  ATILA-CHIEF. 

The  parameters  that  are  required  by  CWCMDL  for  each  run  are  contained  in  a  text  file.  The 
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program  reads  the  parameters  from  the  input  file  and  then  performs  the  necessary  computations.  A 
more  detailed  description  of  CWCMDL  is  given  in  Appendix  B. 

EXPERIMENTAL  SETUP 

To  evaluate  the  performance  of  a  PVDF  CWC  a  prototype  was  constructed  with  40  driven  bands. 
An  experimental  setup  was  designed  and  constructed  that  allowed  operation  of  the  CWC  prototype 
with  water  on  the  inside  and  with  either  air  or  castor  oil  on  the  outside.  It  was  desirable  to  investigate 
both  of  these  configurations  because  a  full  sized  CWC  might  be  more  easily  constructed  if  it  were  air 
backed.  However,  if  pressurized  operation  were  required  (i.e.,  to  simulate  conditions  deep 
underwater),  a  full  sized  CWC  might  be  enclosed  in  a  liquid-filled  pressure  vessel.  Independently 
controlled  electronic  drives  were  supplied  for  up  to  40  bands.  An  electronic  measurement  system  was 
constructed  and  software  written  to  accurately  map  the  acoustic  pressure  within  the  CWC. 

CWC  Array 


A  prototype  CWC  was  constructed  using  a  PVDF  tube  as  the  active  element.  The  CWC  consists 
of  (from  the  inside  to  the  outside,  see  Fig.  8);  (1)  a  urethane  coating,  (2)  an  inner  (ground)  electrode, 
(3)  a  PVDF  tube,  and  (4)  the  outer  electrodes. 

A  piezoelectric  PVDF  tube  was  acquired  from  THORN-EMI.  PVDF  is  a  fluorocarbon  polymer 
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that  is  used  in  a  wide  range  of  applications  from  wiring  insulation  to  chemical  storage  containers.  It 
can  be  easily  extruded  into  sheets  and  tubes.  The  PVDF  has  very  little  piezoelectric  activity  when  it  is 
extruded.  The  piezoelectric  properties  are  dramatically  enhanced  by  first  stretching  the  PVDF  and 
then  applying  a  large  electric  field  across  the  sheet  (or  tube  wall)  under  elevated  temperature  [65], 

The  PVDF  tube  was  supplied  with  thin  copper  electrodes  covering  all  of  the  inside  and  outside 
surfaces.  The  tube  has  an  inner  radius  of  1.55  cm,  an  outer  radius  of  1.625  cm,  and  a  length  of  25.2 
cm.  A  tube  of  these  dimensions  was  used  because  it  has  the  largest  radius  and  length  that  could  be 
supplied  at  the  time.  However,  the  manufacturer  feels  that  there  would  be  no  serious  technical 
problems  involved  in  producing  tubes  several  times  this  size  —  provided  that  there  was  sufficient 
incentive  for  the  investment  in  additional  equipment  Insulator  bands  were  then  formed  on  the  outer 
side  of  the  tube  by  machining  away  narrow  strips  of  the  outer  electrode.  It  was  not  necessary  to 
divide  the  inner  electrode  into  separate  bands  because  it  was  used  as  the  common  ground  electrode  for 
all  bands.  This  process  created  40  independent  outer  electrode  bands  with  band  spacing  of  6  mm  and 
an  additional  grounded  band  at  each  end  approximately  6-mm  wide.  The  ground  bands  are  electrically 
connected  to  the  inner  electrode  with  fine  wires  and  conductive  epoxy.  This  provides  for  electrical 
connections  to  the  inner  (ground)  electrode  from  the  outside  of  the  tube. 

Urethane  (Eccobond  45)  end  collars  and  an  approximately  1.5-mm-thick  layer  of  urethane  were 
then  cast  on  the  inside  of  the  array  for  waterproofing.  The  end  collars  are  provided  so  that  the  CWC 
can  be  fitted  snugly  into  a  Lexan  box  (described  later)  and  sealed  to  prevent  leaks.  A  wiring  harness 
is  attached  to  the  outer  bands  of  the  array  with  conductive  epoxy  providing  electrical  connections  from 
the  drive  electronics. 

Drive  Electronics 

The  electrical  band  drives  for  the  CWC  are  supplied  by  a  drive  system  [66]  that  provides  for  40 
channels  of  simultaneous  sinusoidal  drive  with  independently  controlled  amplitudes  and  phases  at  any 
specified  frequency  between  100  and  2000  Hz.  The  system  is  capable  of  providing  drive  levels  up  to 
1000  V  peak-to-peak.  The  hardware  is  controlled  by  software  on  a  VAX  computer.  The  drive  for 
each  channel  is  produced  by  first  loading  the  drive  system’s  internal  memory  with  data  representing 
one  or  more  sine  waves  at  the  desired  amplitude  and  phase  for  each  channel.  Each  channel  has  a 
separate  digital-to-analog  converter  (DAC).  The  sinusoidal  output  is  produced  by  reading  the  sine- 
wave  information  stored  in  memory  and  passing  the  data  through  each  channel’s  DAC.  The  frequency 
at  which  the  DACs  operate  is  adjusted  so  that  a  continuous  sinusoidal  output  is  sustained.  The  DAC 
output  is  filtered  to  remove  the  DAC  frequency  component  and  smooth  the  desired  sinusoidal  output; 
then  if  applied  to  a  power  amplifier  that  develops  the  high  voltage  drive.  This  process  is  multiplexed 
to  drive  all  the  channels  simultaneously.  In  the  above  mentioned  process,  one  sine  wave  is  loaded  in 
memory  for  each  channel  unless  this  would  force  the  system  to  operate  at  a  higher  clock  rate  than 
allowed.  In  this  case  enough  sine  waves  are  loaded  in  memory  to  reduce  the  required  clock  rate  to  an 
acceptable  value. 

Measurement  Setup 

A  measurement  setup  and  rigging  were  designed  that  allow  accurate  measurement  of  the  acoustic 
pressure  within  the  CWC  with  either  air  or  liquid  outside  the  PVDF  array.  The  measurement  setup 
and  water  tank  used  for  the  measurements  are  shown  in  Figs.  9  and  10.  The  measurement  setup 
consists  of:  (1)  a  LEXAN  box  containing  the  CWC,  (2)  a  hydrophone,  (3)  a  computer-controlled, 
single-axis  positioner  and  yoke  assembly,  (4)  a  water  tank,  (5)  a  rigging  assembly  that  mounted  the 
LEXAN  box  and  positioner  in  proper  alignment  on  the  tank. 
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CWC  Mounting  Box 

The  LEXAN  box  that  was  designed  and  constructed  to  mount  the  CWC  is  shown  in  Fig.  11.  It 
is  constructed  from  6.35-mm-thick  LEXAN  sheets  cemented  together  and  is  31.0-cm  high  by  16.5-cm 
wide  by  25.8-cm  long  to  accommodate  the  length  of  the  CWC.  The  box  was  designed  with  holes  in 
the  two  end  walls  that  are  snug  fits  with  the  urethane  collars  at  the  ends  of  the  CWC.  When  the  CWC 
is  inserted  into  the  holes  in  the  box,  and  the  box  is  partially  submerged  in  water,  the  inside  of  the 
CWC  is  flooded  with  water.  The  box  can  then  either  be  left  empty  to  provide  an  air-backed 
configuration  for  the  CWC  or  filled  with  castor  oil  to  provide  a  liquid-backed  configuration.  Castor 
oil  is  used,  instead  of  water,  because  it  is  a  good  insulator  and  therefore  does  not  cause  shorting  across 
the  outer  electrodes.  The  box  is  equipped  with  mounting  hardware  so  that  the  box,  and  therefore  the 
CWC,  can  be  accurately  aligned  with  the  positioner  assembly. 

Hydrophone 

The  hydrophone  used  was  a  Briiel  &  Kjaer  (B&K)  Miniature  Hydrophone  Type  8103.  This 
hydrophone  has  a  cylindrical  piezoelectric  element  with  an  outer  diameter  of  7  mm  and  a  length  of  8 
mm.  A  somewhat  smaller  hydrophone  that  was  available  (Atlantic  Research  LC-5)  was  used  in  some 


29 


LUKER 


Fig.  10  —  Cross-section  view  of  measurement  setup  and  rigging. 

early  measurements  with  the  CWC  but  was  abandoned  in  favor  of  the  B&K  8103  because  the  8103’s 
higher  pressure  sensitivity  gave  more  usable  signals  when  the  CWC  was  operated  with  low  drive 
voltages.  The  low  drive  voltages  were  deemed  prudent  because  the  first  prototype  CWC  was 
destroyed  when  arcing  occurred  across  the  outer  electrodes.  A  B&K  8103  had  also  been  successfully 
used  in  measurements  with  a  planar  wavenumber  calibrator  having  the  same  electrode  spacing  and 
generating  evanescent  fields  in  the  same  frequency  and  wavenumber  regions  that  are  intended  for  the 
prototype  CWC  [19]. 

The  hydrophone  is  attached  to  a  monofilament  line  that  is  strung  tightly  through  the  center  of  the 
CWC  (see  10).  The  line  is  attached  to  a  yoke  assemble  that  in  turn  is  attached  to  the  traveling  head 
on  the  positioner.  This  allows  positioning  of  the  hydrophone  anywhere  within  the  CWC  with  a 
minimum  amount  of  rigging  within  the  CWC.  The  positioner  can  then  automatically  position  the 
hydrophone  anywhere  along  a  line  running  through  the  CWC  (usually  along  the  axis  of  the  CWC). 

Positioner 

The  positioner  is  a  Klinger  Positioning  Stage  with  a  UE72  Stepping  Motor.  It  is  controllable 
either  manually  with  a  Klinger  CC1.1  Programmable  Stepping  Motor  Controller  or  remotely  by 
commands  sent  to  the  controller  via  a  GPIB-488  interface.  The  positioner  and  yoke  assembly  can 
position  the  hydrophone  throughout  the  length  of  the  CWC  and  up  to  1.5  cm  beyond  the  ends  with 
steps  as  small  as  0.0254  mm  (0.001  in.). 

Water  Tank 

The  water  tank  is  made  of  fiberglass  and  is  reinforced  by  a  steel  and  wood  frame.  The  tank 
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Fig.  11  —  LEXAN  CWC  mounting  box  and 
rigging. 


measures  120x60x60-cm  deep  and  is  supported  on  air  bags  to  minimize  coupling  of  structural 
vibrations.  During  measurements  the  tank  is  typically  filled  to  a  depth  of  45  cm  and  the  CWC  is 
submerged  to  a  depth  of  at  least  15  cm. 

Electronic  Measurement  System 

An  electronic  measurement  system  was  constructed  and  software  written  to  accurately  map  the 
pressure  field  (both  amplitude  and  phase)  within  the  CWC. 

Measurement  System  Hardware 

A  block  diagram  of  the  measurement  system  is  shown  in  Fig.  12.  The  system  is  controlled  by  a 
PC-compatible  microcomputer.  The  PC  can  control  the  motion  of  the  positioner,  and  therefore  the 
hydrophone,  by  issuing  commands  over  the  General  Purpose  Interface  Bus  (GPIB)  to  the  positioner 
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Fig.  12  —  Diagram  of  the  measurement  system. 

controller.  The  signal  from  the  hydrophone  is  amplified  by  an  EG&G-PARC  Model  113  preamplifier 
and  then  captured  on  a  Nicolet  Explorer  III  digital  oscilloscope.  Additional  filtering  is  sometimes  used 
between  the  preamplifier  and  the  digital  oscilloscope.  The  other  channel  of  the  oscilloscope  is  used  to 
monitor  a  reference  signal.  The  data  from  the  oscilloscope  is  transferred  over  the  GPIB  to  the  PC  for 
processing.  The  PC  displays  and  records  the  position  of  the  measurement  and  the  amplitude  and 
relative  phase  of  the  signal  from  the  hydrophone.  The  phase  relative  to  the  reference  signal  is 
determined  by  performing  a  single-frequency  discrete-Fourier-transform  [67,68]  on  each  of  the  signals. 

Measurement  System  Software 

A  QuickBASIC  program  called  CWCMAP  was  written  for  the  PC  to  control  the  positioner  and 
oscilloscope  and  to  determine  the  amplitude  and  phase  of  the  signals  captured  by  the  oscilloscope. 

The  program  first  presents  the  user  with  a  menu  (see  Fig.  13)  allowing  the  entry  of  required 
information  such  as  the  first  and  last  positions  of  the  hydrophone  and  the  distance  between  data  points 
for  mapping  the  pressure  field.  The  user  also  enters  the  drive  frequency  and  the  name  of  the  data  file 
that  is  to  receive  the  position  and  pressure  data. 

The  procedure  used  to  map  the  pressure  field  within  the  CWC  is: 

•  Enter  the  appropriate  parameters  (e.g.,  frequency,  band  drive  amplitudes  and  phases, 
etc.)  with  the  drive  system  software  on  the  VAX. 

•  Instruct  the  VAX  to  start  the  drive  system  and  apply  the  appropriate  drives  to  the 
CWC. 

•  Enter  the  appropriate  mapping  parameters  with  the  measurement  system  software  on 
the  PC 

•  Instruct  the  measurement  system  PC  to  start  the  map  of  the  pressure  field. 

During  the  measurement  process  one  channel  of  the  digital  oscilloscope  is  connected  to  the 
output  of  the  preamplifier  while  the  other  channel  is  connected  to  one  of  the  drive  system  outputs  to 
provide  a  phase  reference  to  determine  the  relative  phase  of  the  pressure  signal.  The  measurement 
software  then  maps  the  pressure  field  by  moving  the  positioner  to  the  specified  measurement  positions 
and  at  each  position: 

•  Captures  the  hydrophone  and  reference  signals  on  the  digital  oscilloscope. 
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Cylindrical  Wavenumber  Calibrator 
Pressure-Field  Mapping  Program 
L.D.  Luker  NRL/USRD 


1.  Position  of  origin: 

2.  First  mapping  position: 

3.  Mapping  step  size; 

4.  Last  mapping  position: 

5.  Frequency: 

6.  Data  path: 

7.  Data  filename: 

8.  Annotation: 
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CWCMAP  data 


Press  the  item  number  to  change  or  a  function  key: 


FI  F2  F3  F4  F5  F6 

HELP  FILES  CAL  MATRIX  MAP 


F7  F8 

PLOT 


F9 


FI  0 
EXIT 


Fig.  13  —  Menu  screen  for  the  CWC  measurement  software. 


Transfers  the  data  to  the  PC 

Determines  the  amplitudes  and  phases  of  the  signals. 

Records  the  hydrophone  position  and  amplitude  of  the  hydrophone  signal  and  its 
phase  relative  to  the  reference  signal.These  data  are  then  written  to  the  specified  data 


file. 


EXPERIMENTAL  RESULTS 

The  goal  of  the  research  presented  in  this  thesis  was  to  investigate  the  characteristics  of  the 
PVDF  CWC  prototype.  This  included  the  calculation  of  sets  of  band  drives  that  could  be  used  to 
produce  evanescent  pressure  fields  within  the  CWC.  Several  methods  were  used  to  determine  the 
transfer  matrix  for  the  system.  The  performance  of  the  CWC  was  then  evaluated  for  each  of  the 
different  methods  of  determining  the  transfer  matrix.  This  section  describes  the  measurements  that 
were  made  and  presents  the  results  of  those  measurements. 

Calibration  of  the  Drive  System 

A  technique  was  devised  to  improve  the  accuracy  of  the  electrical  drives  to  the  CWC  bands  by 
correcting  for  errors  in  the  drive  system.  These  errors  exist  because  of  variations  in  the  performance 
of  the  electronic  modules  from  channel  to  channel.  The  bands  are  first  driven  at  uniform  amplitude 
with  a  10°  phase  increment  between  each  succeeding  band.  The  amplitudes  and  phases  of  the  voltages 
at  the  connections  on  the  bands  are  then  measured  and  recorded  using  the  measurement  system.  The 
10°  phase  increment  between  bands  (rather  than  uniform  phase)  is  used  to  assist  in  detecting  swapped 
connections  in  the  drive  system  or  CWC  wiring  harness  and  misdirection  of  drives  by  the  drive-system 
software.  A  FORTRAN  program  was  written  that  will  read  a  set  of  desired  drives  from  a  file,  correct 
the  drives  to  cancel  out  the  measured  drive  system  errors,  and  write  the  corrected  drives  to  a  file  that 
can  then  be  loaded  directly  into  the  drive  system.  Plots  of  a  typical  set  of  drive  system  errors 
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Fig.  14  —  Amplitude  and  phase  errors  at  the  band  connections  for  800  Hz  drive  with 
and  without  the  drive  system  error  corrections. 


measured  directly  at  the  band  connections  at  800  Hz  (without  and  with  the  corrections)  are  shown  in 
Fig.  14.  The  results  show  a  significant  improvement  in  the  accuracy  of  the  band  drives.  The  small 
errors  remaining  after  the  corrections  give  an  indication  of  the  accuracy  of  the  electronic  measurement 
system.  This  system  is  able  to  reproduce  measurements  with  amplitude  errors  within  ±0.4%  and  phase 
errors  within  ±0.3°.  Although  little  variation  has  been  observed  in  the  corrections  required  at  a  given 
frequency  (unless  some  of  the  electronics  are  replaced),  the  calibration  process  is  repeated  before  each 
measurement  session  to  insure  that  all  the  band  drives  are  functioning  properly.  The  drive-correction 
program  is  then  run  to  correct  and  load  each  set  of  desired  drives. 
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Measured-Transfer-Matrix  Method 

One  method  used  for  determining  the  transfer  matrix  for  the  system  is  to  actually  measure  it 
This  is  probably  the  most  accurate  method,  if  done  carefully,  since  the  only  errors  involved  are 
measurement  errors.  Unlike  the  other  methods  presented,  it  does  not  rely  on  a  model  of  the  an 
idealized  CWC  (i.e.,  in  the  models  all  the  bands  are  assumed  to  exactly  match  the  dimensions 
specified,  have  perfectly  uniform  material  properties,  etc.).  The  measured  transfer  matrix  includes  the 
effects  of  any  imperfections  in  the  actual  array.  For  example,  if,  because  of  manufacturing 
irregularities,  some  of  the  bands  have  dimensional  variations  or  slightly  larger  piezolelectric 
coefficients  than  others,  the  effects  of  this  are  included  in  the  measured  transfer  matrix.  Since  any 
such  imperfections  will  effect  the  evanescent  field,  it  is  advantageous  to  have  them  included  in  the 
transfer  matrix. 

Measurement  of  the  Transfer  Matrix 

The  measurement  program  CWCMAP  was  written  with  an  option  to  measure  the  transfer  matrix 
by  measuring  the  pressure  at  specified  field  points  with  each  band  of  the  CWC  driven  alone  in  turn 
and  write  the  results  to  data  files  in  the  form  of  a  matrix.  The  output  is  actually  to  two  files,  with  one 
containing  the  real  parts  and  one  containing  the  imaginary  parts  of  the  matrix.  The  output  from  the 
DAC  monitor  for  the  drive  to  band  number  one  is  always  used  as  the  phase  reference.  Accordingly, 
the  drive  system  is  always  instructed  to  drive  this  band  at  unity  amplitude  and  zero  phase  shift.  The 
input  to  the  power  amplifier  for  band  number  one  is  then  disconnected  whenever  another  band  is  to  be 
driven  alone. 

Unfortunately,  there  is  no  electronic  communication  between  the  measurement  system  (controlled 
by  a  PC)  and  the  drive  system  (controlled  by  a  VAX).  This  means  that  the  measurement  system 
cannot  automatically  set  up  the  drive  system  or  enable  the  drives.  The  measurement  system  PC  must 
instruct  the  user  to  start  the  drive  system  with  a  specified  set  of  drives.  The  user  then  instructs  the 
VAX  to  operate  the  drive  system  with  the  specified  drives.  When  the  drives  are  achieved,  the  user 
indicates  to  the  PC  that  the  drive  system  is  operating  so  that  the  measurement  can  continue. 

The  procedure  followed  in  performing  the  transfer-matrix  measurements  is: 

•  The  required  parameters  (e.g.,  frequency,  locations  of  the  field  points,  output 
filenames,  etc.)  are  input  to  CWCMAP  on  the  measurement  PC  and  the  program  is 
instructed  to  start  the  transfer-matrix  measurements. 

•  The  PC  then  positions  the  hydrophone  at  the  first  field  point  and  instructs  the  user  to 
drive  band  number  one  alone. 

•  The  drive  system  VAX  is  instructed  to  drive  band  number  one  at  a  constant  amplitude 
and  zero  phase  at  the  desired  frequency  and  to  set  all  the  other  drives  to  zero. 

•  The  PC  is  then  instructed  to  continue  and  measures  the  pressure  at  all  the  field  points. 

•  After  completing  the  measurements  with  band  one  driven,  the  PC  writes  the  data  to  the 
files  and  then  instructs  the  user  to  drive  the  next  band  alone. 

•  This  process  is  continued  until  the  field-point  pressures  have  been  measured  and 
recorded  for  every  band  driven  alone. 

Using  this  procedure,  the  transfer  matrix  for  the  CWC  at  800  Hz  (with  air  outside  the  CWC)  was 
measured  for  70  field  points  located  along  the  axis  of  the  CWC  (i.e.,  every  4  mm  along  the  z  axis 
from  z  =  -13.8  to  z  =  +13.8  cm).  Even  though  each  pressure  measurement  only  took  several  seconds, 
setting  up  the  drives  to  each  band  and  making  the  2,800  measurements  took  several  hours.  A  surface 
plot  of  the  amplitude  of  the  measured  transfer  matrix  is  shown  in  Fig.  15.  In  this  figure  the  70  points 
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for  band  number  one  driven  are  along  the  right-front  edge  running  left  to  right  from  z  equals  -13.8  to 
+  13.8  cm.  The  values  for  each  successive  band  driven  are  represented  by  another  line  on  the  surface. 


Figure  16  shows  the  real  part  of  the  pressure  field  displayed  in  Fig.  15.  These  figures  show  that 
standing-wave  patterns  are  created  in  the  pressure  field  within  the  CWC.  The  pressure  drops  to  a 
relatively  small  value  just  outside  of  the  CWC  no  matter  which  band  is  driven.  (Note:  The  CWC 
extends  from  -12.6  to  +12.6  cm.)  This  frequency  (800  Hz)  is  between  the  first  (half-wavelength)  and 
second  (full-wavelength)  water  column  resonances  within  the  CWC.  The  shape  of  the  standing-wave 
pattern  depends  on  which  band  is  driven.  The  pressure  pattern  is  symmetric  within  the  CWC  (i.e.,  the 
pressure  fields  from  the  upper  20  bands  are  mirror  images  of  the  fields  from  the  lower  20  bands). 

This  is  reassuring  since  the  CWC  has  a  symmetric  design.  This  also  indicates  that  there  is  little  effect 
on  the  pressure  field  from  the  B&K  8103  hydrophone  housing.  The  housing  extends  for 
approximately  3  cm  only  on  the  left  side  (i.e.,  negative  z)  of  the  hydrophone’s  active  element  If  the 
housing  was  having  a  significant  effect  on  the  pressure  field,  this  effect  would  tend  to  show 
differences  when  measuring  pressures  near  the  left  end  of  the  CWC  (where  most  of  the  housing  is 
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Fig.  16  —  Real  part  of  the  measured  transfer  matrix  at  800  Hz. 


outside  the  CWC)  and  pressures  near  the  right  end  of  the  CWC  (where  all  of  the  housing  is  well 
within  the  CWC). 

Pressure  Fields  from  Calculated  Drives 

The  measured  transfer  matrix  was  used  to  calculate  a  set  of  drives  to  produce  an  evanescent  field 
within  the  CWC.  This  was  accomplished  with  the  aid  of  Mathcad.  Mathcad  [69]  is  a  computer 
program  that  allows  the  user  to  perform  complex  mathematical  computation  and  plotting.  A  Mathcad 
document  was  created  to  calculate  CWC  drives  using  the  method  described  earlier  and  to  calculate  and 
plot  the  predicted  field  within  the  CWC  if  these  drives  were  used. 

The  measured  transfer  matrix  values  for  the  46  positions  running  from  -9.0  to  +9.0  cm  were  used 
to  calculate  a  set  of  band  drives  that  should  produce  a  60  m/s  evanescent  field  over  this  region.  The 
predicted  field  for  these  drives  is  shown  in  Fig.  17.  In  this  figure  the  thin  line  indicates  the  extent  of 
the  window  (the  region  in  which  the  evanescent  field  was  being  created).  Any  deviation  from  the 
desired  level  (0  dB)  is  considered  an  error. 
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Fig.  17  —  Predicted  amplitude  errors  for  a  60-m/s  evanescent  pressure  field  at  800  Hz. 

The  measured  transfer-matrix  for  46  positions  running  from  -9.0  to  +9.0  cm  was  used. 

The  CWC  was  then  driven  with  the  calculated  drives  and  the  pressure  field  mapped.  The  results 
are  shown  in  Fig.  18.  Since  this  graph  format  will  be  used  extensively  to  present  plots  of  the  pressure 
field  within  the  CWC,  it  will  be  described  in  some  detail  here.  The  z  axis  coincides  with  the 
centerline  of  the  CWC  with  the  center  of  the  CWC  at  z  =  0.  The  upper  graph  plots  the  amplitude  in 
relative  decibels  and  the  lower  graph  plots  the  phase  in  degrees.  The  predicted  amplitude  and  phase 
are  shown  with  the  thick  lines.  The  thin  line  at  0  dB  again  indicates  the  window  extent.  This  thin 
line  is  almost  completely  obscured  by  the  thick  line  representing  the  predicted  pressure  in  this 
particular  plot.  The  shaded  dots  indicate  the  actual  pressures  measured  along  the  centerline  of  the 
CWC  (the  z  axis). 

The  measured  pressure  field  in  Fig.  18  shows  the  large  variations  of  amplitude  (18  dB)  and  phase 
that  are  characteristic  of  a  standing  wave  field  rather  than  the  desired  evanescent  field.  An 
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Fig.  18  —  Pressures  within  the  CWC  for  the  60-m/s  evanescent  field  at  800  Hz.  The 
measured  transfer-matrix  for  46  positions  running  from  -9.0  to  +9.0  cm  was  used. 


examination  of  the  drives  used  (see  Fig.  19)  shows  large  differences  of  amplitude  and  phase  between 
adjacent  bands.  Although  solutions  of  this  form  can  produce  good  predicted  results  when  inserted  into 
a  mathematical  model  because  of  the  high  precision  available  in  the  computations,  they  place  severe 
demands  on  the  precision  of  the  electronic  drive  system.  This  can  be  demonstrated  if  one  considers  a 
hypothetical  case  where  it  is  desired  to  produce  unity  pressure  at  a  point  in  space  that  is  equally 
effected  by  the  drives  to  only  two  bands.  This  could  be  done  by  driving  both  bands  with  drives  that 
would  produce  a  pressure  of  +0.5  at  the  point  if  each  band  were  driven  alone.  The  desired  pressure 
could  also  be  achieved  by  driving  the  first  band  such  that  its  drive  alone  would  produce  a  pressure  of 
+  11.0  at  the  point  and  driving  the  second  band  such  that  its  drive  alone  would  produce  a  pressure  of  - 
10.0  at  the  point.  Mathematically,  both  of  these  methods  are  equally  acceptable.  However,  if  the  true 
pressure  produced  at  the  point  is  1.0%  larger  in  magnitude  than  expected  when  the  second  band  is 
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Fig.  19  —  Amplitudes  and  phases  of  the  band  drives  used  to  produce  the  pressure 

field  shown  in  Fig.  18. 


driven  alone,  the  total  pressure  in  the  first  case  (+1.005)  will  have  an  error  of  only  +0.5%  while  the 
total  pressure  in  the  second  case  (+0.9)  will  have  an  error  of  -10.0%.  Even  with  no  errors,  the  second 
case  requires  much  larger  drives  to  achieve  the  same  pressure.  For  these  reasons  it  is  desirable  to  have 
drives  that  are  as  uniform  as  possible  in  amplitude  and  phase. 

Pressure  Fields  from  Drives  Calculated  with  Non-Zero  a 

To  compute  more  suitable  drives,  the  technique  described  earlier  was  utilized  by  adding  one- 
thousandth  of  the  average  absolute  amplitude  of  the  elements  of  M'  to  the  diagonal  elements. 

Although  this  shifts  the  solution  from  the  ideal  least-squares  solution,  it  can  significantly  improve  the 
uniformity  of  the  drives.  These  drives  (see  Fig.  20)  show  more  uniform  amplitudes  and  essentially 
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Fig.  20  —  Drives  calculated  by  setting  a  equal  to  one-thousandth  the  average  absolute 

amplitude  of  the  M'  elements. 


uniform  phase  steps  between  interior  bands.  The  predicted  field  for  these  drives  is  shown  in  Fig.  21. 
This  shows  the  degradation  of  the  predicted  pressure  field  incurred  with  this  a.  Although  it  is 
debatable  just  how  uniform  the  evanescent  field  within  an  operational  CWC  should  be,  a  field  that  is 
uniform  to  within  ±1.5  dB  would  certainly  give  useful  information. 

The  CWC  was  then  operated  with  these  drives  and  the  resultant  pressure  field  mapped.  The 
results  (see  Fig.  22)  show  fairly  good  agreement  with  the  desired  and  predicted  evanescent  field  (±1.7 
dB,  ±12°)  throughout  the  window  region. 
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Fig.  21  —  Predicted  amplitude  errors  for  the  drives  shown  in  Fig.  20. 
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CHIEF  Model  with  Rigid  Walls 

To  first  examine  the  CWC  with  CHIEF,  it  was  modeled  as  though  it  was  a  rigid  tube  totally 
emersed  in  water  with  only  band  1  driven.  This  was  done  using  the  CHIEF  model  with  all  the  normal 
surface  velocities  specified.  The  CWC  was  modeled  with  all  the  surface  velocities  set  to  zero  except 
for  the  surfaces  of  band  1.  The  inner  and  outer  surfaces  of  band  1  were  given  normal  velocities  of 
-1.0  and  1.0,  respectively  (in  CHIEF,  normal  surface  velocities  are  always  defined  with  the  normal 
pointing  outward  from  the  radiating  surface).  The  pressure  within  the  CWC  was  then  computed  for  31 
points  spaced  1.0-cm  apart  along  the  axis  of  the  tube  for  a  2772  Hz  drive.  This  frequency  was  chosen 
because  it  is  the  first  resonant  frequency  for  an  open-ended  rigid  tube  with  these  dimensions  emersed 
in  water.  The  results  (see  Fig.  23)  showed  a  sharp  pressure  peak  opposite  the  driven  band  and  no 
signs  of  the  water  column  resonance  that  should  be  there.  The  position  of  band  1  is  indicated  by  the 
short  line  just  above  the  z  axis. 
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Fig.  22  —  Predicted  and  measured  pressures  within  the  CWC  using  the  drives  shown 

in  Fig.  20. 


To  determine  if  the  integration  routines  used  by  CHIEF  were  giving  erroneous  results  because  of 
the  relatively  thin  wall  of  the  CWC  compared  to  the  wavelength  of  sound  at  this  frequency  (see 
Appendix  C),  the  dimensions  of  the  CWC  in  the  CHIEF  model  were  changed.  The  inner  radius  and 
length  of  the  CWC  were  kept  the  same  but  the  outer  radius  was  changed  to  2.15  cm.  This  increased 
the  wall  thickness  from  0.075  to  6.00  mm.  All  the  other  parameters  for  the  CHIEF  model  remained 
unchanged.  The  CHIEF  model  was  then  run  with  the  new  dimensions  and  the  pressures  within  the 
CWC  were  computed.  The  results  (see  Fig.  24)  show  the  eigenfunction  one  would  expect  near  the 
first  resonance  in  the  water  column.  This  confirmed  that  the  problem  with  the  CHIEF  model  using  the 
correct  dimensions  was  caused  by  the  relatively  thin  wall.  Unfortunately,  this  limitation  with  CHIEF 
is  not  mentioned  in  the  CHIEF  Users  Manual  [22]. 
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Fig.  23  —  Predicted  pressures  along  the  centerline  of  a  rigid  CWC  with  a  0.75-mm 
wall.  Pressures  were  computed  using  CHIEF  with  band  1  having  unity  velocity  at 

2772  Hz. 


It  is  believed  that  the  problem  with  CHIEF  is  from  errors  introduced  when  evaluating  the  surface 
integrals  for  points  that  are  very  near  the  integration  surface.  In  particular,  that  errors  are  encountered 
when  numerically  evaluating  the  surface  integral  for  a  point  that  is  on  the  "back"  side  of  the 
integration  surface.  To  make  CHIEF  usable  for  the  0.75-mm-wall  CWC,  the  program  was  modified  to 
perform  these  integrals  using  the  same  approach  that  CHIEF  uses  for  its  self  integrals.  This  essentially 
consists  of  dividing  the  surface  to  be  integrated  into  four  triangular  surfaces  that  meet  at  a  point  in  the 
center  of  the  original  surface  and  changing  the  integration  variable.  A  more  detailed  description  of  the 
modifications  made  to  CHIEF  to  model  the  CWC  are  given  in  Appendix  C.  For  this  modification  to 
CHIEF  to  work,  the  front  and  back  of  the  thin  part  of  the  radiating  body  must  both  be  divided  into  the 
same  arrangement  of  subregions.  To  put  it  more  precisely,  the  front  and  back  surfaces  must  be 
divided  in  such  a  way  that  (in  CHIEF’S  «,v  coordinates)  the  center  of  the  back  subregion  approaches 
the  center  of  the  front  subregion  as  the  thickness  of  the  thin  body  approaches  zero.  Although  this 
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Fig.  24  —  Predicted  pressures  along  the  centerline  of  a  rigid  CWC  with  a  6.00-nun 
wall.  Pressures  were  computed  using  CHIEF  with  band  1  having  unity  velocity  at 

2772  Hz. 


requirement  limits  the  general  application  of  the  technique,  the  requirement  is  satisfied  for  the  CWC 
models. 

The  modified  CHIEF  model  was  then  run  with  the  CWC’s  actual  dimensions  (0.75-mm  wall)  and 
the  pressures  within  the  CWC  were  computed.  Figure  25  shows  a  comparison  of  these  pressures  with 
the  previously  shown  pressures  from  the  unmodified  CHIEF  model  of  a  6.00-mm-wall  CWC.  The 
results  show  that  the  modified  CHIEF  program  gives  excellent  agreement  with  the  6.00-mm-wall 
model.  The  only  differences  one  would  expect  between  0.75-mm-wall  and  6.00-mm-wall  rigid  tubes 
would  be  slightly  different  end  corrections  and  a  small  difference  in  coupling  at  the  ends  between  the 
water  outside  and  inside  the  tube.  The  modified  CHIEF  program  was  used  for  all  further  CHIEF 
models  of  the  0.75-mm-wall  CWC. 
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Fig.  25  —  Predicted  pressures  at  2772  Hz  from  standard  CHIEF  with  a  6.00-mm  wall 
and  CHIEF  modified  for  the  thin  wall  with  a  0.75-mm  wall. 


CHIEF-I  Model 

The  CWC  was  modeled  using  the  modified-CHIEF  model  with  the  interior-wall  impedance 
specified  (CHIEF-I).  This  model  was  described  earlier. 

Comparison  with  a  Single  Band  Driven 

The  CHIEF-I  model  was  used  to  calculate  the  predicted  pressure  within  the  CWC  if  only  band  4 
were  driven.  The  pressure  within  the  CWC  was  then  measured  with  only  band  4  driven.  The  CWC 
was  air-backed  for  this  measurement  Figure  26  shows  a  comparison  of  the  modeled  and  measured 
pressures  along  the  centerline  of  the  CWC  with  band  4  driven  alone  at  800  Hz.  The  position  of  band 
4  (from  z  =  -10.2  cm  to  z  =  -9.6  cm)  is  indicated  by  the  short  line  just  above  the  z  axis.  It  can  be 
seen  that  there  is  fairly  good  agreement  between  the  CHIEF-I  model  and  the  measured  pressures  for 
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Fig.  26  —  Predicted  and  measured  pressures  within  the  air-backed  CWC  with  only 
band  4  driven  at  800  Hz.  The  predicted  pressure  was  calculated  using  the  CHIEF-I 

model. 


this  simple  model.  The  shape  of  the  pressure  Field  standing  wave  pattern  is  determined  by  how  much 
of  the  energy  supplied  by  the  driven  band  is  coupled  into  each  of  the  allowed  modes  within  the  CWC. 
This  pattern  is  different  for  each  driven  band  because  the  bands  are  located  at  different  positions  along 
the  CWC  and  therefore  excite  the  modes  in  different  proportions  (see  Figs.  15  and  16). 

RING-CHIEF  Model 

The  CWC  was  modeled  using  the  RING-CHIEF  model.  This  model  considers  the  bands  of  the 
CWC  thin-walled  piezoelectric  rings  and  includes  the  effects  of  the  mutual  radiation  impedances.  The 
radiation  impedance  matrix  for  this  model  was  calculated  using  CHIEF. 
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Comparison  with  a  Single  Band  Driven 

The  RING-CHIEF  model  was  used  to  calculate  the  predicted  pressure  within  the  CWC  if  only 
band  4  were  driven.  Since  the  radiation  impedance  matrix  computed  by  CHIEF  assumes  that  the 
CWC  is  filled  with  and  surrounded  by  water,  the  LEXAN  mounting  box  was  filled  with  castor  oil. 

The  pressure  within  the  CWC  was  then  measured  with  only  band  4  driven.  Figure  27  shows  a 
comparison  of  the  modeled  and  measured  pressures  along  the  centerline  of  the  CWC  with  band  4 
driven  alone  at  800  Hz.  It  can  be  seen  that  there  is  fairly  good  agreement  between  the  RING-CHIEF 
model  and  the  measured  pressures  for  this  model. 

ATI  LA-CHIEF  Model 

In  order  to  obtain  a  more  complete  model  of  the  CWC,  the  ATILA-CHIEF  model  was  used. 
Because  of  problems  encountered  in  the  first  attempt  to  model  the  CWC  with  the  ATILA-CHIEF 
combination,  an  effort  was  undertaken  to  measure  the  wall  velocities  of  the  CWC  so  that  these  could 
be  compared  with  the  predicted  velocities  from  ATILA.  This  would  allow  the  evaluation  of  the 
performance  of  the  ATILA  model  with  and  without  the  fluid  loading  provided  by  CHIEF  but  without 
the  added  complication  of  then  using  CHIEF  to  calculate  field  pressures  from  the  ATILA  velocities. 

Radiation-Impedance  Matrix 

The  ATILA-CHIEF  model,  like  the  RING-CHIEF  model,  requires  that  CHIEF  be  used  first  to 
determine  the  radiation-impedance  matrix  for  the  radiating  surface.  The  radiation  impedances  must  be 
calculated  by  CHIEF  using  a  surface  division  scheme  that  matches  the  nodal  pattern  for  ATILA.  This 
was  accomplished  using  the  170  CHIEF  surface  divisions  given  in  Table  2.  A  surface  plot  of  the 
imaginary  part  of  the  radiation-impedance  matrix  is  shown  in  Fig.  28.  The  imaginary  part  is  about  an 
order  of  magnitude  larger  that  the  real  part.  Only  every  fourth  element  of  every  fourth  row  is  plotted 
to  improve  clarity.  The  surfaces  run  up  the  outside  of  the  CWC  then  back  down  the  inside  (see  Table 
2).  The  diagonal  elements  are  the  self  impedances.  The  plot  shows  that  there  is  very  little  interaction 
between  the  inside  and  the  outside  of  the  CWC.  It  also  shows  that  there  is  a  relatively  large  amount 
of  interaction  between  elements  inside  the  CWC.  This  would  be  expected  since  the  CWC  is 
essentially  a  duct.  The  shapes  of  the  curves  within  the  CWC  are  similar  to  the  pressure  field  within  a 
rigid  CWC  for  each  band  driven  alone  at  this  frequency. 

Figure  29  shows  the  erroneous  radiation-impedance  matrix  computed  with  the  unmodified  CHIEF 
program  (i.e.,  without  the  thin-wall  modifications  described  in  Appendix  C)  for  this  same  model.  This 
figure  shows  large  interactions  between  the  inside  and  outside  of  each  band  that  fade  away  quickly  as 
the  distance  from  the  band  increases.  It  should  be  noted  that  this  model  divides  the  surface  into  finer 
divisions  (i.e.,  3-mm-wide  subregions  rather  than  6-mm-wide  subregions)  than  were  used  in  the 
previous  CHIEF  models,  but  the  standard  CHIEF  program  still  gives  errors  with  the  0.75-mm  wall. 
These  results  again  demonstrate  problems  that  can  occur  using  the  standard  CHIEF  program  with 
relatively  thin  vibrating  bodies. 

Predicted  and  Measured  Surface  Velocities 

Different  methods  were  considered  for  measuring  the  surface  velocity  of  the  CWC.  An 
accelerometer  was  first  considered  but  was  rejected  for  the  following  reasons: 

•  The  CWC  wall  is  thin  and  made  of  fairly  light  material.  The  motion  of  the  wall  could 
easily  be  effected  by  the  additional  loading  of  an  accelerometer. 

•  The  accelerations  anticipated  were  smaller  than  could  be  effectively  measured  with 
small  accelerometers.  Larger  accelerometers  might  be  sensitive  enough,  but  this  would 
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Fig.  27  —  Predicted  and  measured  pressures  within  the  oil-backed  CWC  with  only 
band  4  driven  at  800  Hz.  The  predicted  pressure  was  calculated  using  the  RING- 

CHIEF  model. 

aggravate  the  loading  problem. 

•  It  was  desirable  to  measure  the  velocity  at  numerous  places  along  the  surface  of  the 
CWC.  This  would  either  require  numerous  accelerometers  or  removal  and 
reattachment  of  a  single  accelerometer. 

To  avoid  these  problems,  it  was  decided  to  measure  the  surface  velocity  with  the  Dantec  55X 
Laser  Doppler  Vibrometer  (LDV)  available  at  USRD  [70].  This  system  measures  the  surface  velocity 
by  detecting  the  Doppler  shift  of  laser  light  reflected  or  scattered  by  the  surface.  For  a  light  beam 
impinging  on  and  reflected  back  in  a  direction  normal  to  the  reflecting  surface,  the  frequency  of  the 
light  is  shifted  upon  reflection  according  to  the  Doppler  equation  [71] 
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Fig.  28  —  Surface  plot  of  the  imaginary  part  of  the  radiation-impedance  matrix 
computed  with  CHIEF  (modified  to  handle  the  thin-walled  CWC)  at  800  Hz. 


A/  =  21  ,  (6T> 

where  A/  is  the  Doppler  frequency  shift,  v  is  the  normal  velocity  of  the  surface,  and  X  is  the 
wavelength  of  the  light.  The  LDV  system  used  has  a  He-Ne  laser  with  a  wavelength  of  632.8  nm. 

The  optical  system  detects  the  Doppler  shift  of  the  reflected  light  by  an  optical  heterodyne  process, 
where  the  received  light  beam  is  mixed  with  a  reference  beam.  Using  an  LDV  to  measure  the  normal 
surface  velocity  did  not  load  the  surface  at  all  and  could  measure  the  surface  velocity  at  numerous 
places  along  the  outer  wall  of  the  CWC  by  just  moving  the  LDV  beam  along  the  wall. 

The  LDV  was  placed  so  that  it  could  be  translated  parallel  to  the  axis  of  the  CWC.  With  this 
arrangement,  the  beam  could  be  aimed  at  any  point  along  a  line  on  the  outer  wall  of  the  CWC  running 
parallel  to  the  axis  while  keeping  the  beam  perpendicular  to  the  CWC  surface.  This  was  required 
because  a  measurement  of  the  normal  velocity  was  desired  and  the  LDV  indicates  the  velocity 
component  parallel  to  the  beam.  Although  this  arrangement  limited  surface  velocity  measurements  to 
the  outer  wall  of  the  CWC,  it  was  felt  that  these  measurements  could  verify  correct  operation  of  at 
least  that  much  of  the  ATILA-CHIEF  model.  The  model  predicts  that  the  outer  and  inner  wall 
surfaces  of  each  band  have  similar  velocities,  which  is  characteristic  of  thin  walled  tubes.  The  CWC 
was  removed  from  the  LEXAN  box  and  suspended  by  its  leads  in  air.  Band  4  was  driven  at  1000  Hz 
while  all  the  other  bands  were  shorted  to  ground.  Figure  30  shows  the  amplitude  of  the  outer-wall 
surface  velocity  as  measured  by  the  LDV.  Also  shown  in  the  figure  is  the  surface  velocity  predicted 
by  ATILA-CHIEF.  For  this  ATILA  run  the  CWC  was  modeled  operating  in  a  vacuum  by  setting  all 
the  radiation-impedance  matrix  values  equal  to  zero.  This  is  a  reasonable  approximation  of  the  CWC 
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Fig.  29  —  Surface  plot  of  the  imaginary  part  of  the  erroneous  radiation-impedance 
matrix  computed  with  CHIEF  (without  the  thin-wall  modifications)  at  800  Hz. 


operating  in  air  because  air  provides  negligible  loading  for  the  CWC.  This  technique  had  the 
advantage  of  removing  CHIEF  from  the  analysis.  Although  the  signal-to-noise  ratio  of  the  LDV 
measurements  is  fairly  low,  it  can  be  seen  that  the  only  surface  velocity  that  stands  out  from  the  pack 
is  on  the  surface  of  the  driven  band,  as  predicted. 

The  CWC  was  then  suspended  by  its  leads  in  a  tank  of  castor  oil  so  that  it  was  filled  and 
surrounded  by  the  oil.  Oil  rather  than  water  was  used  in  the  tank,  because  the  CWC  was  only  lightly 
coated  with  a  sealer  on  the  outside.  It  was  felt  that  water  would  present  too  much  danger  of  shorting 
across  the  outer  electrode  bands.  Castor  oil  was  chosen  because  it  is  a  good  insulator  and  has  a 
density  and  sound  speed  close  to  water.  The  tank  had  clear  glass  walls.  The  LDV  was  set  up  as 
before  with  the  LDV  laser  beam  passing  through  the  tank  wall  and  castor  oil  to  measure  the  CWC 
wall  velocity.  The  CWC  was  driven  at  1000  Hz,  as  before,  and  the  wall  velocity  was  measured. 

Figure  31  shows  the  measured  and  modeled  CWC  outer-wall  surface  velocities  for  this  setup.  The 
results  show  the  velocity  of  the  driven  band  riding  on  a  wall  velocity  profile  that  coincides  with  the 
pressure  field  in  the  fluid  column  within  the  CWC,  as  would  be  expected  for  a  tube  with  fairly  flexible 
walls.  The  second  fluid  column  resonance  is  near  1000  Hz.  This  also  indicates  why  a  frequency  of 
1000  Hz  was  used  for  the  LDV  measurements.  As  the  frequency  moves  away  from  a  resonance,  the 
pressure  amplitude  within  the  CWC  drops  rapidly.  This  causes  the  effect  of  the  wall  velocity  to  drop 
into  the  background  noise  of  the  LDV  velocity  measurement.  The  LDV  results  indicate  that  the 
ATILA-CHIEF  velocity  profile  prediction  shows  fairly  good  agreement  with  the  measured  velocities. 

Comparison  with  a  Single  Band  Driven 

In  order  to  utilize  the  ATILA-CHIEF  model,  the  CWC  was  again  operated  with  castor  oil  rather 
than  air,  inside  the  LEXAN  mounting  box.  This  was  done  because  the  ATILA-CHIEF  model  also 
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Fig.  30  —  Measured  and  modeled  CWC  outer-wall  surface  velocities  with  the  CWC 
operated  in  air.  Band  4  was  driven  at  1000  Hz  with  all  other  bands  shorted  to  ground. 

assumes  that  the  CWC  is  totally  emersed  in  the  liquid.  The  CWC  was  then  operated  with  only  band  4 
driven  at  800  Hz  and  the  pressure  along  the  centerline  was  measured.  Figure  32  shows  the  measured 
pressure  and  the  pressure  predicted  from  the  ATILA-CHIEF  model.  The  pressure  predicted  from  the 
RING-CHIEF  model  with  band  4  driven  alone  is  also  plotted  again.  This  graph  show  that  there  is 
little  difference  between  the  pressures  predicted  by  RING-CHIEF  and  ATILA-CHIEF  at  this 
frequency.  It  also  shows  that  there  is  fairly  good  agreement  between  both  models  and  the  measured 
pressure. 

RING-CHIEF  Model  with  Modified  Radiation-Impedance  Matrix 

The  three  different  models  of  the  CWC  that  have  been  presented  consist  of: 

•  CHIEF-I  —  Specifies  the  interior-wall  impedance  based  on  a  simple  vacuum-backed  tube 
equation  and  specifies  that  the  outer-wall  surfaces  have  zero  velocity. 
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Fig.  31  —  Measured  and  modeled  CWC  outer-wall  surface  velocities  with  the  CWC 
operated  in  castor  oil.  Band  4  was  driven  at  1000  Hz  with  all  other  bands  shorted  to 

ground. 


•  RING-CHIEF  —  Models  the  CWC  bands  as  thin-walled  piezoelectric  rings  and  includes  the 
effects  of  the  fluid  loading  with  the  radiation  impedance  matrix. 

•  ATILA-CHIEF  —  Uses  a  finite-element  of  the  piezoelectric  PVDF  tube  and  includes  the 
effects  of  the  fluid  loading  with  the  radiation  impedance  matrix. 

All  of  these  models  assume  that  the  CWC  is  filled  with  and  surrounded  by  the  same  fluid.  However, 
the  CHIEF-I  model  was  used  to  predict  the  pressures  within  the  CWC  when  it  is  air  backed.  The 
assumption  was  made  that  the  surface  velocities  on  the  outside  wall  of  the  CWC  had  little  effect  on 
the  pressure  within  the  CWC.  The  RING-CHIEF  model  was  used  to  confirm  this  assumption.  Since 
this  model  includes  the  effects  of  the  outer-surface  velocity  through  the  radiation-impedance  matrix,  it 
is  possible  to  remove  these  effects  by  zeroing  the  appropriate  elements  of  the  matrix.  The  radiation- 
impedance  was  modified  by  zeroing  all  elements  that  related  to  surfaces  outside  of  the  CWC. 
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Fig.  32  —  Pressures  within  the  oil-backed  CWC  with  only  band  4  driven  at  800  Hz 
using  the  RING-CHIEF  model  and  the  ATILA-CHIEF  models. 

Referring  to  Fig.  28,  this  would  be  all  elements  except  the  elements  in  the  rear-most  quadrant  of  the 
surface. 

The  RING-CHIEF  model  was  used  to  predict  the  pressure  within  the  CWC  using  the  modified 
radiation-impedance  matrix  at  800  Hz  with  band  4  driven  alone.  Figure  33  shows  a  comparison  of  the 
predicted  pressures  using  the  full  matrix  and  the  partially  zeroed  matrix.  The  measured  pressure 
within  the  air-backed  CWC  is  also  shown  (see  Fig.  26).  The  results  do  show  some  difference  between 
the  pressures  predicted  using  the  full  matrix  and  the  modified  matrix.  However,  the  difference  is  small 
and  is  comparable  to  the  differences  previously  noted  between  the  different  models.  Since  the  results 
with  the  partially  zeroed  radiation-impedance  matrix  do  show  somewhat  better  agreement  with  the 
measured  pressures  in  the  air-backed  CWC,  this  technique  is  used  to  predict  the  pressures  for  the  air- 
backed  configuration. 
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Fig.  33  —  Predicted  pressures  within  the  CWC  from  RING-CHIEF  with  band  4  driven 
at  800  Hz.  [Z]  was  modified  by  zeroing  elements  involving  outer  surfaces. 

Comparison  of  Evanescent-Field  Results 

To  evaluate  the  performance  of  the  CWC  prototype  as  an  evanescent-wave  generator,  drives  were 
computed  to  produce  evanescent  waves  within  the  CWC.  Various  phase  speed  and  window  width 
combinations  were  investigated.  The  pressures  fields  resulting  from  drives  computed  using  measured 
transfer  matrices  and  modeled  transfer  matrices  are  now  presented  and  compared. 

800  Hz  with  Air  Outside  the  CWC 

To  evaluate  the  air-backed  configuration,  the  CWC  was  operated  at  800  Hz  with  air  inside  the 
LEXAN  box.  Drives  designed  to  produce  a  57  m/s  evanescent  wave  were  calculated  using  the  42 
measured  transfer-matrix  points  running  from  -8.2  to  +8.2  cm.  The  CWC  was  then  operated  with 
these  drives  and  the  pressure  field  within  the  CWC  mapped.  Figure  34  shows  the  predicted  and 
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Fig.  34  —  Pressures  within  the  air-backed  CWC  for  a  57-m/s  evanescent  field  at  800 
Hz.  The  measured  transfer-matrix  for  42  positions  running  from  -8.2  to  +8.2  cm  was 

used. 


measured  pressure  fields  within  the  CWC.  Drives  were  then  computed  for  the  same  scenario  using  the 
CHIEF-I  model  to  determine  the  transfer  matrix  rather  than  measured  values.  Figure  35  shows  the 
measured  pressure  field  within  the  CWC  for  this  set  of  drives.  Although  the  measured  phase  does 
show  a  trend  toward  the  desired  phase,  the  phase  error  is  quite  large.  The  amplitude  error  is  certainly 
much  too  large  for  use  as  a  calibrator.  It  is  believed  that  the  demands  placed  on  the  accuracy  of  the 
model  are  more  than  the  model  is  capable  of  producing  and  that  the  deviations  between  the  modeled 
pressures  and  actual  pressures  are  causing  the  poor  performance. 

It  can  be  observed  that  the  measured  pressure  field  within  the  CWC  is  not  dropping  in  amplitude 
outside  of  the  window  but  is  remaining  large  until  forced  to  drop  by  the  low  impedances  encountered 
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Fig.  35  —  Pressures  within  the  air-backed  CWC  for  a  57-m/s  evanescent  field  at  800 
Hz.  CHIEF-I  pressure  values  for  42  positions  running  from  -8.2  to  +8.2  cm  were 

used. 

at  (or,  more  precisely,  just  beyond)  the  ends  of  the  CWC  (see  Figs.  34  and  35).  This  situation  allows 
the  large  reflections  from  the  ends  of  the  CWC  to  have  a  major  effect  on  the  pressure  field  within  the 
CWC.  This  places  stringent  demands  on  the  accuracy  of  the  pressures  predicted  by  the  model  An 
effort  was  undertaken  to  minimize  these  reflections  and,  thereby,  reduce  the  requirements  on  the 
accuracy  of  the  model. 

It  was  felt  that  if  drives  could  be  found  that  reduced  the  pressure  to  a  relatively  small  value  in 
the  regions  between  the  window  and  the  ends  of  the  CWC  there  would  be  little  energy  left  to  reflect 
from  the  ends.  To  accomplish  this,  all  70  points  of  transfer  matrix  were  used  to  compute  the  drives. 
It  should  be  noted  that  previously  only  transfer  matrix  points  that  were  inside  the  window  were  used 
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and  no  demands  were  placed  on  the  solution  outside  of  the  window.  The  new  desired  pressure  field  P 
consists  of  pressures  with  amplitudes  and  phases  matching  the  evanescent  field  inside  the  window  and 
zero  amplitudes  for  the  points  outside  the  window. 

Drives  were  computed  using  this  technique  and  the  measured  transfer  matrix.  Figure  36  shows 
the  predicted  and  measured  pressure  fields  within  the  CWC  for  this  set  of  drives.  The  plot  shows  that 
the  pressure  is  indeed  dramatically  reduced  before  the  ends  of  the  CWC.  It  should  be  noted  that  both 
the  predicted  and  measured  pressures  show  a  decrease  in  uniformity  inside  the  window.  This  is  not 
surprising  since  the  new  technique  is  a  tradeoff  in  that  the  40  complex  drives  to  the  bands  now  have  to 
also  try  to  drive  the  field  to  zero  outside  of  the  window. 

Drives  were  then  computed  specifying  the  same  desired  pressure  field  but  using  the  CHIEF-I 
model  to  determine  the  transfer  matrix  rather  than  measured  values.  Figure  37  shows  the  measured 
pressure  Field  within  the  CWC  for  this  set  of  drives.  It  can  be  seen  that  the  pressures  field  outside  the 
window  is  driven  down  considerably,  and  that  there  is  a  dramatic  improvement  in  the  uniformity  of 
the  pressure  field  inside  the  window.  This  indicates  that  reducing  the  pressure  outside  the  window 
does  make  it  easier  to  achieve  the  desired  pressures  inside  the  window  using  the  model. 

Another  advantage  of  reducing  the  pressure  outside  the  window  becomes  apparent  when  one 
considers  the  use  of  a  full  scale  CWC  in  actual  operation  as  a  calibrator.  If  a  device  was  placed  in  the 
CWC  and  an  evanescent  field  created  in  a  window  within  the  CWC,  it  would  be  more  desirable  to 
have  little  or  no  pressure  field  outside  of  the  window  to  effect  other  portions  of  the  device.  For  these 
reasons  it  was  decided  to  use  the  technique  of  reducing  the  pressure  outside  the  window  for  the  further 
evaluation  of  the  CWC. 

Drives  were  then  computed  specifying  the  same  desired  pressure  field  but  using  the  RING- 
CHIEF  model  to  determine  the  transfer  matrix.  Since  this  is  an  air-backed  case,  the  partially-zeroed 
radiation-impedance  matrix  was  used.  Figure  38  shows  the  measured  pressure  field  within  the  CWC 
for  this  set  of  drives.  The  resultant  pressure  field  is  very  similar  to  that  obtained  using  CHIEF-I.  This 
is  not  surprising  since  the  predicted  pressure  field  from  driving  a  single  band  was  also  very  similar  for 
both  of  these  models. 

800  Hz  with  Oil  Outside  the  CWC 

In  order  to  study  evanescent  fields  within  the  liquid-backed  CWC,  it  was  again  operated  with 
castor  oil  inside  the  LEXAN  box.  The  transfer  matrix  was  then  measured  for  the  oil-backed  CWC  at 
800  Hz  at  the  70  field  points  used  previously  for  the  air-backed  case.  The  measured  transfer  matrix 
was  used  to  calculate  a  set  of  drives  to  produce  a  55-m/s  evanescent  wave  in  a  window  running  from 
-8.5  to  +8.5  cm.  Figure  39  shows  the  predicted  and  measured  pressure  fields  within  the  CWC  for  this 
set  of  drives.  These  results  are  very  similar  to  the  air-backed  case  shown  in  Fig.  36.  This  is  not 
surprising  since  the  previous  analysis  has  shown  that  the  characteristics  of  the  CWC  do  not  change 
very  much  at  800  Hz  whether  the  CWC  is  air  or  oil  backed  and  the  phase  speeds  and  window  widths 
specified  for  these  two  cases  are  almost  the  same.  Any  differences  in  characteristics  of  the  CWC  are 
automatically  taken  into  account  by  the  measured  transfer  matrix. 

Drives  were  then  computed  specifying  the  same  desired  pressure  Field  but  using  the  CHIEF-I 
model  to  determine  the  transfer  matrix  rather  than  measured  values.  Figure  40  shows  the  measured 
pressure  field  within  the  CWC  for  this  set  of  drives.  The  measured  pressure  field  is  similar  to  Fig.  37 
but  shows  slightly  less  uniformity  inside  the  window. 

Drives  were  then  computed  specifying  the  same  desired  pressure  field  but  using  the  ATILA- 
CHIEF  model  to  determine  the  transfer  matrix,  and  predicted  pressure  fields  were  computed.  Figure 
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Fig.  36  —  Pressures  within  the  air-backed  CWC  for  a  57-m/s  evanescent  field  at  800 
Hz  in  a  window  running  from  -8.2  to  +8.2  cm  using  the  measured  transfer  matrix. 


41  shows  the  predicted  and  measured  pressure  Fields  within  the  CWC  for  this  set  of  drives.  These 
results  are  similar  to  those  obtained  with  the  CHIEF-I  model.  Neither  model  performed  as  well  as 
when  using  the  measured  transfer  matrix. 

The  measured  transfer  matrix  was  then  used  to  calculate  a  set  of  drives  to  produce  a  92-m/s 
evanescent  wave  in  the  same  window  (from  -8.5  to  +8.5  cm).  Figure  42  shows  the  predicted  and 
measured  pressure  fields  within  the  CWC  for  this  set  of  drives.  The  uniformity  inside  the  window  is 
not  quite  as  good  as  it  was  for  the  55-m/s  case.  It  is  also  apparent  that  the  pressure  field  was  not 
driven  down  as  low  outside  the  window. 

Drives  were  then  computed  specifying  the  same  desired  pressure  field  but  using  the  CHIEF-I 
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Fig.  37  —  Pressures  within  the  air-backed  CWC  for  a  57-m/s  evanescent  pressure  field 
at  800  Hz  in  a  window  running  from  -8.2  to  +8.2  cm  using  CHIEF-I. 


model  to  determine  the  transfer  matrix  rather  than  measured  values.  These  results  are  shown  in  Fig. 
43.  The  field  shows  much  more  variation  inside  the  window  than  the  results  for  the  measured  transfer 
matrix. 

Finally,  drives  were  computed  specifying  the  same  desired  pressure  field  but  using  the  ATILA- 
CHIEF  model  to  determine  the  transfer  matrix.  These  results,  shown  in  Fig.  44,  indicate  a  significant 
improvement  of  the  uniformity  of  the  field  inside  the  window  compared  with  the  CHIEF-I  model. 

In  the  evanescent-field  results  presented  thus  far,  the  window  sizes  and  evanescent  phase  speeds 
were  chosen  to  make  Eq.  (1)  equal  zero  at  the  effective  acoustic  wavenumber  within  the  CWC.  This 
minimizes  contamination  of  the  evanescent  field  from  acoustic  energy  propagating  axially  within  the 
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Fig.  38  —  Pressures  within  the  air-backed  CWC  for  a  57-m/s  evanescent  pressure  field 
at  800  Hz  in  a  window  running  from  -8.2  to  +8.2  cm  using  RING-CHIEF. 


CWC.  The  window  half-width  w  was  determined  from 


Wn 


nCaCe 


2 f(ca~  O  ’ 


n  =  1,  2,  3,  ... 


(68) 


where  ca  is  the  propagation  speed  for  sound  of  this  frequency  freely  propagating  down  the  inside  of 
the  CWC,  ce  is  the  desired  evanescent  phase  speed,  and  /  is  the  frequency.  Figure  45  shows  a  plot  of 
Eq.  (1)  for  the  92-m/s  case  described  in  Figs.  42  through  44.  A  value  of  one  was  chosen  for  n  to  find 
the  narrowest  window  that  would  satisfy  the  equation.  The  evanescent  wavenumber  ke  =  to /ce  and  the 
effective  acoustic  wavenumber  within  the  CWC  ka  =  i o/ca  are  also  indicated  in  Fig.  45.  Other 
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Fig.  39  —  Pressures  within  the  oil-backed  CWC  for  a  55-m/s  evanescent  field  at  800 
Hz  in  a  window  running  from  -8.5  to  +8.5  cm  using  the  measured  transfer  matrix. 


evanescent  phase  speeds  that  satisfy  Eq.  (1)  for  this  window  width  can  be  found  from 

1 

Cgm  ~ 1  ~  m~  ’ 

~a  +  2/W 

where  larger  values  of  m  give  lower  allowed  evanescent  phase  speeds. 

In  an  attempt  to  improve  the  performance  of  the  CWC,  it  was  decided  to  try  and  eliminate  any 
monopole  component  of  the  acoustic  energy  by  forcing  Eq.  (1)  to  also  equal  zero  for  a  wavenumber  of 


m  =  1,  2,  3,  - 


(69) 
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Fig.  40  —  Pressure  within  the  oil-backed  CWC  for  a  55-m/s  evanescent  field  at  800 
Hz  in  a  window  running  from  -8.5  to  +8.5  cm  using  CHIEF-I. 


zero.  This  can  be  achieved  with  a  window  half-width  computed  from  Eq.  (68)  using  an  evanescent 
phase  speed  equal  to  cj 2.  This  gives  window  widths  that  contain  integral  numbers  of  evanescent 
wavelengths.  Figure  46  shows  the  wavenumber  spectrum  for  the  oil-backed  CWC  ( ca  =  282  m/s) 
satisfying  this  requirement  at  800  Hz  with  m  -  n  =  1.  This  gives  an  evanescent  phase  speed  of  141 
m/s.  Unfortunately,  it  also  gives  a  window  half-width  of  17.6  cm,  which  requires  a  window 
considerably  larger  than  the  prototype  CWC.  This  would  not  be  a  problem  with  a  full  scale  CWC, 
however,  because  it  would  be  several  times  as  long.  Since  the  window  half-width  w  is  inversely 
proportional  to  the  frequency  [see  Eq.  (68)],  it  was  decided  to  reduce  the  window  width  to  a  usable 
size  for  the  prototype  by  increasing  the  operating  frequency. 
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Fig.  41  —  Pressures  within  the  oil-backed  CWC  for  a  55-m/s  evanescent  field  at  800 
Hz  in  a  window  running  from  -8.5  to  +8.5  cm  using  ATILA-CHIEF. 


1600  Hz  with  Air  Outside  the  CWC 

The  frequency  of  operation  of  the  CWC  was  raised  to  1600  Hz  to  permit  the  use  of  a  window- 
width,  evanescent-phase-speed  combination  that  made  Eq.  (1)  equal  zero  when  the  wavenumber 
equalled  both  zero  and  the  ka  within  the  CWC.  Figure  47  shows  a  wavenumber  spectrum  for  the  air- 
backed  CWC  operating  at  1600  Hz.  For  this  case  ca  =  284  m/s,  ce  =  71  m/s,  and  w  =  8.9  cm.  It  can 
be  seen  that  the  criteria  mentioned  above  are  satisfied.  This  window  is  four  evanescent  wavelengths 
wide. 


The  transfer  matrix  was  measured  for  the  air-backed  CWC  at  1600  Hz  at  the  70  field  points  used 
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Fig.  42  —  Pressures  within  the  oil-backed  CWC  for  a  92-m/s  evanescent  field  at  800 
Hz  in  a  window  running  from  -8.5  to  +8.5  cm  using  the  measured  transfer  matrix. 


previously  at  800  Hz.  The  measured  transfer  matrix  was  then  used  to  calculate  a  set  of  drives  to 
produce  a  71-m/s  evanescent  wave  in  a  window  running  from  -8.9  to  +8.9  cm.  Figure  48  shows  the 
predicted  and  measured  pressure  fields  within  the  CWC  for  this  set  of  drives.  The  measured  pressure 
shows  good  agreement  with  the  predicted  pressure. 

The  CWC  was  then  modeled  using  CHIEF-I  with  band  4  driven  alone  at  1600  Hz.  The  predicted 
and  measured  pressures  within  the  CWC  for  this  case  are  shown  in  Fig  49.  These  results  show  that 
the  agreement  between  the  modeled  and  measured  pressures  was  not  as  good  as  it  was  at  800  Hz.  The 
higher  frequency  (1600  Hz)  is  just  above  the  third  water-column  resonance  within  the  CWC.  This 
places  more  demands  on  how  accurately  an  idealized  model  duplicates  the  actual  CWC.  One  example 
of  this  is  the  influence  of  the  air-filled  LEXAN  box.  If  the  box  changes  the  conditions  at  the  ends  of 
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Fig.  43  —  Pressure  within  the  oil-backed  CWC  for  a  92-m/s  evanescent  field  at  800 
Hz  in  a  window  running  from  -8.5  to  +8.5  cm  using  CHIEF-I. 


the  CWC  so  that  the  effective  length  of  the  water  column  within  the  CWC  is  altered,  this  will  change 
the  pressure  field  in  the  water  column.  This  effect  is  more  pronounced  at  higher  frequencies  and 
especially  near  resonances. 

Drives  were  then  computed  specifying  the  same  desired  pressure  field  as  for  Fig.  48  but  using  the 
CHIEF-I  model  to  determine  the  transfer  matrix.  Figure  50  shows  the  measured  pressure  field  within 
the  CWC  for  this  set  of  drives.  These  results  show  poor  performance  in  both  amplitude  and  phase. 
This  is  most  likely  due  to  the  inaccuracy  of  the  model.  Although  the  differences  between  the  model 
and  the  measured  pressures  may  not  appear  to  be  very  great  when  a  single  band  is  driven,  the  results 
from  single-band  drives  are  being  used  here  to  calculate  a  set  of  40  complex  drives  intended  to 
produce  an  evanescent  field  within  the  CWC  that  has  a  wavelength  of  only  4.4  cm,  thereby  increasing 
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Fig.  44  —  Pressures  within  the  oil-backed  CWC  for  a  92-m/s  evanescent  field  at  800 
Hz  in  a  window  running  from  -8.5  to  +8.5  cm  using  AULA-CHIEF. 


the  effects  of  the  errors.  This  situation  is  aggravated  by  the  fact  that  at  this  evanescent  wave  speed 
there  are  only  7.4  bands  per  wavelength.  This  is  below  the  minimum  of  10  bands  per  wavelength 
recommended  for  planar  arrays  [19]. 

Drives  were  also  computed  for  this  case  using  the  RING-CHIEF  model  to  determine  the  transfer 
matrix.  Figure  51  shows  the  predicted  and  measured  pressure  fields  within  the  CWC  for  this  set  of 
drives.  These  data  also  show  poor  performance  as  with  the  CHIEF-I  model. 

To  determine  if  the  performance  of  the  model  would  improve  with  higher  evanescent  wave 
speeds,  drives  were  calculated  for  the  same  window  size  but  with  an  evanescent  wave  speed  of  142 
m/s  using  the  measured  transfer  matrix.  Figure  52  shows  the  predicted  and  measured  pressure  fields 
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Fig.  45  —  Wavenumber  spectrum  [Eq.  (1)]  for  a  92-m/s  pressure  field  at  800  Hz  in  a 
window  from  -8.5  to  +8.5  cm  in  the  oil-backed  CWC. 


within  the  CWC  for  this  set  of  drives.  The  results  show  uniformity  inside  the  window  that  is  almost 
as  good  as  the  71 -m/s  case  for  the  measured  transfer  matrix. 

Drives  were  then  computed  for  this  case  using  the  CHIEF-I  model  and  the  RING-CHIEF  model 
to  determine  the  transfer  matrices.  Figures  53  and  54  show  the  pressure  fields  within  the  CWC  for 
these  two  cases.  These  results  show  general  agreement  with  the  results  from  the  measured  transfer 
matrix  although  the  amplitude  uniformity  inside  the  window  is  not  as  good.  The  greatly  improved 
results  with  the  models  (compared  to  the  71-m/s  case)  do  indicate  that  the  demands  on  the  accuracy  of 
the  models  are  lessened  for  higher  evanescent  wave  speeds. 

1600  Hz  with  Oil  Outside  the  CWC 

To  evaluate  the  oil-backed  configuration  at  1600  Hz,  the  CWC  was  operated  with  castor  oil 
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Fig.  46  —  Wavenumber  spectrum  for  a  141-m/s  evanescent  pressure  field  at  800  Hz  in 
a  window  from  -17.6  to  +17.6  cm  in  the  oil-backed  CWC  (ca  =  282  m/s). 


inside  the  LEXAN  box.  The  transfer  matrix  was  measured  with  this  configuration  and  used  to 
calculate  drives  to  produce  a  66-m/s  evanescent  wave  in  a  window  running  from  -8.3  to  +8.3  cm. 
Figure  55  shows  the  predicted  and  measured  pressure  fields  within  the  CWC  for  this  set  of  drives. 

This  evanescent  phase  speed  corresponds  to  7  bands  per  wavelength.  This  data  shows  that  the 
uniformity  of  the  pressure  inside  the  window  is  degrading  even  when  using  the  measured  transfer 
matrix. 

Radiation-impedance  matrices  for  1600  Hz  were  computed  with  CHIEF  for  use  by  RING-CHIEF 
and  ATILA-CHIEF.  Predicted  pressures  were  then  computed  using  both  of  these  models  with  band  4 
driven  alone.  Figure  56  shows  the  pressures  predicted  by  RING-CHIEF  and  ATILA-CHIEF  and  the 
measured  pressures  within  the  CWC  for  this  case.  This  data  shows  that  ATILA-CHIEF  definitely  does 
a  better  job  of  modeling  the  CWC  at  1600  Hz  than  RING-CHIEF  but  still  shows  significant  variation 
from  the  measured  pressure. 
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Fig.  47  —  Wavenumber  spectrum  for  a  71-m/s  evanescent  pressure  field  at  1600  Hz  in 
a  window  from  -8.9  to  +8.9  cm  in  the  air-backed  CWC  (ca  =  284  m/s). 


Drives  were  then  computed  to  produce  a  66  m/s  evanescent  wave  using  the  ATTLA-CHIEF 
model  to  determine  the  transfer  matrix.  Figure  57  shows  the  predicted  and  measured  pressure  fields 
within  the  CWC  for  this  set  of  drives.  This  shows  very  poor  uniformity  as  with  the  air-backed  case  at 
71  m/s  for  the  same  reasons. 

To  reduce  the  number  of  bands  per  evanescent  wavelength,  drives  were  computed  for  an  88-m/s 
evanescent  wave  in  the  same  window  using  the  measured  transfer  matrix.  This  corresponds  to  9.2 
bands  per  evanescent  wavelength.  Figure  58  shows  the  predicted  and  measured  pressure  field  within 
the  CWC  for  this  set  of  drives.  This  shows  much  better  uniformity  than  for  the  66-m/s  field. 

Drives  were  computed  for  this  case  using  the  ATILA-CHIEF  model  to  determine  the  transfer 
matrix.  Figure  59  shows  the  pressure  fields  within  the  CWC  for  this  set  of  drives.  The  results  show 
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Fig.  48  —  Pressures  within  the  air-backed  CWC  for  a  71-m/s  evanescent  field  at  1600 
Hz  in  a  window  running  from  -8.9  to  +8.9  cm  using  the  measured  transfer  matrix. 


much  better  performance  than  for  the  ATILA-CHIEF  model  at  66-m/s  but  still  not  as  good  as  with  the 
drives  computed  using  the  measured  transfer  matrix. 

Radial  Variation  in  Pressure 

On  numerous  occasions,  at  various  evanescent  wave  speeds,  attempts  were  made  to  measure  any 
variation  in  the  pressure  within  the  CWC  as  a  function  of  radial  position.  This  was  done  by  deflecting 
the  monofilament  line,  and  therefore  the  hydrophone,  from  the  centerline  to  the  inner  wall  of  the  CWC 
at  a  fixed  z  position  while  measuring  the  amplitude.  No  significant  variations  in  amplitude  were  ever 
observed.  This  uniformity  was  made  possibly  because  of  the  relatively  slow  natural  propagation  speed 
within  the  PVDF  CWC.  Although  the  lowest  evanescent  speed  for  which  data  is  presented  (55-m/s)  is 
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Fig.  49  —  Predicted  and  measured  pressures  within  the  air-backed  CWC  with  only 
band  4  driven  at  1600  Hz.  The  predicted  pressures  were  calculated  using  the  CHIEF-I 

model. 


1/27  of  the  sound  speed  in  open  water,  this  only  requires  a  wavenumber  multiplier  of  about  5  for  the 
PVDF  CWC.  As  can  be  seen  from  Fig.  2,  the  pressure  would  only  be  expected  to  drop  by  about  0.2 
dB  at  the  center  of  the  CWC  relative  to  the  wall  for  this  wavenumber  multiplier.  If  a  rigid-walled 
CWC  were  used,  the  pressure  drop  near  the  center  would  be  about  6  dB. 

CONCLUSIONS  AND  RECOMMENDATIONS 

An  investigation  of  a  prototype  CWC  has  been  performed.  The  study  included  both  theoretical 
and  experimental  investigations  of  the  CWC.  Various  methods  for  determining  the  transfer  matrix  for 
the  CWC  were  explored.  The  results  of  the  investigations  show  that  reasonably  uniform  nonacoustic- 
wavenumber  pressure  fields  can  be  created  within  the  CWC. 
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Fig.  50  —  Pressures  within  the  air-backed  CWC  for  a  71-m/s  evanescent  pressure  field 
at  1600  Hz  in  a  window  running  from  -8.9  to  +8.9  cm  using  CHIEF-1. 


The  analysis  of  the  pressure  field  inside  an  infinite-length  cylindrical  shell  with  a  specified 
normal  velocity  showed  the  radial  dependence  of  the  evanescent  pressure  field  inside  the  shell.  From 
this  it  was  concluded  that  the  wavenumber  multiplier  should  be  kept  small  to  keep  the  pressure  within 
the  CWC  from  dropping  rapidly  as  the  distance  from  the  wall  is  increased.  This  can  be  accomplished 
if  a  CWC  is  constructed  with  a  low  natural  propagation  speed  in  the  water  within  the  CWC.  The 
analysis  of  an  infinite-length,  cylindrical,  elastic  tube  filled  with  a  fluid  and  surrounded  by  another 
fluid  showed  that  a  PVDF  CWC  can  be  constructed  that  has  propagation  speeds  in  the  water  within 
the  CWC  that  are  less  than  one-fifth  the  open-water  sound  speed.  These  results  were  verified  by 
calculating  the  propagation  speed  within  the  CWC  using  measured  resonant  frequencies. 

The  measured  pressure  fields  within  the  CWC  with  a  single  band  driven  were  compared  with  the 
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Fig.  51  —  Pressures  within  the  air-backed  CWC  for  a  71 -m/s  evanescent  field  at  1600 
Hz  in  a  window  running  from  -8.9  to  +8.9  cm  using  RING-CHIEF. 


pressure  fields  predicted  by  the  CHIEF-I,  RING-CHIEF,  and  ATILA-CHIEF  models.  With  a  drive 
frequency  of  800  Hz  the  pressures  predicted  by  the  models  are  very  similar  and  show  fairly  good 
agreement  with  the  measured  pressure.  With  a  drive  frequency  of  1600  Hz  the  ATILA-CHIEF  model 
does  a  better  job  of  matching  the  measured  pressures  than  RING-CHIEF  but  still  shows  significant 
differences. 

Evanescent  pressure  fields  were  created  within  the  CWC  with  drives  that  were  computed  using 
transfer  matrices  that  were  either  measured  or  determined  using  the  models.  The  technique  of  adding 
a  small  real  quantity  to  the  diagonals  of  the  reduced  transfer  matrix  resulted  in  drives  with  more 
uniform  amplitudes  and  phases  that  were  more  usable  by  the  drive  system  to  produce  uniform 
evanescent  fields.  It  was  also  found  that  the  uniformity  of  the  evanescent  field  was  poor  using 
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Fig.  52  —  Pressures  within  the  air-backed  CWC  for  a  142-m/s  evanescent  field  at 
1600  Hz  in  a  window  running  from  -8.9  to  +8.9  cm  using  the  measured  transfer 

matrix. 


modeled  transfer  matrices  if  the  pressure  within  the  CWC  was  allowed  to  remain  large  out  to  the  ends 
of  the  CWC.  The  performance  was  greatly  improved  if  drives  were  computed  that  forced  the  pressure 
down  in  the  region  between  the  window  and  the  ends  of  the  CWC. 

With  800-Hz  drive,  the  uniformity  of  the  evanescent  pressure  field  inside  the  window  was  about 
as  good  with  any  of  the  models.  However,  none  of  the  models  produced  evanescent-field  uniformity 
that  was  as  good  as  that  achieved  from  measured  transfer  matrices. 

The  technique  of  minimizing  any  monopole  energy  by  forcing  Eq.  (1)  to  equal  zero  when  the 
wavenumber  equals  zero  seemed  to  improve  the  performance  of  the  CWC  slightly  when  using 
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Fig.  53  —  Pressures  within  the  air-backed  CWC  for  a  142-m/s  evanescent  field  at 
1600  Hz  in  a  window  running  from  -8.9  to  +8.9  cm  using  CHIEF-I. 


measured  transfer  matrices.  The  agreement  with  predicted  pressures  throughout  most  of  the  window 
was  generally  within  ±1.5  dB.  Unfortunately,  the  shift  to  higher  frequency  (required  because  of  the 
short  length  of  the  prototype)  also  brought  about  more  differences  between  measured  pressures  within 
the  CWC  and  pressures  predicted  by  the  models.  Using  transfer  matrices  computed  from  models  at 
1600  Hz,  the  uniformity  of  the  evanescent  pressure  field  inside  the  window  was  definitely  better  with 
the  more  accurate  ATILA-CHIEF  model.  This  difference  was  greatest  for  higher  evanescent 
wavenumbers  (i.e.,  lower  evanescent  wave  speeds).  The  performance  at  1600  Hz  with  drives 
computed  using  measured  transfer  matrices  was  considerable  better  than  the  results  using  transfer 
matrices  computed  from  models. 


It  was  found  that  the  criterion  of  a  minimum  of  ten  bands  per  evanescent  wavelength  that  was 
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Fig.  54  —  Pressures  within  the  air-backed  CWC  for  a  142-m/s  evanescent  field  at 
1600  Hz  in  a  window  running  from  -8.9  to  +8.9  cm  using  RING-CHIEF. 


previously  recommended  for  planar  arrays  is  about  right  for  the  CWC.  Degradation  of  the  evanescent 
field  was  obvious  at  seven  bands  per  wavelength  with  measured  transfer  matrices  and  was  evident 
below  about  nine  bands  per  wavelength  when  using  the  models  to  determine  the  transfer  matrices. 

Overall,  the  performance  of  the  CWC  with  measured  transfer  matrices  was  always  as  good  as  or 
better  than  with  the  models.  In  many  cases  (particularly  at  the  higher  frequency)  the  results  with 
measured  transfer  matrices  were  much  better  than  with  the  models.  This  should  not  be  surprising 
since  the  models  are  idealized  approximations  of  a  perfect  CWC  while  the  measured  transfer  matrices 
determine  the  true  interactions  for  the  actual  CWC  in  its  operating  environment.  This  situation  would 
be  even  worse  with  an  invariably  more  complex  full-size  CWC  than  with  the  simple  prototype.  This 
leads  to  the  conclusion  that  a  practical  CWC  should  be  capable  of  measuring  the  transfer  matrix  rather 
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Fig.  55  —  Pressures  within  the  oil-backed  CWC  for  a  66-m/s  evanescent  field  at  1600 
Hz  in  a  window  running  from  -8.3  to  +8.3  cm  using  the  measured  transfer  matrix. 


than  relying  on  a  model  of  the  system. 

An  operational  CWC  would  have  a  diameter  several  times  as  large  as  the  prototype  and  be 
several  times  as  long  in  order  to  accommodate  towed-sonar-array  modules.  It  must  have  narrower 
bands  in  order  to  produce  higher  evanescent  wavenumbers.  This  fact,  combined  with  the  additional 
length,  means  that  there  would  have  to  be  many  times  the  number  of  driven  bands  used  in  the 
prototype.  If  each  drive  must  be  independently  controlled,  there  must  be  hundreds  (perhaps  thousands) 
of  sets  of  drive  electronics  (i.e.,  power  amplifiers,  etc.).  It  is  possible  that  this  number  could  be 
dramatically  reduced,  however.  It  was  noted  that  the  drives  to  the  interior  bands  (bands  that  were  not 
near  the  ends  of  the  CWC)  tended  to  have  relatively  uniform  amplitudes  and  a  constant  phase 
increment  that  coincided  with  the  phase  of  the  evanescent  wave.  It  is  therefore  possible  that  many  of 
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Fig.  56  —  Pressures  within  the  oil-backed  CWC  with  only  band  4  driven  at  1600  Hz. 
Predicted  pressures  were  calculated  using  RING-CHIEF  and  ATTLA-CHIEF  models. 

the  interior  bands  could  be  driven  in  parallel  by  a  single  set  of  drives  that  exploits  this  periodicity  (one 
period  per  evanescent  wavelength).  This  would  have  to  be  confirmed  by  additional  experimentation, 
however,  with  a  CWC  that  is  longer  and  has  more  bands  than  the  current  prototype. 
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Fig.  57  —  Pressures  within  the  oil-backed  CWC  for  a  66-m/s  evanescent  field  at  1600 
Hz  in  a  window  running  from  -8.3  to  +8.3  cm  using  ATILA-CHIEF. 
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Fig.  58  —  Pressures  within  the  oil-backed  CWC  for  a  88-m/s  evanescent  field  at  1600 
Hz  in  a  window  running  from  -8.3  to  +8.3  cm  using  the  measured  transfer  matrix. 
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Fig.  59  —  Pressures  within  the  oil-backed  CWC  for  a  88-m/s  evanescent  field  at  1600 
Hz  in  a  window  running  from  -8.3  to  +8.3  cm  using  ATILA-CHIEF. 
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PROPAGATION  SPEED  WITHIN  AN  ELASTIC  TUBE 


To  determine  the  predicted  propagation  within  the  PVDF  CWC  prototype,  a  study  of  the  propagation 
speed  in  the  fluid  within  an  infinite-length  elastic  tube  surrounded  by  another  fluid  was  undertaken. 
Figure  A1  shows  a  cross-section  view  of  the  tube  and  the  fluid  regions.  Region  I  is  filled  with  a  fluid 
with  density  p4-  and  sound  speed  c4-.  Region  II  consists  of  an  elastic  tube  with  density  pp  dilatational 


Region  III  (fluid) 

p  ,o 

's'  s 


Fig.  A1  —  Infinite-length,  cylindrical,  fluid-elastic- 
fluid  arrangement. 
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propagation  speed  cd,  and  transverse  propagation  speed  cr  Region  III  is  the  fluid  outside  the  tube 
extending  to  infinity  with  density  p5  and  sound  speed  cs. 

General  Solutions 


Region  I 

The  fluid  in  region  I  can  be  described  by  the  displacement  potential  function  O,  with  [72] 


2  1  3  °l 

V20,  =  _ _ i 

1  2  a#2 


Ci 


dt‘ 


(Al) 


Assuming  radial  symmetry  in  cylindrical  coordinates  (r,8,z),  and  exp(ico()  time  dependence,  Eq.  (Al) 
becomes 


1 

r 


_3_ 

Hr 


30, 

~Hr 


3zO, 

dz2 


*  -  0  , 


(A2) 


where  |  kt  \  =  co/c-.  The  radial  and  axial  displacements  and  pressure  in  region  I  are  related  to  the  potential 
function  by 


30, 


U‘r  dr 


(A3) 


30, 


““  dz 


(A4) 


320 


Pi  =  “Pi 


d2t 


(A5) 


Separating  O,  into  the  product  of  a  function  of  r  only  R\(r)  and  a  function  of  z  only  Z,(z),  one  has 

3 2/?  i  1  dR  |  o 

_ I  +  1  _ L  +  k;R,  =  0 


(A6) 


dr 


2  r  3r 


and 


3  Z,  ,2 


3z' 


I  +  *z2Z,  =  0  , 


(AT) 


where  kir  and  kz  are  the  r  and  z  components  of  kr  Equation  (A6)  has  the  solution 
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Hr)  =  AlJ0 (kirr)  +  BlY0 M  »  (A8) 

where  Aj  and  B,  are  arbitrary  constants.  Requiring  that  the  solution  be  finite  at  r  =  0  forces  Bj  to  equal 
zero.  Therefore, 


(A9) 

Equation  (A7)  has  the  solution 

Z,(z)  =  Cj  exp (~ikzz)  .  (A10) 

The  solution  for  <f>j  can  then  be  written  as 

=  A  J0(kirr)  exp[i(o)f  ~kzz)\  ,  (All) 

where  A  is  an  arbitrary  constant  combining  Aj  and  Cv 
Region  II 

The  elastic  solid  in  region  II  can  be  described  [73]  by  the  scalar  potential  and  the  vector  potential 
W.  Assuming  axial  symmetry,  the  radial  and  axial  displacements  in  region  II  are  related  to  these  potential 
functions  by 


dS’lI  _  dW 
dr  dz 


(A12) 


and 


ao>n  j 
u  = -  +  — 

z  dz  r 


d(rW) 


dr 


where  W(r^)  is  the  azimuthal  component  of  *1'. 

The  potential  functions  On  and  W  must  be  solutions  to 


VzOn  =  — 


1 


2  at2 
cd  dt 


(A13) 


(A14) 


and 


v2  w  -  _i_  w  =  —  -  — 
r2  c,2  dt 2 


(A15) 


respectively,  where  cd  is  the  dilatational  wave  speed  and  ct  is  the  shear  wave  speed  in  the  elastic  solid. 
Assuming  solutions  of  Eqs.  (A14)  and  (A15)  in  the  form  of  waves  traveling  in  the  z  direction,  one  has 
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=  [CJ0(kdrr)  *  DY0(kdrr]  cxp[i(cot  -kzz)}  (A16) 

and 

W  =  [FJ^k'S)  +  GYx{klrr ]  exp[i(co/  -kzz)}  ,  (A17) 


where 


and 


(A18) 


(A19) 


Region  III 

The  fluid  in  region  III  can  be  described  by  the  potential  function  <J>m,  which  is  identical  with  the 
function  for  region  I.  Following  the  same  procedure  as  in  region  I  one  arrives  at  the  solution 

^m(r)  =  Ani^o(^s/)  +  ^lll^o(^rrr)  »  (A20) 


where  k ^  is  the  r  components  of  ks.  Assuming  an  outgoing  wave  and  combining  constants  gives 


Rm(r)  =  BH^\ksrr)  , 


(A21) 


where  Hq®  is  a  Hankel  function  of  the  second  kind.  The  solution  for  region  III  can  then  be  written  as 


=  BHc !2>(V)  exp[i(o)t  -  kzz)\  . 


(A22) 


Boundary  Conditions 

The  boundary  conditions  relating  the  pressures,  stresses,  and  displacements  at  the  region  boundaries 
are: 

1.  The  negative  of  the  pressure  in  region  I  must  equal  the  normal  stress  in  region  II  at  r  =  a. 

2.  The  normal  component  of  the  particle  displacement  in  region  I  must  equal  the  normal  component 

of  the  surface  displacement  in  region  II  at  r  =  a. 

3.  The  tangential  stress  in  region  II  must  be  zero  at  r  =  a. 

4.  The  tangential  stress  in  region  II  must  be  zero  at  r  =  b. 

5.  The  negative  of  the  pressure  in  region  III  must  equal  the  normal  stress  in  region  II  at  r  =  b. 

6.  The  normal  component  of  the  particle  displacement  in  region  III  must  equal  the  normal  component 
of  the  surface  displacement  in  region  II  at  r  =  b . 

The  following  relationships  for  this  axially-symmetric  situation  are  useful  for  expressing  these 
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boundary  conditions  mathematically. 
The  non-zero  strains  are  [74]: 


dUr 

dr 


eBB  ~ 


ezz  = 


dUz 

dz 


e  -  e  = _ 

rr  2r  2 


dur  duz 
dz  dr 


A  =  e  +  eaa  +  e  , 

rr  09  zz  » 

where  A  is  the  dilatation. 

The  stress-strain  relationships  are  [75]: 

xrr  =  XA  +  2 fierr  , 

xee  =  ^  +  2fieBB  > 
xzz  =  XA  +  2/rezz  , 

Xr2  =  2Perz  » 

where  X  and  /i  are  Lame’s  constants  for  the  elastic  material  in  region  II  related  to  cd  and  c t  by 


2  _  X  +  2jM  2  -  P 

- >  ci  — 


= 


P/ 


The  boundary  conditions  can  then  be  expressed  mathematically  as: 
B.  C.  #  1 


Trr|  +  Pi\  =  0  * 

rr\ r-a  *1  r-a 


B.  C.  #  2 


Uir\  ~  Ur\  "  0  > 
"lr=a  ^ir=a 


b.  c.  #  3 


(A23) 

(A24) 

(A25) 

(A26) 

(A27) 

(A28) 

(A29) 

(A30) 

(A31) 

(A32) 
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(A33) 


B.  C.  #  4 


(A34) 


B.  C.  #  5 


(A35) 


B.  C.  #  6 


Propagation  Speeds 

By  inserting  the  appropriate  items  into  the  boundary  condition  equations,  taking  the  required 
derivatives,  and  considerable  algebraic  manipulation  (carried  out  with  the  aid  of  SMP  [76]),  one  can  arrive 
at  the  6x6  coefficient  matrix  shown  in  Fig.  A2.  After  expressing  all  the  wavenumbers  in  terms  of  kz,  the 
values  of  kz  that  make  the  determinant  of  the  coefficient  matrix  equal  to  zero  give  the  propagation  speeds 
for  the  modes  that  are  supported  at  that  frequency.  A  FORTRAN  computer  program  was  written  to 
evaluate  the  determinant  and  search  for  the  zeros  at  each  specified  frequency.  The  program  automatically 
switches  from  Bessel  functions  to  the  appropriate  modified  Bessel  functions,  whenever  the  function 
arguments  become  imaginary,  using  the  relationships  [77,78]: 


Hf\x)  =  JJx)  -iY„(x)  , 


(A37) 


Jo(ix)  =  W  > 


(A38) 


J^ix)  =  il^x)  , 


(A39) 


Y0(ix)  =  -Ik 0(x)  *  il0(x) 

K 


(A40) 
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Fig.  A 2  —  Coefficient  matrix  for  the  fluid-elastic-fluid  problem. 


APPENDIX  B 


CWCMDL  PROGRAM 


Description  of  the  Program 

A  FORTRAN  computer  program  called  CWCMDL  was  written  to  calculate  sets  of  electrical  band 
drives  that  could  be  used  to  generate  evanescent  pressure  fields  within  the  CWC.  The  program  computes 
the  complex  band  drives  using  transfer  matrices  determined  from  either  CHIEF-I,  RING-CHIEF,  or 
ATILA-CHIEF  models  of  the  CWC.  The  program  can  also  compute  the  predicted  pressure  field  within 
the  CWC  from  a  specified  set  of  band  drives. 

Data  for  the  CWC  Modeling  Program  (CWCMDL) 


CA0807A 

!(CHR)  Name  of  the  run  (run_name) 

3 

!(I)  Function  (l=calc  drives,  2=map,  3=both)  (oper_mode) 

3 

!(I)  Model  type  (l=rigid,2=imped.,3=vel.  file)  (mod_type) 

1 

!(I)  Printout  level  (prntjvl) 

55,  92 

!(R)  Phase  speeds  of  evanescent  waves  (m/s)  (phasjspd) 

1  2 

!(I)  Numbers  of  phase  speeds  to  map  (mapjspd_num) 

800 

!(R)  Frequency  in  (Hz)  (freq) 

1000 

!(R)  Density  of  fluid  (kg/m3)  (rho) 

1485 

!(R)  Speed  of  sound  in  fluid  (m/s)  (snd_spd) 

.0155 

!(R)  CWC  inner  radius  (m)  (r_inner) 

.01625 

!(R)  CWC  outer  radius  (m)  (r_outer) 

.006 

!(R)  Width  of  grounded  band  (m)  (gm_bnd_wid) 

.006 

!(R)  Width  of  active  band  (m)  (act_bnd_wid) 

40 

!(I)  Number  of  active  bands  (num_bnds) 

alibands 

!(CHR)  File  containing  the  driven  band  numbers  (bnd^file) 

70 

!(I)  Number  of  points  for  the  transfer  matrix  (num_pts) 

-.138 

!(R)  Minimum  z  for  transfer  matrix  points  (m)  (min_zp) 

.138 

!(R)  Maximum  z  for  transfer  matrix  points  (m)  (max_zp) 

.1 

!(R)  Percent  average  added  to  matrix  diagonals  (alpha) 

-.085 

!(R)  Minimum  z  for  window  (m)  (min_zc) 

.085 

!(R)  Maximum  z  for  window  (m)  (max_zc) 

0 

!(R)  x  coordinate  for  map  (m)  (x_map) 

0 

!(R)  y  coordinate  for  map  (m)  (y_map) 

-.135 

!(R)  Minimum  z  to  map  (m)  (z_map_min) 

.005 

!(R)  Delta  z  for  mapping  (m)  (z_map_delta) 

.135 

!(R)  Maximum  z  to  map  (m)  (z_map_max) 

CA0800.VEL 

!(CHR)  File  containing  the  single-drive  velocities  (vel  file) 

Fig.  B1  —  Sample  parameter- input  file  for  CWCMDL. 


The  parameters  that  CWCMDL  uses  in  each  run  are  contained  in  a  text  file.  A  sample  of  the  input- 
parameter  file  is  given  in  Fig.  Bl.  The  program  reads  the  parameters  from  the  input  file  and  then 
performs  the  necessary  computations. 
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Flow  of  Processing 

The  flow  of  processing  in  CWCMDL  is: 

1.  Call  the  subroutine  INPUT,  which  opens  the  input-parameter  file  and  reads  the  parameters  for  the 
run. 

2.  Call  the  subroutine  LOADBNDS,  which  reads  the  driven  band  numbers  from  the  specified  file.  If 
the  file  name  "allbands"  is  specified,  no  driven-band-number  file  is  accessed  and  all  the  bands  are 
used. 

3.  If  calculation  of  drives  is  requested,  call  the  subroutine  CALCDRV,  which: 

a.  Calls  either  CHFSUB  or  CHZSUB  to  set  up  CHIEF  and  generate  the  CHIEF  matrices. 
CHFSUB  sets  up  CHIEF  for  the  case  where  all  the  normal  surface  velocities  are  specified  and 
CHZSUB  sets  up  CHIEF  for  the  case  where  interior-wall  impedances  are  specified. 

b.  Calculates  the  columns  of  the  M  matrix  by  utilizing  CHIEF  (with  or  without  ATILA)  to 
determine  the  pressure  at  each  field  point  while  each  band  is  driven  alone.  If  the  RING-CHIEF 
or  ATILA-CHIEF  models  are  being  used  the  resultant  normal  surface  velocities  for  each  surface 
when  a  single  band  is  excited  are  read  from  the  specified  file.  CALCDRV  checks  for  and 
utilizes  any  symmetry  of  field  points  to  reduce  the  number  of  unique  columns  that  have  to  be 
calculated. 

c.  Fills  any  remaining  non-unique  columns  with  the  appropriate  results  from  the  unique  columns. 

d.  Calculates  the  desired  pressures  P .  Multiple  P  vectors  can  be  specified  (e.g.,  different  phase 
speeds  and  window  widths)  since  solving  for  the  multiple  drives  with  back  substitution  requires 
very  little  additional  computation. 

e.  Call  the  subroutine  DCXLSTSQ  (Double-precision  Complex,  LeaST  SQuares)  to  collapse  the 
m  equations  to  n  equations.  The  solution  of  the  new  system  is  the  least-square  solution  of  the 
original  system. 

f.  Calculates  a  as  a  percentage  of  the  average  value  of  the  absolute  amplitudes  of  the  elements 
of  the  reduced  matrix  Af  and  adds  it  to  the  diagonal  elements. 

g.  Perform  lower-upper  decomposition  on  M1. 

h.  Obtains  each  set  of  drives  by  inserting  the  appropriate  desired  pressures  and  using  back- 
substitution. 

i.  Writes  the  resulting  drive  levels  to  files. 

4.  If  maps  of  the  predicted  pressure  fields  are  requested,  calls  the  subroutine  MAPFLDS  to  map  the 
resultant  pressure  fields  as  follows: 

a.  If  necessary  sets  up  CHIEF  by  calling  CHFSUB. 

b.  For  each  set  of  drives  to  be  mapped: 

i.  Calls  LOADDRV,  which  loads  the  drives  from  a  specified  file.  If  the  RING-CHIEF  or 
ATILA-CHIEF  models  are  specified,  the  normal  surface  velocities  are  computed  by 
adding  the  complex  contributions  to  each  surface  as  each  band  is  excited  with  the 
specified  drive. 

ii.  Computes  the  pressure  at  each  position  to  be  mapped  by  calling  the  subroutine  CHFNFP. 

iii.  Writes  the  map  of  the  pressure  field  to  a  file. 

FORTRAN  Code 

Listings  of  the  more  significant  routines  used  by  CWCMDL  follow: 

CWCMDL. FOR  (Main  program) 

program  CWCMDL 

c  This  program  models  the  Cylindrical  Wavenumber  Calibrator  (CWC). 
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c  It  also  calculates  the  band  drives  required  to  produce  an 

c  evanescent  field  inside  the  CWC. 

c  L.  D.  Luker,  NRL/USRD,  24-SEP-1990.  Last  modified  on  9-MAY-1992. 

implicit  none 

include  'CWCMDL. INC /list ' 

character*3  symtyp 
character*4  fldtyp 
integer* 4  i,  j,  nblks 

c  Open  the  main  input  file 

open(unit  =  inp_un,  file  =  inp_file,  status  =  'OLD' ) 

c  Input  the  run  parameters 
call  input(l,  45) 
close  (unit  =  inp_un) 

c  If  desired,  read  the  driven  band  numbers  and  calc,  and  save  drives 
if  ((oper_mode  .  eq.  1)  .or.  (oper_mode  .  eq.  3))  then 
call  loadbnds 
call  calcdrv 
end  if 

c  If  desired,  map  the  pressure  fields 

if  ((oper_mode  .eq,  2)  .or.  (operjmode  .eq.  3))  then 
call  mapflds 
end  if 

c  All  done. 

close  (unit  =  outrun) 

end 


CWCMDLJNC  ( Included  file ) 


c  This  file  (CWCMDL. INC)  is  included  in  all  of  the  Cylindrical 
c  Wavenumber  Calibrator  (CWC)  routines  to  declare  common  variables. 


c  L.  D.  Luker,  NRL/USRD,  24-SEP-1990.  Last  modified  on  12-DEC-1992. 


character 

1  date_time*20,  I Date  and  time 
1  run_name*15  I Name  of  the  run 

character*20 

main  input  data 
main  output  data 
numbers  of  driven  bands 
band  drives 

single  drive  velocities 
pressure  map 

parameter  (inp_file  =  9 CWCMDL . INP ' ) 


1 

inp_file. 

lFile 

containing 

the 

1 

out_file, 

IFile 

containing 

the 

1 

bnd_f ile, 

lFile 

containing 

the 

1 

drv_f ile. 

IFile 

containing 

the 

1 

vel_file, 

lFile 

containing 

the 

1 

map_f ile 

IFile 

containing 

the 

integer*4 

1  inp_un,  iUnit 
1  outrun,  iUnit 
1  bnd_un,  IUnit 
1  drv_un,  iUnit 


number  for  main  input  file 
number  for  main  output  file 
number  for  driven  band  file 
number  for  band  drives 
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c 

c 

c 

c 


1  vel_un, 
1  map_un, 
1  out  len 


iUnit  number  for  single  drive  vel.  file 
l  Unit  number  for  map  output 
I Record  length  for  main  output 


parameter 

parameter 

parameter 

parameter 

parameter 

parameter 

parameter 


( inp_un  =  1 ) 
(bnd_un  =  2) 
(drv_un  =  3) 
(vel_un  =  4) 
(map_un  =  8) 

( out_un  =  9 ) 
(out_len  =  1024) 


integer* 4 
1  max_bnds , 
1  max_pnts , 
1  max_spds 


IMaximum  number  of  bands 

IMaximum  number  of  field  points  or  map  points 
IMaximum  phase  speeds 


parameter  (max_bnds=50,  max_pnts=151 ,  max_spds=9 ) 


integer* 4 
1  oper_mode , 

1  prnt_lvl, 

1  num__bnds , 

1  num_drv_bnds 
1  num__pnts , 

1  num_spds, 

1  num_map_spds 
1  drv_bnd ( max_ 
1  map_spd_num( 
1  mod_type 


lOperating  mode  ( l=calc  drvs, 2=map, 3=both) 

I Printout  level 
INumber  of  bands 
INumber  of  driven  bands 

INumber  of  random  points  for  transfer  matrix 
INumber  of  evanescent  phase  speeds 
INumber  of  evanescent  phase  speeds  to  map 
bnds),  1 Array  containing  numbers  of  driven  bands 

max_spds),  INumbers  of  phase  speeds  to  map 
1 Model  type: 

1  ->  Rigid  tube 

2  ->  Interior  impedance 

3  ->  CHIEF  &  ATILA 


real*4 

1  pi, 

1  rp, 

1  zp, 

1  z_bnd ( max_bnds ) / 
1  max_drv_amp 


IGood  old  pi 

Ir  coordinate  of  matrix  point 
1 z  coordinate  of  matrix  point 
Iz  coordinate  of  band  edges 
IMaximum  normalized  drive  amplitude 


real*4 

1  freq,  IFrequency 

rhO/  IDensity  of  fluid 

snd_spd,  I Speed  of  sound  in  fluid 

phas_spd(max_spds ) /  lPhase  speed  of  evanescent  wave 


runner; 
r_outer f 
gnd_bnd_wid, 
act_bnd_wid/ 
z_arr_min, 
z_arr_max/ 
z_bnds_min, 
z_bnds_max, 
min_zp/ 
max_zp, 
alpha, 
min_zc(max_spds) , 
max_zc (max_spds ) 


I Inner  radius  of  array 
l Outer  radius  of  array 
IWidth  of  grounded  band 
IWidth  of  active  band 
I Minimum  z  coordinate  of  array 

IMaximum  z  coordinate  of  array 

IMinimum  z  coordinate  of  driven  bands 

IMaximum  z  coordinate  of  driven  bands 


IMinimum  point  region 
IMaximum  point  region 


z  coordinate 
z  coordinate 
IValue  added  to  matrix  diagonals  to  force  sol. 
IMinimum  control  region  z  coordinate 
IMaximum  control  region  z  coordinate 


real*4 
1  x_map , 

1  y_map , 

1  z_map_min, 

1  z_map_delta, 
1  z_map_max 


lx  coordinate  to  map 
ly  coordinate  to  map 

IMinimum  z  to  map 
I  Delta  z  for  mapping 
IMaximum  z  to  map 
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complex* 16 

1  wall_imp,  iCWC  inner-wall  impedance 

1  drive ( max_bnds ) 

parameter  (pi=3. 1415927 ) 

common  /cwccom/ 

1  date_time,  run_name, 

1  bnd_file/  drv_file,  vel_file,  map_file, 

1  oper_mode,  prnt_lvl, 

1  num_bnds,  num_drv_bnds ,  num_pnts,  num_spds,  num_map_spds , 
1  map_spd_num,  drv_bnd,  mod_type,  rp,  zp, 

1  z_bnd,  max_drv_amp, 

1  freq,  rho,  snd_spd,  phas_spd, 

1  r_inner,  r_outer, 

1  gnd_bnd_wid,  act_bnd_wid, 

1  z_arr_min,  z_arr_max, 

1  z_bnds_min,  z_bnds_max, 

1  min_zp,  max_zp,  min_zc,  max_zc, 

1  alpha,  x_map,  yjmap, 

1  z_map_min/  z_map_delta,  z_map_raax,  wall_imp,  drive 
c  CHIEF  section 
integer*4 

1  mxars,  IMax  surface  subdivisions  in  sym.  block 

1  mxnfp  IMax  number  of  nearfield  points 

parameter  (mxars  *=  86,  mxnfp  =  max_pnts ) 

integer*4 

1  numfpn,  INumber  of  field  points 

1  infnrm,  1 Field  point  to  normalize  to 

1  mdsize 


parameter  (mdsize  =  mxars  *  max(mxars,  mxnfp)) 


real*4 
1  rmnnrm, 

1  nfpxs(3,  mxnfp),  INearfield  point  coordinates 

1  pnrmnf (mxnfp) 


complex* 8 
1  vel( mxars), 

1  sp(mxars), 

1  nfp(mxnfp), 

1  pmatx (mdsize ) , 
1  vmatx(  mdsize) 


I  Surface  velocities 
I  Surface  pressures 
INearfield  pressures 
1  Pressure  matrix  for  CHIEF 
1  Velocity  matrix  for  CHIEF 


common  /nfinp/  numfpn,  nfpxs 
common  /velsps/  vel,  sp 

common  /nfvals/  nfp,  pnrmnf,  infnrm,  rmnnrm 
common  /chfcl/  pmatx,  vmatx 


INPUT.FOR  (Reads  parameters  from  the  input  file) 

subroutine  input ( inp_col,  out_col) 

c  This  subroutine  reads  the  input  data  from  unit  Minp_un"  starting 
c  at  "inp_col"  and  writes  the  data  to  the  unit  "out_un"  starting 
c  at  "out_col"  preceeded  by  the  data  description. 

c  L.  D.  Luker,  NRL/USRD,  26-SEP-1990.  Last  modified  on  12-DEC-1992. 
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d 


implicit  none 

include  'CWCMDL. INC /list ' 


character  dummy* 7 8 

integer* 4  data_mode,  inp_col,  out_col,  vals 
logical  ask 


data__mode  =  1 
ask  =  .true. 

read(inp_un,  '(a)')  dummy 

call  getchr (data_mode,  ask,  inp_un,  inp_col,  -1,  out_col, 
1  'Name  of  the  run',  run_name) 

open(unit  =  outrun,  file  =  run_name  //  '.OUT', 

1  reel  =  out_len,  status  =  'NEW') 


date_time  =  '  ' 

call  date ( date_time ( 1 : 9 ) ) 

call  time(date_time(12:20) ) 

write ( out_un ,  '(x,''"CWC  Modeling  Program  ('', 

1  a,'')  Run  Name:  #',a/)')  date_time ( 1 : 19 ) ,  run_name 

write(out_un,  ' (x, '  9 " ' 9 ,a,/) ' )  dummy 


map_file  =  run_name 
drv_file  =  run^name 

call  geti ( data_mode ,  ask,  inp_un,  inp_col,  out_un,  out_col,  1, 

1  ' "Operating  mode  ( l=calc,2=map, 3=both) ' ,  oper_mode,  vals) 

call  geti(data_mode,  ask,  inp_un,  inp_col,  out_un,  out_col,  1, 

1  '"Model  type  (l=rigid,  2=imped.,  2=vel.)',  mod_type,  vals) 


call  geti(data_mode,  ask,  inp_un,  inp_col,  out_un,  out_col,  1, 
1  '"Printout  level',  prnt_lvl,  vals) 


call  getr ( data_mode ,  ask,  inp_un,  inp_col,  outrun, 

1  out_col,  max_spds, 

1  '"Phase  speeds  of  evanescent  waves  (m/s)',  phas_spd, 

1  num_spds ) 

call  geti(data_mode,  ask,  inp_un,  inp_col,  out_un, 

1  out_col,  max_spds, 

1  '"Numbers  of  phase  speeds  to  map',  map_spd_num, 

1  num_map_spds ) 


call  getr (data_mode,  ask,  inp_un,  inp_col,  out_un,  out_col,  1, 
1  '"Frequency  (Hz)',  freq,  vals) 


call  getr ( data_mode ,  ask,  inp_un,  inp_col,  out_un,  out_col,  1, 
1  '"Density  of  fluid  (kg/m3)',  rho,  vals) 


call  getr ( data_mode ,  ask,  inp_un,  inp_col,  outrun,  out_col,  1, 
1  '"Sound  speed  in  fluid  (m/s))',  snd_spd,  vals) 


call  getr (data_mode,  ask,  inp_un,  inp_col,  out_un,  out_col,  1, 
1  '"CWC  inner  radius  (m)',  r_inner,  vals) 

call  getr (data_mode,  ask,  inp_un,  inp_col,  out_un,  out_col,  1, 
1  '"CWC  outer  radius  (m) ' ,  r_outer,  vals) 

call  getr (data_mode,  ask,  inp_un,  inp_col,  out_un,  out_col,  1, 
1  '"Width  of  grounded  band  (m) ' ,  gnd_bnd_wid,  vals) 

call  getr (data_mode,  ask,  inp_un,  inp_col,  outrun,  out_col,  1, 
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1  ' "Width  of  active  band  (m)',  act_bnd_wid,  vals) 

call  geti ( data_mode ,  ask,  inp_un,  inp_col,  out_un,  out_col,  1, 
1  ' "Number  of  active  bands',  num_bnds,  vals) 

call  getchr (data_mode,  ask,  inp_un,  inp_col,  outrun,  out_col, 

1  '"File  containing  driven  bands',  bnd_file) 
if  (bnd_file(l:l)  .  eq.  bnd_file  =  run_name  //  bnd_file 


call 

1 

geti (data_mode,  ask,  inp_un,  inp_col,  out_un,  out_col,  1, 
'"Number  of  points  for  transfer  matrix',  num_j?nts,  vals) 

call 

1 

getr (data_mode,  ask,  inp_un,  inp_col,  out_un,  out_col,  1, 

' M Minimum  z  for  matrix  points  (m) ' ,  min_zp,  vals) 

call 

1 

getr ( data_mode ,  ask,  inp_un,  inp_col,  outrun,  out_col,  1, 
'"Maximum  z  for  matrix  points  (m)',  max_zp,  vals) 

call 

1 

getr ( data_mode ,  ask,  inp_un,  inp_col,  out_un,  out_col,  1, 
'"%  average  added  to  matrix  diagonals',  alpha,  vals) 

call 

1 

1 

getr ( data_mode ,  ask,  inp_un,  inp_col,  out_un, 
out_col,  max_spds, 

' "Minimum  z  for  controlled  region  (m) ' ,  min_zc,  vals) 

call 

1 

1 

getr (data_mode,  ask,  inp_un,  inp_col,  out_un, 
out_col ,  max_spds , 

'"Maximum  z  for  controlled  region  (m) ' ,  max_zc,  vals) 

call 

1 

getr ( data_mode ,  ask,  inp_un,  inp_col,  out_un,  out_col,  1, 
'"x  coordinate  for  map  (m)  '  ,  x__map,  vals) 

call 

1 

getr (data_raode,  ask,  inp_un,  inpcol,  out_un,  out_col,  1, 
'"y  coordinate  for  map  (m)',  y_map,  vals) 

call 

1 

getr ( data_mode ,  ask,  inp_un,  inp_col,  out_un,  out_col,  1, 
'"Minimum  z  to  map  (m) ' ,  z_map_min,  vals) 

call 

1 

getr (data_mode,  ask,  inp_un,  inp_col,  outrun,  out_col,  1, 
'"Delta  z  for  mapping',  z_map_delta,  vals) 

call 

1 

getr (data_mode,  ask,  inp_un,  inp_col,  outrun,  out_col,  1, 
'"Maximum  z  to  map  (m) ' ,  z_map_max,  vals) 

call 

1 

getchr (data_mode,  ask,  inp_un,  inp_col,  outrun,  out_col, 
'"File  containing  single  drive  vel.',  vel_file) 

z_bnds_max  =  (num_bnds  *  act_bnd_wid)  /  2. 
z_bnds_min  =  -  z_bnds_max 
z_arr_max  =  z_bnds_max  +  gnd_bnd_wid 
z_arr_min  =  -  z_arr_max 

wall_imp  =  cmplx(0.0,  -4.0e9  *  0.497  /  freq) 

return 

end 

CALCDRV.FOR  ( Calculates  band  drives) 
subroutine  calcdrv 

c  This  subroutine  calculates  the  band  drive  required  for  the  CWC. 

c  This  is  accomplished  by  constructing  the  set  of  simultaneous  eqns 
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c  M  X  =  Y  (1) 

c  —  - 

c  where  (for  each  equation)  the  values  for  the  elements  in  M  and  Y  are 
c  calculated  for  a  point  in  the  region  where  the  desired  sound  field  is 
c  to  be  created.  This  gives  m  equations  with  n  unknowns  (m  greater 
c  than  or  equal  to  n)  where  n  is  the  number  of  driven  bands  in  the  CWC. 

c  This  set  of  equations  is  then  collapsed  to  form  the  new  set  of  eqns: 

c  M'  X  =  Y'  (2) 

c  =  " 

c  consisting  of  n  equations  with  n  unknowns.  Equation  2  is  arrived  at 
c  by  multiplying  both  sides  of  Eq.  1  by  the  conjugate  transpose  of  M. 
c  Equation  2  is  then  solved  to  yield  the  required  drive  values  X. 
c  The  solution  to  Eq.  2  is  the  least-squares  solution  to  Eq.  1. 

c  L.  D.  Luker,  NRL/USRD,  25-OCT-1990.  Last  modified  on  17-JUL-1992. 

implicit  none 

include  ' CWCMDL. INC /list ' 

character  name* 30 

integer*4  i,  jf  bnd, 

1  uniq_col,  ! Number  of  unique  columns  in  M 

1  rows_used,  ! Total  rows  used  (some  zeros) 

1  perm_yec(max_bnds)  IPermutation  vector  for  dcxludc 

real*4  angle,  beta,  del_zp,  col_max,  mat_max 

real* 8  work_vec (max_bnds )  l Working  vector  for  dcxludc 

complex* 16 

1  M_mat(max_pnts,max_bnds) ,  ITransfer  matrix. 

1  M_star(max_bnds,max_pnts) ,  iConjugate  transpose  of  Mjmat 

1  M_prime(max_bnds,max_bnds) ,  ICollapsed  matrix. 

1  Y_mat ( max_pnts , max__spds ) ,  l Right-hand-side  of  the  equation. 

1  Y_jprime(max_bnds,max_spds) ,  ICollapsed  right-hand-side. 

1  det 

d  “  write ( out_un ,  ' ( / , ' '  <CALCDRV>  Entered . ' ' , / ) ' ) 

c  Set  up  CHIEF  and  generate  the  matrices 
if  (mod_type  .eq.  2)  then 
call  chzsub 
else 

call  chfsub 
endif 

numfpn  =  num_pnts 

del_zp  =  (max_zp  -  min_zp)  /  (num_pnts  -  1) 

c  Load  the  point  coordinates 
do  50  i  =  1,  nuin_pnts 
nfpxs(l,  i)  =0. 
nfpxs (2,  i)  =  0. 

nfpxs ( 3 ,  i)  =  min_zp  +  (i  -  1)  *  del_zp 
d  write (out_un,  '(,#  nfpxs(3,'',  i3,  ' ' )  =,f,  f8.4)') 

d  1  i,  nfpxs ( 3, i) 

50  continue 

if  (prnt_lvl  .eq.  4)  then 

open (unit  =  1,  file  =  r  CWCMDL . TMP ' , 

1  reel  =  out__len,  status  =  ' OLD 9 ) 

call  dcxin(M_mat,  name,  max_pnts,  num_pnts, 
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d  1  num_drv_bnds,  1,  out_len,  34  +  num_pnts) 

d  close (unit  =  1) 

d  goto  200 

d  endif 

c  Fill  the  M  matrix 

c  Determine  the  number  of  unique  columns 

if  (min_zp  .  eq.  -max_zp)  then  ICan  use  symmetry 
uniq_col  =  (num_drv_bnds  +  1)  /  2 
else  1 Can't  use  symmetry 

uniq_col  =  num_drv_bnds 
endif 

c  For  each  unique  column 

matjmax  =  0. 
do  100  i  =  1,  uniq_col 

call  filled (i,  M_mat,  col_max) 
mat_max  =  max  (matjmax,  col_max) 

100  continue 

c  Fill  any  remaining  columns 

do  150  i  =  uniq_col  +  1,  num_drv_bnds 
do  150  j  =  1,  num_pnts 

M_mat(j,i)  =  M_mat( num_pnts  +  1  -  j,  num_drv_bnds  +  1  -  i) 
150  continue 

c  Write  the  transposed  M  matrix  as  real  and  imaginary 
d  open(unit  —  1,  file  *  run_name  //  '.RE', 

d  1  reel  =  out_len,  status  =  'NEW') 

d  do  160  i  =  1,  num_drv_bnds 

d  write ( 1 , ' ( 85  (x, lpell . 4 ) ) ' ) 

d  1  (dreal(M_mat( j , i) ) ,  j = 1 , num_pnts ) 

160  continue 

d  close ( unit=l ) 

d  open (unit  =  1,  file  =  run_name  //  '.IM', 

d  1  reel  =  out_len,  status  =  'NEW') 

d  do  170  i  =  1,  num_drv_bnds 

d  write ( 1 , ' ( 85 ( x, lpell . 4 ) ) ' ) 

d  1  ( dimag (M_mat ( j , i) ) ,  j=l  ,num_jDnts ) 

170  continue 

d  close(unit=l ) 

c  Calculate  the  right-hand  sides 
do  200  j  =  1,  num_spds 

if  ((max_zc(j)  -  min_zc(j))  .eq.  0.)  then 
min_zc ( j )  =  min_zc ( j  -  1 ) 
max_zc(j)  =  max_zc(j  -  1) 
endif 

do  200  i  =  1,  nuinjpnts 

if  ((nfpxs(3,i)  .ge.  (min_zc(j)  -  .0001)) 

1  .and.  (nfpxs(3,i)  .le.  (max_zc(j)  +  .0001)))  then 

angle  =  ((2.  *  pi  *  freq)  /  phas_spd(j))  *  nfpxs(3,  i) 
d  type  * ,  i,  j,  phas_spd(j),  nfpxs(3,i),  angle 

Y_mat(i,j)  =  cmplx ( cos ( angle) , sin ( angle ) ) 
else 

Y_mat(i,j)  =  cmplx(0.0,0.0) 
endif 

200  continue 

rows_used  =  num  pnts 

c  Collapse  the  equations  for  the  least-square  solution, 
call  dcxlstsq(M_mat,  max_pnts, 

1  Y_mat,  max__pnts# 

1  Mjjrime,  max_bnds, 

1  Y_prime/  max_bnds. 
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1  M_star,  max_bnds, 

1  rows_used,  num_drv_bnds,  num_spds) 

*  l  num_j?nts,  num_drv__bnds ,  num_spds) 

d  call  dcxout ( M_mat ,  '<CALCDRV>  M  '//name, 

d  1  maxjnts,  rows_used,  num_drv_bnds ,  out_un,  out_len) 

*  1  max_pnts,  num_pnts,  num_drv_bnds ,  out_un,  out_len) 

d  call  dcxout ( Y_raat ,  '<CALCDRV>  Y', 

d  1  max_pnts,  rows_used,  num_spds,  out_un,  out_len) 

*  1  max__pnts,  num__pnts,  num_spds,  out_un,  out_len) 

d  call  dcxout ( M__prime ,  '<CALCDRV>  M'", 

d  1  max_bnds,  num_drv_bnds ,  num_drv_bnds ,  out_un,  out_len) 

c  Calculate  beta  as  alpha  percent  of  the  average  | M ' |  value, 
if  (alpha  .ne.  0.0)  then 
beta  =  0.0 

do  300  i  =  1,  num_drv_bnds 
do  300  j  =  1/  num_drv_bnds 

beta  =  beta  +  abs(M_prime( j,i) ) 

300  continue 

beta  =  alpha  *  beta  /  (100.  *  num_drv_bnds  **  2) 

c  Add  beta  to  the  diagonal  elements  of  M'. 
do  400  i  =  1,  num_drv_bnds 

M_prime(i,i)  *=  M_prime(  i,  i)  +  beta 
400  continue 

write (out_un,  '(/,''  "beta  =  ' ' , lpgll . 4/ ) ' )  beta 
endif 

d  call  dcxouttMjjrime,  '<CALCDRV>  M' '  +  beta' , 

d  1  max__bnds,  num_drv_bnds ,  num_jirv_bnds ,  outrun,  out_len) 

d  call  dcxout ( Y_prime ,  '<CALCDRV>  Y' ' ' , 

d  1  max_bnds,  num__drv_bnds ,  num_spds,  outrun ,  out_len) 

c  Calculate  the  drives  using  LU  decomposition 
call  dcxludc(M_j?rime,  max_bnds, 

1  num_drv_bnds,  work_vec,  perm_vec,  det) 

c  For  each  phase  speed 

do  600  j  =  1,  nura_spds 

c  Load  the  RHS  vector 

do  500  i  =  1,  num_drv_bnds 
drive ( i )  =  Y_prime ( i / j ) 

500  continue 

c  Solve  for  the  drives  with  back-substitution 

call  dcxlubk(M_prime,  max_bnds, 

1  num_drv_bnds ,  perm_vec,  drive) 

d  call  dcxout (drive,  '<CALCDRV>  DRIVE ' , 

d  1  max_bnds,  num_drv_bnds ,  1,  out_un,  out_len) 

c  Save  the  drives 

call  savedrv(j) 

600  continue 

return 

end 
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LOADBNDS.FOR  (Loads  band  drives  from  a  file) 

subroutine  loadbnds 

c  This  routine  loads  the  active  band  numbers  from  the  file  specified 
c  by  bnd_file. 

c  L.  D .  Luker,  NRL/USRD,  09-JAN-1991.  Last  modified  on  10-APR-1991. 
implicit  none 
include  ' CWCMDL. INC /list 1 
integer* 4  i 

d  type  *,  ' <LOADBNDS>  bnd_file  =  ',  bnd_file 

if  ( (bnd_file( 1 : 8 )  .eq.  'allbands')  .  or. 

1  (bnd_file{ 1:8)  .eq.  'ALLBANDS'))  then 

do  100  i  =  1,  num_bnds 
drv_bnd ( i )  =  i 

d  type  *,  '  drv_bnd(',  i,  ')  =  drv_bnd(i) 

100  continue 

else 

open(unit  =  bnd_un,  file  =  bnd_file,  status  =  'OLD') 
do  200  i  =  1,  max_bnds 

read(bnd_un,  *,  end  =  999)  drv_bnd(i) 
d  type  *,  '  dr  v_bnd ( ' ,  i,  ')  ■  drv_bnd(i) 

200  continue 

end  if 

999  num_drv_bnds  =  i  -  1 

d  type  *,  ' <LOADBNDS>  num_drv_bnds  =' ,  num_drv_bnds 

close  (unit  =  bnd_un) 
return 
end 


MAPFLDS.FOR  (Computes  predicted  pressure  fields) 

subroutine  mapflds 

c  This  routine  maps  the  pressure  fields  inside  the  CWC. 

c  L.  D*  Luker,  NRL/USRD,  27-SEP-1990.  Last  modified  on  29-MAY-1992. 

implicit  none 

include  ' CWCMDL. INC /list ' 

character*20  filename 
integer* 4  i,  j 

real*4  angle,  mag,  phase,  phasl80,  map_spd 

c  Set  up  CHIEF  if  necessary 

if  ((oper_mode  .eq.  2)  .and.  (mod_type  .eq.  2))  then 
call  chzsub 

else  if  ((oper_mode  .eq.  2)  .or.  (mod_type  .eq.  2))  then 
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call  chfsub 
endif 

c  Determine  the  number  of  points  to  map 

numfpn  =  nint(l  +  (z_map_max  -  z_map_min)  /  z_map_delta) 

c  Load  the  map  coordinates 
do  100  i  =  1/  numfpn 
nfpxs(l,  i)  =  x_map 
nfpxs(2,  i)  =  y_map 

nfpxs(3,  i)  =  z_map_min  +  (i  -  1)  *  z_map_delta 

100  continue 


c  For  each  phase  speed  to  map 

do  300  j  =  1 ,  num_map_spds 

c  Load  the  drives 

call  loaddrv(map_spd_num( j ) ,  map_spd) 

c  Calculate  the  nearfield  pressures  with  CHIEF 

if  ((oper_mode  .eq.  2)  .and.  (mod_type  .eq.  2))  then 

call  chznfp(drv_bnd( 1 ) ) 
else 

call  chfnfp(pmatx,  vmatx,  mdsize) 
endif 

c  Write  the  map  data 

filename  =  map_file ( 1 : index (map_file/ '  ')-l)  // 

1  '.MA'  //  char ( j  +  48)  //  char(O) 

open(unit  *  map_un,  file  =  filename,  status  =  'NEW') 

write ( map_un, f (' '  "Pressure-Field  Map  for'', 


200 


300 


f 7 . 2 , ' '  m/s  Phase  Speed.  Run  name:  '',a,i,//, 
99  "  z  (cm)  amp  pha  (deg)'', 

''  amp  err  (dB)  pha  err  (deg)'')') 
map_spd ,  run^name ,  j 
do  200  i  =  1,  numfpn 

angle  =  ((360.  *  freq)  /  map_spd)  *  nfpxs(3,  i) 
write (map_un, ' ( f 8 .2 , tl3, lpgl3 . 6 , t27 , Opf 7 .2, 

1  t39 , Opf 7 . 2 , t54 , Opf 7 . 2 ) ' ) 

1  nfpxs (3,  i)  *  100.,  mag(nfp(i) ) ,  phasl80( phase (nfp( i) )) , 

1  20.  *  alogl0(raag(nfp(i) ) ) ,  phasl80 (phase(nfp( i) )  -  angle) 

continue 

close  (unit  =  map_un) 
continue 
return 
end 
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MODIFICATIONS  TO  CHIEF  FOR  THE  CWC 


Description  of  the  Problem 

The  CHIEF  program,  as  currently  distributed,  can  give  erroneous  results  when  modeling  thin 
radiating  bodies.  Although  the  problem  with  thin  bodies  is  not  mentioned  in  the  CHIEF  user’s  manual 
it  has  been  observed  previously  where  it  was  shown  that  the  results  do  improve  if  a  finer  division  of  the 
surface  is  utilized  [79].  Unfortunately,  dividing  the  surface  into  very  small  subdivisions  would  increase 
the  size  of  the  CHIEF  model  of  the  CWC  to  the  point  that  it  could  not  be  run  on  the  available 
computational  facilities.  The  problem  with  thin  radiating  bodies  seems  to  arise  from  errors  introduced 
when  evaluating  the  surface  integrals  for  points  that  are  very  near  the  integration  surface.  In  particular, 
errors  are  introduced  for  thin  bodies  when  evaluating  the  surface  integral  for  a  point  that  is  on  the  "back" 
side  of  the  integration  surface.  The  limiting  case  of  this  integration  problem  (i.e.,  infinitely-thin  body) 
is  the  same  problem  encountered,  and  solved,  by  CHIEF  when  integrating  the  Helmholtz  self  integrals. 

Solution  Implemented  for  the  CWC 

One  solution  to  the  problem  of  thin  bodies  is  to  evaluate  the  troublesome  integrals  using  the  same 
technique  that  CHIEF  uses  to  evaluate  the  Helmholtz  self  integrals  [80],  This  consists  of  dividing  the 
surface  to  be  integrated  into  four  triangular  subregions  (by  cutting  the  original  rectangular  region  with 
diagonals)  and  performing  a  change  of  variables.  This  allows  integration  of  a  function  that  grows  rapidly 
as  the  center  of  the  original  region  is  approached  (such  as  a  singularity  at  the  center).  For  this  technique 
to  work,  the  front  and  back  of  the  thin  body  must  both  be  divided  into  the  same  arrangement  of 
subregions.  To  put  it  more  precisely,  the  front  and  back  surfaces  must  be  divided  in  such  a  way  that  (in 
u,v  coordinates)  the  center  of  the  back  subregions  approach  the  center  of  the  front  subregions  as  the 
thickness  of  the  thin  body  approaches  zero.  Fortunately,  this  condition  is  satisfied  for  the  CWC  models. 

This  technique  has  been  implemented  in  CHIEF  by  modifying  two  subroutines.  The  subroutine 
SURMAT  has  been  modified  to  call  a  new  integration  routine  (CLSINT)  for  each  surface  area  on  the  first 
symmetry  block.  The  subroutine  CLSINT  is  a  combination  of  the  original  subroutines  MATGEN  and 
SLFINT.  This  routine  performs  all  the  surface  integrals  for  the  first  symmetry  block  using  the  algorithm 
that  was  formerly  only  used  for  the  self  integrals. 

The  new  technique  was  tested  by  using  CHIEF’S  point-source  check  for  a  model  of  the  CWC.  The 
CWC  was  modeled  as  a  cylindrical  tube  25.2-cm  long  with  an  inner  radius  of  1.55  cm  and  a  wall 
thickness  of  0.75  mm.  The  tube  was  modelled  using  18-fold  rotational  symmetry.  The  inner  and  outer 
surfaces  of  the  tube  were  divided  into  bands  that  were  6.0-mm  wide.  A  single  point  source  was  centered 
on  the  axis  of  the  tube  (i.e.,  at  z  =  0).  The  point-source  check  was  run  at  800  Hz  using  CHIEF  with  and 
without  the  modifications.  Figure  Cl  shows  the  percent  errors  of  the  surface  pressures  computed  by 
CHIEF  relative  to  the  pressures  from  the  point-source  model  for  the  outer  band  surfaces.  The  results  show 
that  the  errors  are  reduced  dramatically  by  using  the  modified  version  of  CHIEF. 
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Fig.  Cl  —  Pressure  errors  for  the  outer  band  surfaces  using  the  point-source  check 
from  CHIEF  with  and  without  the  modifications. 


The  disadvantages  of  this  technique  include  the  requirement  placed  on  the  layout  of  the  surface 
subregions  mentioned  above  and  the  fact  that  the  more  complex  integration  routine  takes  longer  to  run. 
The  additional  time  required  could  be  reduced  by  only  using  the  more  complex  integration  routine  for  the 
self  integrals  and  the  point  that  is  opposite  the  area  being  integrated,  but  this  would  require  considerable 
knowledge  about  how  the  user  had  subdivided  and  ordered  the  surface  areas. 
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FORTRAN  Code 
SURMAT.FOR 


SUBROUTINE  SURMAT ( FREQ , SYMTYP , NBLKS , PMATX , VMATX , MTXDIM ) 
C  GENERATION  OF  P  AND  V  SURFACE  MATRICES 


* 


PARAMETER 

PARAMETER 

PARAMETER 

PARAMETER 

PARAMETER 

PARAMETER 

PARAMETER 

PARAMETER 


(MXSREG=500) 

(MXIPS=20 ) 

(MXARS=500) 

(MXGAUS=64) 

(MXQPTS=512) 

(MXBLKS=100) 

(MXFFP=361) 

(MXNFP=361) 


PARAMETER  (MXFPS=520) 

PARAMETER  (MXFPS=MAXO (MXARS+MXIPS ,MXFFP,MXNFP ) ) 
*** 


COMPLEX  PMATX (MTXDIM ) ,  VMATX ( MTXDIM ) 

CHARACTER* 3  SYMTYP 

COMMON/CONST/RHO, C 
COMMON/PRTCOM/NUNPRT, NUNERR 

COMMON/NDASG/NDQPTS , NDPMXS , NDVMXS , NDDECM, NDVELS , NDSPS , 

*  NDPMXF , NDVMXF , NDPMXN , NDVMXN , NDPSSP , NDEXPR , NDCOMV , 

*  NDTEMP , NDZRDB , NDPATB 
COMMON/SVALS/NSREG , NSEQNS ( MXSREG ) , SUL ( MXSREG ) , SUU ( MXSREG ) , 

*  SVL ( MXSREG ) , S VU ( MXSREG ) , NSU ( MXSREG ) , NSV ( MXSREG ) , 

*  CCS ( 1 0 , MXSREG ) , TRNSS ( 3 , MXSREG ) , I Z AX ( MXSREG ) , 

*  IORDU ( MXSREG ) , IORDV ( MXSREG ) , NCCEQS 
COMMON / IPTS /NUMIPS , IPXS ( 3 , MXIPS ) 

REAL  IPXS 


COMMON /QPTS /QPS ( 7 , MXQPTS ) 

COMMON/GAUSNS/WKGAUS (MXGAUS ) , ROOTU(MXGAUS ) , ROOTV(MXGAUS ) , 

*  WTUV ( MXGAUS+1 , 2 ) 

COMMON / FLDPTS / FPBK 1 ( 4 , MXFPS ) , FPC ( 4 , MXFPS ) 

COMMON /WORK /VEC1 ( MXARS ) , VEC2 ( MXBLKS ) , VEC3 ( MXBLKS ) , VEC4 ( MXBLKS ) , 

*  VEC5 (MXFPS ) ,VEC 6 (MXFPS) 

COMPLEX  VEC2, VEC3,VEC4,VEC5,VEC6 
DIMENSION  FACTRS( MXARS) 

EQUIVALENCE  ( FACTRS ( 1 ) , VEC1 ( 1 ) ) 

COMMON /PRGVLS/NDIMPV, NUMARS , NUMSFP , NUMFFP , NUMNFP, NWDVEC 
COMMON / SURARS /AREAS ( MXARS ) 

COMMON /TAPREC/RECRD( 10 ) , IRECRD( 30 ) 

COMMON /TAPRC1/ARECRD( 10 ) 

***************  impedance  coating  modification  ************* 

COMMON  /  COAT  ING  /  2  COAT ,  NUMZ 
COMPLEX  Z COAT (MXARS) 

************************************************************ 


CHARACTER* 4  ARECRD 
COMPLEX  JRHOWP 


SURMAT  ERROR  CHECKS 
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IERR=0 

IF  ( NSREG . GT . MXSREG . OR . NSREG . LE . 0 )  THEN 
IERR=2 

WRITE ( NUNERR ,1001) NSREG 

1001  FORMAT ( / , '  ERROR  IN  SURMAT  -  NSREG=',I6) 

END  IF 

IF  ( NBLKS . GT . MXBLKS . OR . NBLKS . LE . 0 )  THEN 
IERR=1 

WRITE ( NUNERR ,1002) NBLKS 

1002  FORMAT ( / , '  ERROR  IN  SURMAT  -  NBLKS=',I6) 

END  IF 

IF  ( SYMTYP . NE . ' REF ' .AND . SYMTYP . NE. ' ROT ' )  THEN 
IERR=1 

WRITE (NUNERR, 1003) SYMTYP 

1003  FORMAT ( / , '  ERROR  IN  SURMAT  -  SYMTYP=  ',A3) 

END  IF 

IF  (NUMIPS.GT.MXIPS.OR.NUMIPS.LT. 0)  THEN 
IERR=1 

WRITE ( NUNERR ,1004) NUMIPS 

1004  FORMAT ( / , '  ERROR  IN  SURMAT  -  NUMIPS=',I6) 

END  IF 

IF  (IERR.EQ.2)  STOP 
NAR=0 

DO  1110  1=1, NSREG 
NAR=NAR+N  SU ( I ) *NSV ( I ) 

IF  ( NSEQNS ( I ) • GT . NCCEQS . OR . NSEQNS ( I ) . LE . 0 )  THEN 
IERR=1 

WRITE ( NUNERR , 1006)1, NSEQNS(I) 

1006  FORMAT ( / , '  ERROR  IN  SURMAT  -  REG#,NSEQNS= ' , 216 ) 

END  IF 

IZ ABS=IABS ( IZ AX (I) ) 

IF  ( IZABS.GT. 3 .OR. IZABS .EQ. 0)  THEN 
IERR=1 

WRITE (NUNERR, 1007) I, IZAX( I) 

1007  FORMAT ( / , '  ERROR  IN  SURMAT  -  REG#, IZAX=' , 216 ) 

END  IF 

IF  ( NSU ( I ) *NSV ( I ) . LE . 0 )  THEN 
IERR=1 

WRITE (NUNERR, 1008 ) I ,NSU ( I ) ,NSV( I ) 

1008  FORMAT ( / , '  ERROR  IN  SURMAT  -  REG#,NSU,NSV=' , 316 ) 

END  IF 

IORDRU=IORDU ( I ) 

IORDRV=IORDV( I ) 

IF  ( IORDRU . GT . MXGAUS . OR . IORDRV . GT . MXGAUS . OR . IORDRU* IORDRV 
*  . GT.MXQPTS. OR. IORDRU * IORDRV. LE.0)  THEN 

IERR=1 

WRITE (NUNERR, 1009 ) I , IORDRU, IORDRV 

1009  FORMAT ( / , '  ERROR  IN  SURMAT  -  REG#, IORDRU, IORDRV= ', 316 ) 
END  IF 

1110  CONTINUE 

IF  ( NAR . GT . MXARS . OR . NAR . LE . 0 )  THEN 
IERR=1 

WRITE ( NUNERR ,1111) NAR 

1111  FORMAT ( / , '  ERROR  IN  SURMAT  -  #OF  AREAS=',I10) 

END  IF 

NFP=NAR+NUMIPS 

IF  ( NFP . GT . MXFPS . OR . NFP . LE . 0 )  THEN 
IERR=1 

WRITE ( NUNERR , 1 1 1 2 ) NFP 

1112  FORMAT ( / , '  ERROR  IN  SURMAT  -  #OF  FIELD  POINTS= ' , 110 ) 

END  IF 

MTXT=NAR*NFP 
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IF  ( MTXT . GT . MTXDIM . OR . MTXT . LE . 0 )  THEN 
IERR=1 

WRITE ( NUNERR ,1113) MTXT 

1113  FORMAT ( / , '  ERROR  IN  SURMAT  -  ',110,'  WORDS  NEEDED  FOR  P,V' ) 

END  IF 

IF  (IERR.NE.0)  STOP 

PI=ACOS ( — 1 • 0 ) 

TWOPI=2 . 0*PI 
FOURPI=4 . 0*PI 
ROMEGA=TWOP I *  FREQ 

JRHOWP=CMPLX (0.0, ROMEGA*RHO ) /FOURPI 

SMALLK=ROMEGA/C 

NDIMPV=MTXDIM 

COMPUTE  QUADRATURE  POINTS  FOR  EACH  SURFACE  SUBDIVISION (AREA) 
STORE  ON  FILE/UNIT  NDQPTS 

COMPUTE  FIELD  POINT  FOR  EACH  SURFACE  AREA  AND  INTERIOR  POINT 
STORE  IN  ARRAY  FPBK1 


CALL  IOSUB ( NDQPTS, 10, QPS,0) 

ISUB=0 

DO  10  1=1 , NSREG 

CALL  GAUSS ( IORDU ( I ) , ROOTU , WTUV (1,1), WKGAUS ) 

WTUV ( MXGAUS+1 , 1 ) = IORDU ( I ) 

CALL  GAUSS ( IORDV( I ) , ROOTV , WTUV (1,2), WKGAUS ) 

WTUV ( MXGAUS+1 , 2 ) =IORDV ( I ) 

DELSU= ( SUU ( I ) — SUL( I ) ) /NSU ( I ) 

DELSV= ( SVU ( I ) -SVL ( I ) ) /NSV ( I ) 

HDELSU=0.5*DELSU 
HDELSV=0 . 5*DELSV 
FACTI=HDELSU*HDELSV 
UAVE=HDELSU+SUL ( I ) 

DO  8  J=1 , NSU ( I ) 

VAVE=HDELSV+SVL ( I ) 

DO  6  K=1 , NSV ( I ) 

ISUB=ISUB+1 
FACTRS ( I SUB ) =FACTI 

CALL  QPTGEN ( UAVE , VAVE , HDELSU , HDELSV, NSEQNS ( I) , CCS (1,1), 

*  TRNSS (1,1), IZ AX ( I ) , IORDU ( I ) , IORDV ( I ) , ROOTU , ROOTV , QPS ) 

CALL  IOSUB (NDQPTS, 1 ,WTUV, 2* (MXGAUS+1 ) ) 

CALL  IOSUB ( NDQPTS , 1 , QPS , 7  * IORDU ( I ) * IORDV ( I ) ) 

CALL  CCUNMD ( UAVE, VAVE, NSEQNS ( I ) , CCS (1,1), TRNSS (1,1), IZAX( I ) , 

*  FPBK1 ( 1 , ISUB ) , FPBK1 ( 2 , ISUB ) , FPBK1 ( 3 , ISUB ) , D1 , D2 , D3 , DD J , 0 ) 

6  VAVE=VAVE+DELSV 

8  UAVE=UAVE+DELSU 

10  CONTINUE 

NUMARS=ISUB 

NUMFPS=ISUB 

NWDVEC=2*NUMARS 

IF ( NUMIPS . NE . 0 )  THEN 
DO  20  J=l, NUMIPS 
ISUB=ISUB+1 
DO  15  1=1,3 
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1 5  FPBK1 ( I , ISUB ) =IPXS ( I ,  J ) 

20  CONTINUE 

END  IF 

NUMFPS=NUMFPS+NUMIPS 
NUMSFP=NUMFPS 
MTXTOT=NUKARS *  NUMFPS 

GENERATION  OF  SURFACE  MATRICES 
P  MATRICES  TO  FILE/UNIT  NDPMXS 

V  MATRICES  TO  FILE/UNIT  NDVMXS 

INTTYP=1 

CALL  IOSUB ( NDPMXS ,10,  PMATX , 0 ) 

CALL  IOSUB ( NDVMXS ,10, VMATX , 0 ) 

DO  50  IBLK=1 , NBLKS 

C  ***************  SYMMETRY  CONSIDERATIONS  ******************* 

IF ( SYMTYP .  EQ .  ' REF ' )  THEN 

CALL  FPREFC ( IBLK  ,  NUMFPS , INTTYP , FPBK1 , FPC ) 

END  IF 

IF  ( SYMTYP . EQ . ' ROT ' )  THEN 
GAMMA=TWOPI /NBLKS 

CALL  FPROTC ( IBLK , GAMMA, NUMFPS , INTTYP , FPBK1 , FPC ) 

END  IF 

Q  *********************************************************** 

C ============ =======  Modified  for  CWC  (thin  wall)  ============= 

IF ( IBLK . EQ . 1 )  THEN 

CALL  CLSINT ( PMATX , VMATX , NUMARS , NUMFPS , SMALLK , INTTYP , 

*  ROOTU , ROOTV , WTU V , WKGAUS , MXG AU S , 

*  NSREG , NSEQNS , CCS , TRNSS , IZ AX , SUL , SUU , SVL , SVU , 

*  NSU,NSV,NDQPTS,QPS, FPC, AREAS) 

ELSE 

CALL  MATGEN ( NUMARS , NUMFPS , FACTRS , INTTYP , WTUV , MXGAUS , 

*  NDQPTS,QPS, FPC, SMALLK, PMATX, VMATX) 

END  IF 

DO  30  I=l,MTXTOT 

D  IF  ( IBLK.EQ. 1 )  PRINT  *,  'I,P,V',  I,  PMATX(I),  VMATX(I) 

VMATX ( I ) = VMATX ( I ) * JRHOWP 
30  PMATX ( I ) =-PMATX ( I ) /FOURPI 

IF ( IBLK.EQ. 1)  THEN 
DO  40  1=1, NUMARS 
K=( 1-1 ) *NUMFPS+I 
40  PMATX (K)=PMATX(K) +0.5 

END  IF 

C***************  impedance  coating  modification  ************* 

CALL  MTXSWZS ( PMATX , VMATX , NUMARS , NUMZ , NUMIPS , ZCOAT ) 
c************************************************************ 


CALL  IOSUB ( NDPMXS , 1 , PMATX , 2  *MTXTOT ) 
CALL  IOSUB ( NDVMXS , 1 , VMATX , 2  *MTXTOT ) 

50  CONTINUE 

RE CRD ( 1 ) =FREQ 
RE CRD ( 2 ) =RHO 


114 


NRL  MEMORANDUM  REPORT  7385 


RE CRD ( 3 ) =C 
IRECRD(1)=NBLKS 
IRECRD ( 2 ) =MTXDIM 
I RE CRD ( 3 ) =N UMARS 
IRECRD ( 4 ) =NUMSFP 
ARECRD ( 1 ) =SYMTYP 

RETURN 

END 


CLSINT.FOR 


SUBROUTINE  CLSINT ( PMATX , VMATX , NUMARS , NUMFPS , SMALLK , INTTYP , 

*  ROOTU , ROOTV , WTUV , WKGAUS , MXGAUS , 

*  NSREG , NSEQNS ,  CCS, TRNSS, I Z  AX , SUL f  SUU , SVL , SVU , 

*  NSU , NSV, NDQPTS, QPS , FPC , AREAS ) 

C  COMPUTES  P  AND  V  MATRIX  ELEMENTS  FOR  CLOSE  SURFACES  AND  SELF  INTEGRALS 


COMPLEX  PMATX ( NUMFPS , NUMARS ) ,  VMATX ( NUMFPS , NUMARS ) 

DIMENSION  ROOTU ( MXGAUS )  , ROOTV (MXGAUS ), WTUV (MXGAUS+1 ,2) 

DIMENSION  WKGAUS (MXGAUS) 

DIMENSION  AREAS ( NUMARS ) , QPS ( 7 , * ) , FPC ( 4 , NUMFPS ) 

DIMENSION  NSEQNS ( *  )  , SUL ( * ) , SUU ( * ) , SVL( * ) , SVU( * ) ,NSU( * ) ,NSV( * ) , 

*  CCS (10,*), TRNSS (3,*) ,IZAX(*) 

COMPLEX  PI I,  VII 

CALL  IOSUB ( NDQPTS, 1 0, QPS, 0) 

DO  100  N=l, NUMARS 

CALL  IOSUB (NDQPTS, 2, WTUV, 2* (MXGAUS+1) ) 

IORDRU=WTUV( MXGAUS+1 , 1 ) 

IORDRV=WTUV( MXGAUS+1 , 2 ) 

CALL  IOSUB ( NDQPTS , 2 , QPS , 7  * IORDRU* IORDRV ) 

11=0 

DO  100  1=1, NSREG 

CALL  GAUSS ( IORDRU , ROOTU , WTUV( 1,1), WKGAUS ) 

CALL  GAUSS ( IORDRV, ROOTV , WTUV (1,2) , WKGAUS) 

DELSU= ( SUU ( I ) -SUL ( I ) ) /NSU ( I ) 

DELSV= ( SVU ( I ) -SVL ( I ) ) /NSV(I) 

HDELSU=0 • 5*DELSU 
HDELSV=0.5*DELSV 
FACTI=HDELSU*HDELSV 
UAVE=HDELSU+SUL ( I ) 

DO  80  J=1 , NSU ( I ) 

VAVE=HDELSV+SVL ( I ) 

DO  60  K=1 ,NSV ( I ) 

11=11+1 

CALL  HELMSI (UAVE, VAVE , HDELSU, HDELSV, NSEQNS ( I ) , CCS ( 1 , I ) , 
TRNSS (1,1), IZAX ( I ) , 

INTTYP , SMALLK , FACTI , IORDRU , IORDRV, FPC ( 1 , N ) , 

ROOTU , ROOTV , WTUV (1,1), WTUV ( 1 , 2 ) , PI I , VI I , AREAVL ) 
PMATX ( N, II ) =PII 
VMATX ( N , I I ) =VI I 
AREAS (II) = AREAVL 
VAVE=VAVE+DELSV 
UAVE=UAVE+DELSU 
CONTINUE 
RETURN 
END 
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